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Abstract: The goal of this work, motivated by the desire to understand causal- 
ity in classical and quantum gravity, is an in depth investigation of causality in 
classical field theories with quasilinear equations of motion, of which General 
Relativity is a prominent example. Several modern geometric tools (jet bundle 
formulation of partial differential equations (PDEs), the theory of symmetric hy- 
. perbolic PDE systems, covariant constructions of symplectic and Poisson struc- 

tures) and applies them to the construction of the phase space and the algebra 
of observables of quasilinear classical field theories. This construction is shown 
to be diffeomorphism covariant (using auxiliary background fields if necessary) 
using categorical tools in a strong parallel with the locally covariant field theory 
(LCFT) formulation of quantum field theory (QFT) on curved spacetimes. In 
this context, generalized versions of LCFT axioms become theorems of classical 
field theory, which includes a generalized Causality property. Considering defor- 
mation quantization as the connection to QFT, a plausible conjecture is made 
about the Causal structure of quantum gravity. In the process, conal manifolds 
are identified as the generalization of the causal structure of Lorentzian geome- 
try to quasilinear PDEs. Several important concepts and results are generalized 
from Lorentzian to conal geometry. Also, the proof of compatibility of the Peierls 
formula for Poisson brackets and the covariant phase space symplectic structure 
for hyperbolic systems is generalized to now encompass systems with constraints 
and gauge invariance. 



1. Introduction 



Great strides have been made in the past decade in terms of the definition and 
perturbative construction quantum field theories (QFTs) in the algebraic frame- 
work [3] . Notable successes include the extension of perturbative renormalization 
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of non- linear field theories to arbitrary, curved (though non-dynamical), glob- 
ally hyperbolic spacetimes and the treatment of gauge theories [64','47','95| . These 
ingredients are, for instance, sufficient to study perturbative general relativity 
(GR) in this framework. Unfortunately, deeply rooted in the formalism of QFT 
on curved spacetimes is the reliance on a non-dynamical background metric to 
supply the causal structure that defines the singularity and support structures 
of various n-point functions, time ordered products, and retarded products. 

The objects in classical field theory that are most closely related to these 
QFT objects are the retarded and advanced Green functions of the dynamical 
field equations of motion linearized at some particular solution (the dynamical 
linearization point). But we know that, in classical gravity, the singularity and 
support structures of these classical Green functions change as a function of the 
metric that provides the dynamical linearization point. In perturbative GR, the 
dynamical metric can be expressed as g = g + nh, where 5 is a fixed background 
metric and h is the dynamical perturbation, is proportional to Newton's grav- 
itational constant, which is to be treated as the formal perturbation parameter. 
As already mentioned, the singularity structure of the Green function of the 
linearized gravitational wave equation depends on the full, dynamical lineariza- 
tion point g = g + nh. But, to order 0{k'^), this singularity structure is fully 
determined by the non-dynamical background metric g. 

Clearly, the necessary dependence of this singularity structure on the dynam- 
ical perturbations h will appear at higher orders in k. However, at least in the 
case of QFT n-point functions, it is not a priori clear what form this modifica- 
tions would take or how they are to be recognized in the result of some higher 
order perturbative calculation. Given the tight connection between the singu- 
larity structure of Green functions and what is commonly referred to as causal 
structure, it is fair to say that the causal structure of gravity is dynamical. There- 
fore, one might say that in the current state of affairs one does not know how 
to recognize a dynamical causal structure in perturbative QFT nor even how 
to interpret the meaning of causal structure in a non-perturbative QFT, whose 
classical version has dynamical causal structure. 

Unfortunately, currently, QFT of GR is only accessible to us perturbatively. 
Eventually, though, we would like to move to a non-perturbative formulation. 
However, at the non-perturbative level, the axiomatic framework of locally co- 
variant field theory (LCFT), both classical and quantum, encompasses only the- 
ories whose equations of motion are systems of PDEs with principal symbols that 
only depend on a non-dynamical background metric. This explicitly excludes so- 
called quasilinear equations like GR, hydrodynamics, elasticity and non-linear 
<T-models. Quasilinear systems are defined by the property that their principal 
symbols (roughly, the coefficients of the highest order derivatives, which play a 
major role in the properties of the previously mentioned Green functions) de- 
pend on the dynamical fields. Fortunately, at the classical level, GR and other 
quasilinear systems can be studied non-perturbatively using classical PDF the- 
ory. Using these tools, the causal structure can be described and understood 
in fairly direct terms, closely related to the geometry of domain 0/ dependence 
theorems. 

This work has multiple goals. In broad outline, it is concerned with obtaining 
a deep, non-perturbative understanding of an algebraic formulation of causality 
in classical field theories with dynamical causal structures, in a way that parallels 
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the well understood geometric formulation of well-posedness in hyperbolic PDE 
systems. Once such an understanding has been obtained another aim is to formu- 
late a conjecture as to how it would translate to quantum field theory. Of course, 
at this point, one can only expect to make conjectures in this direction, owing to 
the difficulties in dealing with quantum field theories non-perturbatively. Once 
formulated, these notions could be adapted to a formal perturbative context, 
where they could be used in conjunction with practical calculations, for instance 
in perturbatively quantized GR. 

More specifically, another goal is to present in a unified way all the necessary 
geometric aspects of the construction of a classical field theory. These aspects 
include an intrinsic geometric formulation of PDE systems using jet bundles, an 
application of the theory of symmetric hyperbolic PDE systems to the construc- 
tion of the phase space, and the use of the covariant phase space method and 
Peierls formula to endow it with symplectic and Poisson structures. Each of these 
aspects is well known in some segments of the mathematical analysis and math- 
ematical physics communities, but certainly is not familiar to the broad high 
energy theory and theoretical relativity communities. Some of these aspects are 
presented here also in potentially novel ways, which highlight how they are used 
together. In a sense, this paper can be seen as a follow that supplements with 
symplectic geometric and Poisson algebraic aspects the older work of Geroch [5D] 
who essentially applied the theory of symmetric hyperbolic systems to a unified 
construction of the spaces of solutions of relativistic classical field theories. 

Finally, this paper also aims to apply the category-theoretic methods that 
have proven so useful in formulating the algebraic structure of QFT on curved 
spacetimes (now known as LCQFT) to both the algebraic and dual geometric 
formulation of classical field theory. Essentially, these tools are used to draw a 
strong parallel or duality between the notion of well-posedness from classical 
PDE theory and the algebraic properties of classical field theory, so that both of 
which could be formulated as mutually dual, succinct theorems. These theorems 
could later be used as a starting point for an improved axiomatization of classical 
and quantum field theories. 

1.1. Outline of the work. Sect. [5] recalls the axiomatic framework of LCFT and 
emphasizes that it currently excludes quasilinear systems like GR. 

Sect. 13] introduces the geometric formulation of PDE systems in terms of jet 
bundles. Sect. 13. 1[ with reference to Apdx. |^ for background on jets. It then 
focuses on the notion of a quasilinear, hyperbolic system of PDEs, with special 
emphasis on first order, symmetric hyperbolic systems. Sect. 13.21 The precise for- 
mulation of the condition of symmetric hyperbolicity given there differs slightly 
from the standard one in order to fit better in the subsequent geometric context. 
It concludes with a discussion of prolongation and equivalence of PDE systems, 
needed for the reduction of variational equations of motion to hyperbolic form, 
Sect. [321 

Sect. 131 examines the role of characteristics of a PDE system in its causal 
structure. Sect. 14.11 identifies chronal and spacelike cone bundles as essential 
abstractions of the resulting causal order. Sect. 14.21 groups solutions together 
into slow patches by the property of being slower than a certain fixed chronal 
or spacelike cone bundle. Sect. 14.31 uses the thus far developed notions of causal 
structure to state some standard theorems of PDE theory, which will be used as 
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the basic tools in the construction of classical field theories. Sect. 14.41 presents 
some more specialized results on linear systems, including global existence and 
properties of Green functions. The presented material is standard, but may be 
formulated in a slightly novel way to fit better with the geometric context of the 
Peierls formula for the Poisson bracket appearing later. 

Sect. [5] constitutes roughly half of the technical bulk of this work. It collects 
the mathematical tools presented thus far and applies them in a unified way to 
construct the phase space of a classical field theory and its algebra of observ- 
ables. Sect. 15.11 reviews the covariant phase space method for constructing the 
symplectic form of a variational PDE system. Since the focus of this work is more 
geometrical than functional analytical, we avoid most details of the treatment of 
infinite dimensional manifolds, of which the space of solutions of a PDE system, 
to be identified with the phase space of a classical field theory, is an example. 
The purpose of Sect. 15.21 is to introduce enough technical machinery to formally 
with the tangent and cotangent spaces of the space of solutions, while keeping 
with the spirit of the preceding remark. Sect. 15.31 identifies the space of solu- 
tions with the phase space, by endowing it with (formal) symplectic and Poisson 
structures. The Poisson structure is given by a generalized Peierls formula. This 
Peierls formula is shown to be equivalent to the Poisson bracket obtained from 
the covariant phase space symplectic form, a result that, compared to previous 
forms of this equivalence, now also encompasses hyperbolic PDE systems with 
constraints. As an application of this construction, a refined version of classical 
microcausality is proven both for algebras of observables on a given slow patch 
as well as on the global phase space. Sects. 15.3.51 and 15.3.61 Also of note is a 
discussion in Sect. 15.3.51 of the relation between the on-shell formalism used in 
this work and the off-shell formalism used in other related literature. 

Sect, ini makes a connection between the constructions thus far and the co- 
variant formalism of LCFT via the notion of natural bundles and PDE systems, 
relying on some background material on category theory from Apdx. [BJ For 
convenience. Sect. 16.21 systematically summarizes the notation for various con- 
cepts and constructions appearing in this work, while also remarking on their 
categorical and functorial properties. 

Sect. [7] constitutes the other half of the technical bulk of this work. It gen- 
eralizes the axioms of LCFT to quasilinear field theories taking into account 
their dynamical causal structures, and proves them as theorems of classical field 
theory. Unfortunately, at this point the proofs rely on a number of sufficient 
technical hypotheses. The necessity of these hypotheses remain to be examined 
on a case by case basis or in the context of deeper investigation of the relevant 
functional analytical details. In particular, the generalization of the Causality 
axiom provides a sought deeper understanding of causality in quasilinear clas- 
sical field theories and its algebraic formulation. An important check, Sect. 17.41 
is that these theorems reduce to the standard axioms of LCFT in the case of a 
well-posed semilinear classical field theory, as is to be expected. 

Sect, m recalls the notion of deformation quantization and remarks that it is 
the leading candidate for the modern formulation of what it means to quantize a 
classical mechanical system. This notion of quantization is then used to make a 
conjecture about the translation of the algebraic formulation of Causality from 
classical to quantum field theory. In particular, this conjecture answers an old 
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question about the structure of the commutator of two quantum metric field 
operators in the quantum field theory of GR. 

Finally, Sec. IHl concludes with a discussion of the results presented in this 
work and lists several ideas, conjectures and hypotheses that had arisen in the 
process that would be fruitful to investigate further. 

2. Free, Interacting and Perturbative Locally Covariant Field 
Theories 

The currently leading axiomatic framework in the algebraic framework for quan- 
tum field theory is called locally covariant field theory (LCFT) and rests on the 
axioms originally proposed by Brunetti, Fredenhagen and Verch |25j . These ax- 
ioms rest on simple physical principles, namely locality, causality, and the ex- 
istence of a dynamical law, and have a convenient mathematical formulation. 
Moreover, they are flexible in that they can and have been adapted to the clas- 
sical {locally covariant classical field theory abbreviated again LCFT), quantum 
{locally covariant quantum field theory or LCQFT), and both classical and quan- 
tum perturbative contexts (respectively abbreviated pLCFT or pLCQFT) [5]. 

2.1. Brunetti-Fredenhagen-Verch axioms for LCFT. To state the axioms in a 
succinct way, we need to appeal to some notions of category theory. Some basic 
information on categories and functors can be found in I10,9pij . A classical me- 
chanical system, at a bare minimum is described by a Poisson algebra over K, 
say {F, {}), called the algebra of observables. Actually, we are likely to be inter- 
ested only the algebras F that correspond to smooth functions on some phase 
space manifold, with the bracket {} defined by a Poisson tensor. However, as we 
will not need it, we do not consider a detailed characterization of these algebras. 
Such Poisson algebras form the category *Poiss, with Poisson homomorphisms 
as morphisms. On the other hand, it will be important below to restrict possible 
homomorphisms to injective ones. Thus we define the following categories. 

Definition 1. Let *Poiss be the category of Poisson algebras over M as objects 
and Poisson homomorphisms as morphisms. Let *PoisSj be the subcategory where 
the morphisms are restricted to injective homomorphisms. 

A field theory assigns a classical mechanical system to a region of spacetime, 
that is, a Poisson algebra of observables whose spacetime support^ overlap with 
that region. A field theory is locally covariant, if these assignments are coherent 
with respect to embedding smaller spacetimes into larger ones, in a sense to be 
specified below. First, however, we must define what constitutes a spacetime. The 
notion of a LCFT was developed in the context of quantum and classical field 
theory on curved spacetimes. Thus it is natural to consider Lorentzian manifolds 
as spacetimes. Overwhelming experience from the physics literature suggests 
that it is sensible to restrict our attention to globally hyperbolic0 spacetimes. 
The spacetime morphisms are restricted to ensure compatibility between the 
metric and causal structures on the source and target manifolds. 

^ The relevant notion of spacetime support will be discussed in detail in Sect. 15.3.51 
^ A globally hyperbolic Lorentzian manifold can be identified in several equivalent ways. Per- 
haps the simplest definition is the existence of a Cauchy surface, which is a surface intersected 
exactly once by every incxtcnsiblc timclikc curve I106II49| . 
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Definition 2. Given two oriented and time oriented Lorentzian n-manifolds 
(M,g) and {M',g'), a smooth map x- M ^ M' is called a causal isometric 
embedding if (i) x '^'^ open embedding preserving orientation, (ii) it is an 
isometry, x*9' = 5; preserving time orientation, and (Hi) if any two points 
^tU ^ x(-^) can be joined by a causal curve in M' then they can be joined by a 
causal curve in x(M) (causal compatibility). 

Let &{obS)X)p denote the category of oriented and time oriented globally hyper- 
bolic Lorentzian spacetimes as objects and isometric embeddings as morphisms. 
Let ©[obi3t)p^ denote the subcategory where the morphisms are restricted to 
causal isometric embeddings. 

To express the causality property, it is convenient to introduce a tensor prod- 
uct on Poisson algebras and the notion of independent subsystems. First, though, 
a few words about tensor products in categories. The necessary formal setting is 
that of a symmetric monoidal category [112j . Such a category is equipped with 
a self-bifunctor {A, B): A® B called tensor product (also monoidal structure), 
as well as an identity object. This product is required to satisfy some identities, 
expressed as commutative diagrams, which guarantee that associativity holds, 
that the identity object as is an identity for the product, and that there existence 
of a canonical isomorphism A®B = B®A. Categorical products and coproducts 
(Sect.|B|) are well known examples of tensor products. In fact, we do not give the 
detailed list of the above axioms because the tensor products that will be used in 
this paper are all constructed by putting some extra structure on an underlying 
(co)product in a way that makes the needed axioms manifest. Like the underly- 
ing (co)products, our tensor products will be equipped with canonical inclusion, 
A,B ^ A® B, or projection. A® B — )■ A,B, morphisms. An important dis- 
tinction of a tensor product from a (co)product is the lack of universality. That 
is, the existence of a pair of morphisms A — )■ C and B ^ C (or C —?' A and 
C ^ B) docs not guarantee the existence of a canonical morphism A® B ^ C 
{C A (E) B), unlike for a coproduct (or product). However, since our tensor 
products are based on underlying (co)products, which such a morphism exists, 
it is canonical. A functor between two tensor categories that the tensor product 
structure is called a tensor functor (also a symmetric monoidal functor). 

Definition 3. Let (F, {}f) o,nd (G, {}g) be Poisson algebras, their independent 
subsystems (tensor) product is defined by 

{F,{}f)®{G,{}g) = {F®G,{}), (1) 

where F®G is the tensor product (coproduct) of commutative algebras to which 
the Poisson bracket is extended by the rule {F ® \,\ ® G} = 0. The product 
Poisson algebra is equipped with the canonical inclusion morphisms 

{F, Uf) > [F, Uf) ® (G, Ug) < (G, Ug) , (2) 

which are given respectively by fi i-^ fi ® 1 and f2 ^ ^ ® f2- 

Definition 4. Two subalgebras (^i,{}i) ^ (-F, {}) i — 1,2, of a larger Poisson 
algebra are said to be independent subsystems if the inclusion morphisms fac- 
tor through the independent subsystems according to the following commutative 
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diagram: 

(Fi,{}i) >(Fi,{}i)®(F2,{}2)^ (i^2,{}2) , (3) 




[F,{)) 

where the horizontal morphisms are the canonical inclusions, the diagonal ar- 
rows are the given subalgebra inclusion morphisms, and the vertical dotted line 
morphism is fi® f2 ^ /1/2 Bind is a canonical injective Poisson homomorphism. 

The axioms classical field theory can now be succinctly stated as follows. 

Definition 5 ( [5j ) . A locally covariant classical field theory satisfies the fol- 
lowing axioms: 

Isotony It is a covariant functor, T : ^Xobf^X))? ^ *PoisSj. 
Time Slice The image J-(x) of x - M ^ M' is a Poisson isomorphism whenever x(-^^) 
contains a Cauchy surface of M' . 
Causality The images of morphisms J-ixi) ■ J' (Mi) — > J-{M) are independent subsys- 
tems of T{M), whenever the images of morphisms Xi- Mi — > M, i = 1,2, are 
spacelike separated in M . 

The Causality axiom may also be rephrased in terms of tensor products. We 
can introduce the disjoint union (Mi, (71) U (A/2,52) = {Mi U M2,gi © 32) as 
a tensor product in (&lobS}r)p^. It is based on the disjoint union (coproduct) or 
manifolds, such that disconnected components are considered spacelike separated 
from each other. As for the independent subsystems product on Poisson algebras, 
the canonical dotted line morphism in the following diagram does not always 
exist: 

(Ml , gi ) > (Ml , .91 ) U (M2 , .92 ) < ( A/2 , 32 ) ■ (4) 




{M,g) 

It exists only if the images of Mi are spacelike separated in M. Thus, the Causal- 
ity axiom simply states that the LCFT functor J^: &iobSj[)p^ — > *PoisSj is a tensor 
functor. 

It should also be clear by now that the above formalism explicitly ties the 
causal structure of a field theory to a fixed (non-dynamical) Lorentzian metric. A 
Lorentzian metric provides a notion of a causal relation, which appears crucially 
in each of the named axioms. What about field theories that have a dynamical 
Lorentzian metric or have no Lorentzian metric naturally associated to them at 
all? It is worth taking a step back and examining why a Lorentzian metric is 
important in the first place. 



2.2. Linear and semilinear PDEs. We start by considering examples of field 
theories that can be shown to actually satisfy the LCFT axioms. There are 
two prominent families of examples: free waves field and globally well-posed. 
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interacting wave fields without derivative couplings. Prototypical representatives 
of each family are the free scalar field, obeying the Klein-Gordon equation 

U<j) - m^(j) = 0, (5) 

and the (/)^-interacting scalar field, obeying the semilinear equation 

□0 - TO^0 - A<?!)^ = 0. (6) 

By wave field, we mean a field theory, defined on a Lorentzian manifold, whose 
equations of motion have the corresponding d'Alambertian wave operator □ as 
the principal symbol. The principal symbol is roughly the differential operator 
consisting of the highest derivative order terms of the PDE (cf . Eq. [TT|) . For the 
wave operator □, the principal symbol is essentially g'^'^df^di^, where g'*" is the 
inverse Lorentzian metric, which is contracted with two coordinate derivatives. 
As is usual, the term free means that the corresponding equations of motion 
are linear. The absence of derivative couplings is taken to mean the absence in 
the equations of motion of non-linear terms with derivatives. Such equations are 
termed semilineai^ Finally, well-posedness is taken to mean, in particular, that 
solutions corresponding to arbitrary smooth initial data (obeying appropriate 
boundary conditions, of course) do not form singularities in finite time. Well- 
posedness of linear equations can be proven under quite general assumptions, 
which is not the case for semilinear ones. 

These field theories depend on the background Lorentzian metric through 
the principal symbol of their equations of motion. The exclusion of interactions 
with derivative couplings, ensures that this principal symbol is always the same 
as that of the wave operator □. Since the metric is Lorentzian, these partial 
differential equations are special cases of hyperbolic PDE systems. Roughly, hy- 
perbolicity is an algebraic and geometric property of the principal symbol, which 
is exploited by standard PDE theory to prove analytical results like (a) local-in- 
time well-posedness of the Cauchy initial value problem (this entails the existence 
of solutions for arbitrary smooth initial data, uniqueness of solutions, continu- 
ous dependence of the solution on the initial data) , (b) finiteness of the speed of 
propagation of disturbances. In particular, as a generic result, the maximal speed 
of propagation in any direction happens to be given by the null directions of the 
background Lorentzian metric. Global well-posedness can also be proven under 
generic circumstances for linear PDEs and under some additional hypotheses for 
semilinear ones. 

These results are important because, roughly, global well-posedness and fi- 
nite propagation speed jointly imply the validity of the Isotony, Time Slice and 
Causality axioms for the corresponding classical field theory. The Causality ax- 
iom, and hence the notion of causal structure, is particularly closely related to 
finite propagation speed. We shall see in Sect. l7.4l in more detail how these results 
from PDE theory can be used to construct LCFTs and verify these axioms. 

Following the logic of the last paragraph, it makes sense to look for a gen- 
eralization of the notion of causal structure, and corresponding generalizations 
of LCFT axioms, by studying the way PDE theory establishes the properties of 
well-posedness and finite propagation speed in more general hyperbolic systems. 



^ Actually, semilinear equations allow non-linear terms with derivatives, as long as they are 
not of the highest present order. 
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Of particular interest are PDE systems where these properties are not directly 
tied to a background, non-dynamical Lorentzian metric. As will be shown in 
the following sections, such generalizations can be found in the PDE theory of 
quasilinear hyperbolic systems, of which GR is a special case. 



3. Quasilinear Hyperbolic Systems 

In this section we formulate several aspects the theory of systems of partial 
differential equations (PDEs) in a geometrical way, namely in terms of jets. This 
approach is not completely standard in the analysis or mathematical physics 
literatures, but it does have some advantages. 

Eventually, we would like to make contact with locally covariant field theory, 
which assigns algebras to spacetime regions in a functorial way (a diffeomorph- 
ism-covariant way). We would like to construct these algebras as algebras of 
functions on spaces of solutions of some PDEs. Thus we would like to assign 
PDE systems to manifolds in a functorial way as well. It turns out that this 
is more conveniently done using natural bundles (to be introduced in Sect. [5]) 
and jet bundles. Moreover, it is also more convenient to discuss integrability 
conditions and equivalence of PDE systems in the language of jets. Finally, the 
variational bicomplex on jet space is the natural setting for conservation laws 
and the covariant symplectic structure (to be discussed in Sect. [5]). 

Within this framework, we define PDE systems of quasilinear hyperbolic type. 
We rely mostly on the notion of symmetric hyperbolicity of first order systems 
(following |50p. but also briefly comment on the related notion of regular hyper- 
bolicity. 

From this point, we will be discussion finite dimensional spacetime manifolds, 
bundle manifolds and infinite dimensional manifolds of sections and solutions, as 
well as morphisms between them. For later convenience and to fix some notations 
it is helpful to define the following categories. 

Definition 6. Let 9Jlan denote the category of smooth manifolds (of possibly in- 
finite dimension) and COtarie the subcategory where morphisms are restricted to 
open embeddings. The morphism are always C°° maps. Denote also by 9Jtan" 
and 9Jtan" the respective subcategories of manifolds of fixed dimension n — 
0, 1, 2, . . . , oo. The finite dimensional manifolds are taken to be oriented. 

Let *8nc)[ denote the category of smooth, finite dimensional bundles over mani- 
folds and ®nc)[e be the subcategory where morphisms are restricted to open embed- 
dings of the total space. These categories are fihered over manifolds, 11: *Bnc)[ — > 
aton, with H: [E ^ M) ^ M. Let <Bn£)l(A/) denote the subcategory o/ <Bnt)[ 
of bundles over the base manifold M with base fixing morphisms (projection to 
the base is always id : M M). Thus, in an obvious way, we have defined a 
functor from manifolds to categories, *Bn5[: 9Jlan — > Cat. 

Let 23Snc)I denote the category of vector bundles, a subcategory o/*Bni)[, and 
QJ^nOlg — 2J*BnJ)(n*Bn5[e. We similarly introduce the notations 2J*Bn5[(M) and 
^^nUeiM). 

We shall adopt 11 as the generic notation for forgetful functors. More gener- 
ally, from time to time, we shall reuse functor names as long as they can be 
distinguished by the domain category and are similar in purpose. 
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We shall not be specific about the precise definition of the class of infinite 
dimensional manifolds in question. We presume that such details can be added 
when applying the current formalism to specific situations. Such functional an- 
alytical questions have been recently tackled in earnest in (cf. also [19]). 

A prominent example of an infinite dimensional space that we will treat as a 
manifold is the space r{E) of all smooth sections of the vector bundle E — M, as 
well as some infinite dimensional subspaces thereof, like the space of all solutions 
of a PDE system. 

3.1. Jet bundles and systems of PDEs. This section outlines the description of 
PDEs as submanifolds of the jet bundle. Jet bundles are briefly introduced in 
Sect. El where also notation is fixed (not all of it being completely standard) and 
standard literature references are given. Such a description of PDEs is more in- 
trinsic than than the usual one in terms of equations, but is essentially equivalent. 
This approach is well known in the geometric and formal theory of differential 
systems |1021l?7] . It will be later combined with the notion of natural bundles to 
define locally covariant field theories. The discussion of natural PDE systems in 
connection with natural bundles will be postponed to Sect. [6] Below wc consider 
only fixed base manifolds. 

From now on, fix M to be finite dimensional manifold and let n = dimM. 
Also fix a vector bundle F -> M. We refer to M as the spacetime manifold and 
to F as the field bundle. 

We restrict our attention to regular PDEs in the following sense. 

Definition 7. A PDE system £ of order k is a smooth, closed sub-bundle (in 
the «8nO[ sense) of J^F M, £ C J^F. 

Note that £ need not be a vector sub-bundle of J^F. The above definition may 
seem unfamiliar to some, but can be cast in more recognizable form using the 
following 

Proposition 1. Given a PDE system £ of order k, there exists a vector bundle 
E — > M , a smooth sub-bundle E' C_ E containing the zero section of E ^ M , 
and a smooth base fixing bundle morphism f: J^F — >■ E (in the *BnI)[ sense) 
such that the image of f is contained in E' , the image of f is transverse in E' 
to the zero section of E and £ is precisely the preimage of the zero section, that 
is, £ satisfies f — 0. 

The proof follows from basic differential topology, up to a global topological ob- 
struction ^SBl §7]. Clearly, the equation form is not unique. For instance, applying 
any invertible transformation to the equations / = gives another equation form 
/' — 0, which describes exactly the same PDE system. 

We refer to — >■ M as the equation bundle and to / or the pair (/, E) as the 
equation form of the PDE system £. A section 0: M ^ F, also referred to as a 
field configuration, is said to satisfy the PDE system £ if the fc-jet prolongation 
of (f) is contained in £, j''(j){x) & £x C J^{F,M). Equivalently, is a section 
oi £ ^ M. We denote the space of all solution sections by S{F) C T{F) or 
Ss (F) when the PDE system needs to be mentioned explicitly. Using the above 
proposition, we can equivalently say that is a solution of the PDE system £ if 



(7) 
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Expressing the fc-jct in local coordinates, j''(j){x) = {x, (jf" {x) , di(t)°' {x) , . . .), it is 
clear that f{x, 4>°'{x), di<j)°'{x), . . .) = is a system of partial differential equations 
in the usual sense of the term. Starting with a PDE system in the usual sense, its 
geometric form as a sub-bundle of the jet bundle can be obtained by a converse 
of the above lemma. At this point, the regularity assumptions on both £ and 
/ become important. Namely, the transversality properties of / ensure that the 
zero set of / = is a submanifold of j'^F and vice versa. 

The linear and affine structures on J^'F give us the possibility of defining the 
notion of linear and quasilinear PDE systems. 

Definition 8. A PDE system £ C J''F is called linear if £ ^ M is a vector 
sub-bundle of the vector bundle J'^F — ^ M. The PDE system is called quasilinear 
if £ ^ J''~^F is an affine sub-bundle of the affine bundle J^F J''~^F. 

The connection to the usual meanings of these terms can be seen through 
adapted equation forms. 

Lemma 1. The PDE system £ C J^F is linear iff it has an equation form 
{f,E), where f: J^F E is a morphism of vector bundles over M. 

The PDE system £ C J''F is quasilinear iff it has an equation form {f,E), 
where f : J'^F ^ E is a morphism of affine bundles, which fits into the commu- 
tative diagram 

jkp — [ ^p 

(8) 

jk-ip ^jyj- ^ 

where the vertical maps define the affine bundles, with the vector bundle E ^ M 
naturally considered an affine one. 

The proof is immediate. Alternatively, the quasilinear case can be cast into the 
form of a base fixing aSine bundle morphism /: J''F — ^ {E)'^~^, where both 
bundles are over J'^~^F. Such equation forms are called adapted. 

In the more common language of adapted local coordinates, the conditions 
of linearity and quasilinearity are expressed as follows. Consider adapted local 
coordinates {x^,va) on the equation bundle E, {x'^, u°') on the field bundle F, 
and the corresponding {x'^,uf) on the fc-jct bundle J^F. Let (j): M F he a 
field configuration, then its fc-jet in local coordinates is j'^(j){x) = {x^ ,di(j)°'{x)). 
The above lemma asserts the existence of an equation form that looks like 

fL{^)dir{x) = o. (9) 

Note that this equation is linear in ^(a;) and its derivatives and that the coeffi- 
cients /^a(a;), with multi- indices /, depend only on the base space coordinates 
X. On the other hand, for a quasilinear equation, the lemma asserts the existence 
of an equation form that looks like (with |7| = k) 

fL{x,3^-^<t>{x))dir{x) + = 0. (10) 

Note that, in the linear case, the fact that the coefficients of /: J^{F, M) E 
only depend on the base space coordinates x is captured by the requirement 
that it is a morphism of vector bundles over M . In the quasilinear case, the 
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coefficients of / can obviously depend on both x, 4>{x) as well as all derivatives 
di(j){x) up to order \I\ — k — 1, which is captured by allowing /: J'^F — > (E)''^^ 
to be a (base fixing) bundle morphism over J^~^F. It is worth remarking that 
any linear PDE system is also naturally quasilinear. 

Recall that the affine bundle J''F J^~^F is modeled on the vector bundle 
{S^T*M®M Ff-'^ J^-^F. Therefore, an adapted equation form (/, E) of a 
quasilinear PDE system £ C J^F naturally singles out a section 

/: J^-^F ^ {E®M F* ®M S^TMf-^. (11) 

In local coordinates, / corresponds to the coefficient /^^ of the highest derivative 
term di<j}{x) with |/| = fc in Eq. (llOp . This section / is called the principal symbol 
of the given equation form of £ . 

Below, we will be mostly concerned with first or second order PDE systems. 



3.2. Symmetric and regular hyperbolicity. The notion of hyperbolicity for a PDE 
system is strongly linked to the ability to formulate it as an initial value problem. 
Locally such a formulation is constrained by the existence of so called charac- 
teristic surfaces or characteristic covectors, which are defined in detail below. 
Various notions of hyperbolicity are then stated either in terms of geometric con- 
ditions on the locus of characteristic covectors or in terms of equivalent algebraic 
conditions on the principal symbol of a special equation form of the system. 

An initial value formulation consists of converting the PDE system into an 
infinite dimensional ODE as follows. Suppose we are given an n-dimensional 
manifold M that can be smoothly factored as M = K x S', where 5* is a manifold 
of dimension n— 1, and the projection on the first factor is denoted hyt: M — ^ R, 
the time function. Then the space of restrictions of field configurations to a level 
set of t (always a codimension-1 surface diffeomorphic to S) forms an infinite 
dimensional space, which can be considered as a fiber of a smooth, infinite di- 
mensional bundle over R. If the PDE system can be turned into an ODE system 
on this bundle, one can use Picard iteration and a Gronwall-type lemma to prove 
local-in-i existence and uniqueness of solutions. Various notions of hyperbolicity 
essentially correspond to sufficient conditions under which the above construc- 
tion can be carried out. Some of these conditions are local (referred to simply 
as hyperbolicity) and some global (referred to as global hyperbolicity). Local con- 
ditions put restrictions on the principal symbol of a PDE system and on the 
germs of the leaves of the 5'-foliations of M to which the above construction 
would be applicable. The global conditions require M to factor in a way similar 
to above, with the ^-leaves satisfying the corresponding local conditions. Strictly 
speaking, when the 5'-leaves are non-compact or have a boundary, boundary con- 
ditions may have to be supplied, on top of those already mentioned, to ensure a 
well-posed initial value problem. Such boundary conditions will not be discussed 
in this work. 

To obtain some necessary conditions on the geometry of initial value surfaces, 
we need to first look at the linear and affine geometry of fc-jet space in the 
presence of a preferred codimension-1 subspace C T^M at a fixed x € M . 
This subspace can be thought of as tangent to a putative initial value surface S 
passing through x. For convenience of notation, we promote Tx to be the fiber of 
a codim-1 vector sub-bundle r C TM (also called a tangent plane distribution). 
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However, all subsequent calculations will be purely local and go through equally 
well if T is only defined in a neighborhood of x, or just at x itself. In particular, 
the distribution t need not be integrable. 

Let V <Z T* M he the 1-dimensional vector sub-bundle of covectors annihilat- 
ing r (the conormals to the putative surface S). In turn, the l-dimensional vector 
bundle ly singles out the 1-dimensional vector sub-bundle v^'' = v®mv®m • • • C 
S''T*M, and finally the vector sub-bundle N ^ i^'^'' (8}m F C S^T*M ®m F. 
Recall that the affine bundle J^F J^~^F is modeled on the vector bundle 
{S^T*M®M Ff-^. Note that the fibers of J^F J'^-^F are foliated by affine 
planes parallel to the fibers of the vector bundle {N)^~^ . We define J^~^'^ to 
be the leaf space of this foliation. In other words, the short exact sequence of 
vector bundles 

^ ^ S''T*M ®M F {S^T*M ®m F)/N ^0 (12) 
induces the following affine bundle projections 

jkp ^ jk-i,±p ^ jk-ip^ (13) 

where the second projection defines an affine bundle with fibers modeled on the 
vector bundle {S^T*M (S)m F/N)^~^ and the first projection defines an affine 
bundle with fiber modeled on the vector bundle [N)^'"^'^ ^ with the latter nota- 
tion meaning the pullback of the bundle N — > M along the canonical projection 
jk-i,± _^ j^j^ ^ more in depth discussion of the bundle J^'^F can be found 
in [ina Ch.6]. 

The above constructions are easily illustrated in local coordinates. Consider 
local coordinates x^ = {t, s*) on a neighborhood oi x G M and local adapted 
coordinates {t,s'',u°') on the corresponding neighborhood of F^ C F. Suppose 
that the t coordinate is chosen such that the level set t = t{x) is tangent to the 
plane C T^M, in other words dt{x) & C T*M. The k- and (fc - 1, -L)-jets 
of a section (p: M F can be represented by 

j''<Pix) = {x\djr{x)), \J\<k, (14) 

= {t,s\didir{x)), \i\ + i<k, (15) 
= {x\dj,r{x),dtdj,r{^),---,dt'di,_,ri^),d'^r{^)), m 

\Ii\ <k~l, 

where the multi-indices in ,/ range over all indices of x^ , while the multi-indices 
of / and // range only over the indices of s* and not t, and by 

/-i>^^«(x) = {t,s\du,r{x),dtdi,r{x),...,dt'di,_,r{x)), (i?) 

<k-i. 

In other words, the affine bmidle projection j^4){x) i-> i^~'^^^(j){x) simply dis- 
cards the last component, which is the highest order derivative d^4f'{t^s) along 
the i-direction, which is transverse to the level set of t passing through x. 

In order to convert the PDE system into an ODE with respect to the t 
coordinate, still working locally at a; G Af, S*, given a section (f): M ^ F that 
satisfies the PDE system, i^(f){x) G £x, we need to uniquely determine the highest 
order t-derivative of the dynamical part of (j) at function of all other 
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derivatives of equal or lower ordefl So, geometrically, each fiber of the affine 
bundle {J^F)''~^ ^jk-i,±p^jk-i ^gg^jg intersect £ exactly once. On the 
other hand, to express this condition algebraically, we need to pick an adapted 
equation form (/, E) for this quasilinear PDE system, with principal symbol 
/. Pick sections u: M J''F, : M ^ J'^'^'-^F and u: M J'^^^F such 
that M and u i-^ u under the affine bundle maps J'^'F — > J^^^'^F and 

J^F — >■ J^~^F respectively. Then by the affine structure of the /c-jet space, any 
other fc-jet section v that projects to u-^ is given hy v — u + p'^'' €5 ip, with 
sections ■0: M — >■ _D and p: M v. li we demand that v is actually a section of 
£ C J^F, then a simple calculation shows that 

f{v) = Mp'"' ® ^) + /(^) = (18) 
=^ {h-p^'')^ = -f{u). (19) 

The last equality is a linear equation for -0. For it to have a unique solution, 
as we demanded above, the linear map /a • p'^^ : F E needs to be invertible. 
This discussion motivates the following definition. 

Definition 9. Consider a vector bundle F M and a quasilinear PDE sys- 
tem, £ C J''F, with adapted equation form {f,E) and principal symbol f. Given 
u e j'^^^F, a covector p £ (T*M)g^^ is called u-non-characteristic (or just 
non-characteristicj if the contraction of the principal symbol with p®*^, (/^ • 
p^ky^ (i^)|~^ — >• (_E)g~^, is invertible. Otherwise, p is said to be u-characteristic 
( or just characteristic J. 

Recall that in local coordinates on F, a (fc - l)-jet u e J''-^F IS rep- 

resented as M = (x'^,Uj), \I\ < k — 1. For linear and semilinear systems, the 
principal symbol depends only on the base space coordinates x\ This implies 
that for such systems a covector p <£ T*M can be decided to be characteristic 
without looking at the field value or higher jet components w". This is the 
usual situation in relativistic field theory with a fixed background metric, where 
characteristic covectors coincide with null covectors of the background metric. 
On the other hand, for quasilinear PDE systems, the principal symbol fu may 
depend on as well as higher jet components u'j up to order fc — 1. In this 
sense, the notion of a characteristic covector becomes field dependent. 

As we shall see later on, the geometry of the locus of characteristic covectors 
of a PDE system is closely related to the causal structure of the corresponding 
classical field theory, in particular to the domain of dependence and finite prop- 
agation speed results. Compare now the equations of motion of relativistic field 
theory (say the Standard Model) with a fixed background metric and coupled 
with GR, which provides a dynamical metric. As remarked above, characteristic 
vectors and hence the causal structure of the theory on a fixed background is 
field independent, since the equations of motion constitute at most a semilin- 
ear PDE system. On the other hand, coupled to GR, relativistic field theory 
becomes quasilinear, since all principal symbols depend on the metric, which 
is dynamical. Thus, GR with any of the Standard Model matter coupled to it 
constitutes a system with field dependent causal structure. This statement will 
be made more precise in the next section. At the very least, we expect the causal 



* In the classical PDE literature this is known as putting the equation in Cauchy- 
Kovalevskaya form |42II34| . 
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structure of GR to be significantly different (and more complicated) than that 
of other relativistic field theories. 

Characteristic covectors are obstructions to converting a PDE system into 
ODE or Cauchy-Kovalevskaya form. An initial value formulation with data on a 
codim-1 surface t: S C M can be achieved if the highest S'-transverse derivatives 
of the unknown section could be solved for in terms of data ip: S ^ l*J'^~^F. 
For this it is necessary that a no n- vanishing conormal section p: 5 — > l*T*M is 
everywhere non-characteristic with respect to ip. Such a pair [S. 93) is referred to 
as non- characteristic initial data. 

However, being non-characteristic is not a sufficient property to set up a well- 
posed initial value problem for a given initial data set (5*, (p) . This is where the 
notion of hyperbolicity comes in. There are several different notions of hypcr- 
bolicity. We will only discuss two, which are sufficient for our purposes. First, 
we define symmetric hyperbolicity and then make some comments about regular 
hyperbolicity. 

For a first order quasilinear system, the geometry of the principal symbol 
simplifies. It can be expressed as a morphism 

f:{F®MT*Mf^{Ef (20) 

of vector bundles over J^F = F. Given a section p: F ^ {T*M)° and ip: F ^ 
(F)", we get a section ^ = f{p®i)) = {f ■p)tp of {E)° — )• F. In local coordinates 

{x^,u°') on F and (x*,w^) on E, we have 

U{x,u) = fX^{x,u)p'{x,u)^''{x,u). (21) 

For the next definition, we introduce some new notation. Recall that yl"M — )• 
M is the bundle of densities on M. Given any vector bundle F — >■ M, we denote 
the bundle of V -valued densities hyV = A^M^mV . The bundle of dual densities 
is denoted by V* = A^M ®m V* ■ Densitized symmetric powers will be denoted 
by S^V — A"-M (E)m S''V. The defining property of an orientable manifold is 
the existence of nowhere vanishing sections oi A"M M. An oriented manifold 
chooses a privileged class of positive densities, which are nowhere vanishing and 
such that any two positive densities are related through multiplication by an 
everywhere positive scalar function. 

Definition 10. Consider a first order, quasilinear PDE system, on F ^ M , with 
M oriented and with adapted equation form (/, F*) and principal symbol f . This 
PDE system is said to be symmetric hyperbolic if for each u G F there exists a 
covector p e (T*M)? such that 

fu-pe {S^F*)l c {A^M ®M F* 0M F*)l, (22) 

that is, fu-p is a symmetric bilinear form on the fiber (F)^, and moreover that 

fu -p is positive definite with respect to the orientation on M. Such a covector p 
called ■u-spacelike (or just spacelike j and future oriented, while —p is ■u-spacelike 
but past oriented. 

Decoding the above definition in local coordinates (x^u") on F, we have that 
Piflij{x, u)xl>°-'il:'' > (as an element of (/1"M)°^ with respect to the orientation 

on M) with p e (T*M)0^_^) and for all V S {F%,uy 
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Note that positive definiteness is stronger than invertibihty. So every spacehke 
covector is also non-characteristic, but not every non-characteristic covector is 
necessarily spacelike. An initial data set {S, tp), with l: S G M, is called spacelike 
if a non- vanishing conormal section p: S l*T*M is everywhere </3-spacelike. 
It is for spacelike initial data that the theory of symmetric hyperbolic PDEs 
establishes local well-posedness. 

Remark 1. Our definition of a symmetric hyperbolic system is slightly different 
from the standard one, where the principal symbol is valued in S'^F* rather than 
the densitized S^F* . For a manifold with a fixed orientation, as we are consid- 
ering, the difference consists of tensoring with a positive density. Densitizing the 
equation form turns out to be more convenient in the current setting, in particu- 
lar in Sect.[SJ However, densitization is also natural in the original context where 
symmetric hyperbolicity is used, the construction of an energy norm by integrat- 
ing over a future oriented, spacelike surface t: E C M, = i* tr /(-!/', V"); 
where the trace naturally converts the principal symbol into a current density 
(the single contravariant index of / is contracted with one of its n antisymmet- 
ric covariant ones), which can be naturally integrated when pulled back to the 
codim-1 surface S. 

Remark 2. Note also that our usage of the word spacelike differs from the stan- 
dard one one in Lorentzian geometry, where our spacelike covectors would be 
referred to as "timelike". However, this terminology is consistent with the lit- 
erature on hyperbolic PDEs. Moreover, outside pseudo-Riemannian geometry, 
there is no natural identification between vectors and covectors. On the other 
hand, covectors are still naturally identified with codimension-1 tangent planes. 
As expected a spacelike plane consists of spacelike vectors we merely extend this 
terminology to a corresponding covector. The terminology introduced for vectors 
below is standard in either literature. 

It is convenient to introduce the following geometric notions as well. 

Definition 11. Consider u £ F. A tangent vector v e (TA/)? is called causal 
and future directed if p ■ v > for every spacelike and future oriented p € 
(T*M)?. // the inequality is reversed, p ■ v < 0, then v is called causal and past 
directed. 

As noted above, the notion of symmetric hyperbolicity applies only to first 
order, quasilinear systems. This is not too strong of a restriction in practice, 
as often higher order PDE systems can be reduced to first order ones by in- 
troducing extra fields (increasing the dimension of the field bundle fibers). In 
particular, this exercise was carried out for GR, all Standard Model fields, as 
well as relativistic hydrodynamics in |50] . On the other hand, many equations 
obtained as the Euler-Lagrange equations of an action are not directly in sym- 
metric hyperbolic form. Such equations can actually be conveniently treated 
directly, without reduction to first order form. The relevant notion which comes 
with a well-posedness theory is regular hyperbolicity , which was developed much 
more recently by Christodoulou in f3T| . It would be ideal if the class of regularly 
hyperbolic systems, upon reduction to first order form, were included in the class 
of symmetric hyperbolic ones. Unfortunately, this does not appear to be the case 
due to subtle differences between the two definitions [TS]. On the other hand, the 
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similarities between symmetric and regular hyperbolicity include similar defini- 
tions of characteristic and spacelike covectors and similar conal properties of the 
geometry of their loci, which are explored in the next section. These properties 
are all we need for the purposes of this paper. Therefore, in the sequel, we refer 
only to results for symmetric hyperbolic systems, even though similar results 
can be obtained for regular hyperbolic ones. 

At this point it is worth making a few comments on the so-called energy 
methods used to establish well-posedness for symmetric and regular hyperbolic 
systems of order k = 1, 2. A central role is played by a family of local, horizontally 
conserved, positive definite, coercive energy current densities s'^ , parametrized 
by r > 0. Locality means that it is a section of the horizontal (n — l)-form 
bundle over J'^F, : J^^^F ^ {A'''-^Mf-^. Equivalently, e ^2"^i^0(i^) 
and is projectable to J^~^F . Horizontally conserved means that d(j'^0)*e'^ = 
(j'^(/))*dh£^ = 0, when (/>: M — >■ f is a solution of the given PDE system. This 
conservation criterion can actually be relaxed up to lower order terms, that is, we 
allow (j*'0)*dhe'^ = l{j'^~^4>) +Tl^{j^4') to be a function of (j) and its derivatives. 
Here, both terms on the right hand side are of lower order: 7 because it depends 
only derivatives up to order k—1 and /? because it is proportional to r, which can 
be made arbitrarily small. Positive definite means that the puUback L*[{j'^(j))*e'^] 
onto a future oriented, spacelike codim-1 surface l: U C M gives a positive 
density, so that the energy integral E'^[(j)\ = f^i'*[{j''(f>)*£^] is also positive. 
Coercivity means that the algebraic structure of e"^ is such that it allows E"^ [ip] 
to dominate || some i^-norm on the restrictions of (fc — l)-jets to U, 

in order to prove an inequality of Gronwall type, 

llj'-Vlli;. <Ce«||j^-Vko, C>0, (23) 

where the surfaces St are spacelike surfaces with common boundary, whose 
interiors foliate a domain in M. Such domains are called lens-shaped and, as 
we shall see later on, they essentially determine domains of dependence. This 
inequality is crucial in the proofs of well-posedness. 

The main difference between the symmetric and regular hyperbolic equa- 
tions is in how the energy current density e'^ is obtained. In the symmetric 
hyperbolic case, e"^ is constructed directly from the principal symbol, namely 
e'^[(j)\ = e^*/'^tr/0(^, 0), where i is a time function whose spacelike level sets 
foliate a lens-shaped domain. In the regular hyperbolic case, Christodoulou ex- 
tends Nother's theorem to convert a timelike vector field into a corresponding 
conserved (up to lower order terms) "symmetry current density" e, which satis- 
fies all the desired properties of an energy current density. See Sect. 5.0 of |31j . 
as well as [K], for a more detailed comparison of the two cases. 

From the above discussion, it follows that it is not so much the structure 
of a special equation form that defines the PDE system that matters directly. 
Rather what is important to establish well-posedness using energy methods is 
the presence of (almost) conserved, local energy current densities. The algebraic 
structure of the PDE system is then important in so far as it gives rise to such 
energy current densities. In the geometric theory of PDE systems (or equiv- 
alently of so-called differential systems) the study of such (almost) conserved 
currents (or conservation laws) has been named characteristic cohomology |28l 
iniin] (this meaning of the overloaded term characteristic is distinct from its use in 
characteristic covector). For some PDE systems, these conservation laws can be 
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classified exhaustively, in an intrinsic manner (that is, independent of the equa- 
tion form used to define the PDE system) |20|. This observation leads one to 
speculate that a more general notion of hyperbolicity can be defined in terms of 
intrinsic invariants of a PDE system given by its characteristic cohomology, such 
that both symmetric and regular hyperbolicity become special cases thereof. 

3.3. Prolongation, integrability, hyperbolization. One reason to discuss PDE sys- 
tems as submanifolds of a jet bundle is independence of a particular equation 
form. Any two equation forms are equivalent if they define the same PDE system 
manifold. We should specify our notion of equivalence. 

Definition 12. Consider two field bundles Fi — > M, i = 1,2, and two PDE 
systems £i C j'^'Fi. Denote the corresponding spaces of smooth solution sections 
by Si{Fi). The PDE systems £i and £2 are said to be equivalent if there exist 
bundle morphisms eij : ,P^Fi Fj , i ^ j , such that 

(f>i S Si{Fi) and (pj = eij o implies £ Sj{Fj), (24) 

as well as that ei2 o j'l and 621 o j''^ are mutual inverses when restricted to the 
solution spaces Si{Fi) and S2{F2). 

We can easily extend the notion of equivalence to equation forms of PDE systems. 
In that case two different equations forms that define the same PDE system 
manifold are trivially equivalent. Note that neither the field bundles nor the 
orders of the PDE systems need to be same for equivalence to hold. 

Let us restrict to the case that will be of importance in a later section, namely 
oi Fi = F2 ^ F and ei2 and 621 respectively equal to the canonical projections 
jhp ^ F and J'^F F, which are in a sense trivial. In this case, it is 
certainly sufficient that £1 — £2 for equivalence to hold, but it is not necessary. 
In fact £1 and £2 could be of different orders. To obtain necessary conditions for 
equivalence, we need to consider prolongation of PDE systems and the possible 
resulting integrability conditions. 

A discussion of these notions in the setting of the jet bundle description of 
PDE systems can be rather technical. On the other hand, the theory of equiva- 
lence of PDE systems formulated in these terms has become quite mature and 
has yielded some important results. The technical details of this theory can be 
found elsewhere P71ll02j . Below we give a brief non-technical introduction to 
this theory and state some simplified results relevant for hyperbolic systems. 

The step by step derivation and inclusion of integrability conditions into a 
PDE is called prolongation. It is easiest to define prolongation in equation form 
and in local coordinates. Consider an equation form (/, E) of a PDE system 
£ C J^F, as well as local coordinates [x^ ,u°') on F and (a;*,^;^) on E. If the 
section (p: M ^ F satisfies the PDE system, we have the following system of 
equations holding in local coordinates 

fA{x\djr) = 0- (25) 

These equations hold for each point a; G M, therefore when both sides are 
differentiated with respect to the coordinates on M, the resulting equations are 
still satisfied, 

djA{x',djr) = idjA){x',d.jr) = 0, (26) 
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where difA are functions on J'^^^F obtained by pulling back the functions /a 
from J'^F to J^~^^F and applying the horizontal vector field di. These new func- 
tions fiA = difA, together with the old /a ones, constitute the local coordinate 
expression for the equation form {p^ f, J^E), where is the 1-prolongation de- 
fined in Scct.[^ We call the corresponding PDE system £^ — f C J'^^^F the 
first prolongation of £ or also its prolongation to order k + 1. Prolongations to 
any higher order, (p'/, ,PE) and C J^^^F, are defined itcratively. 

Let Pi : J^^^F — ^ J'^F be the canonical jet projection, which restricts to 
pi: £'' — >■ £. Notice that we necessarily have pi{£^) C f, since the prolonged 
system contains the original one as a subsystem. We have just shown that 
sections satisfying £ automatically satisfy £\ and vice versa. In other words, 
SsiF) — '^g! (F) and the two PDE systems are equivalent. However, the inclusion 
Pi£^ C £ may be strict, which would mean that there exist non-trivial integra- 
bility conditions. There exists an equation form {f(Bg,E(BG) for pi{£^) C J^F, 
where {g, G) is an equation form for the integrability conditions. These obser- 
vations provide another sufficient condition for the equivalence of two PDE sys- 
tems, namely that there exists an order I > fci,fc2 such that £[~''^ ~ ^2 '^^ as 
subsets of ,PF . 

Prolongation can be iterated indefinitely. Taking this process to its limit, 
we obtain the infinite order prolongation C J°°F from the equation form 
{p°° f, J°°E), which takes all possible integrability conditions into account. One 
can then show that the equality f f° = £^ , as subsets of J°°F, is both a necessary 
and a sufficient condition for the equivalence of two PDE systems. It is a deep 
theorem of the geometric theory of PDE systems [102,27,53^ that, for any given 
PDE system, there exists a finite order I such that prolongations above that order 
introduce no new integrability conditions. Therefore, this restricted version of 
the equivalence problem can be decided in finitely many steps. Finally, it can 
also be shown that any PDE system is equivalent to one of first order, though 
usually defined on a different field bundle. 

It is a well known fact that many PDE systems of mathematical physics are 
not given directly in symmetric hyperbolic equation form, though for somewhat 
different reasons. The Klein-Gordon equation, though regularly hyperbolic, in 
its variational form is not first order. The Dirac equation, though first order 
in its variational form, does not have a symmetric principal symbol. Maxwell 
equations, in terms of the vector potential, cannot be hyperbolic because of 
local gauge invariance (which spoils uniqueness in the Cauchy problem). The 
Proca equation, again in its variational form. Maxwell equations, in terms of the 
field strength or in terms of a gauge fixed vector potential, are only equivalent 
to a hyperbolic system with additional constraints. 

Consider a quasilinear PDE system £ C J^F with an equation form (/ © 
c,F* © E), where (/, F*) is an adapted, first order, quasilinear, symmetric hy- 
perbolic equation form. We refer to {f,F*) as the hyperbolic subsystem and to 
(c, E) as the constraints subsystem. The constraints are said to be (symmetric ) 
hyperbolically integrable if there exists a first order, linear PDE system £' C J^E 
with adapted (symmetric) hyperbolic equation form {h, E*), the consistency sub- 
system, such that the following identity is satisfied: h o c — q o h, with some 
differential operator q: J^F* — > E* satisfying <?(0) — 0. In other words, for any 
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section : M ^ F we have 

h[cm^q[mi (27) 

and /i[c[(/)]] = if /[</)] = 0. The consistency subsystem {h,E*) is linear in the 
sense that /i[c[(^]] depends linearly on c[(/)], but it may depend non-linearly on 
(j), such that the compound system with equation form (/ © h, F* © E*) is first 
order, quasilinear, symmetric hyperbolic. A PDE system with such an adapted 
equation form (/ © c, F* (BE) is called (symmetric) hyperbolic with constraints. 
When the constraints are hyperbolically integrable, by studying the properties 
of the compound hyperbolic PDE system on F (Bm E, we shall see in Sect. 14.31 
that all the well-posedness results for symmetric hyperbolic systems also apply 
in the presence of constraints. 

We call a PDE system (symmetric) hyperbolizable if it is equivalent to one 
that is (symmetric) hyperbolic with constraints. We call the equivalence map 
that brings a PDE system into such a form a (symmetric) hyperbolization. 



3.4- Examples. Many examples of reductions of relativistic field theories are 
given in Appendix A of f50'. Another source of examples of field theories de- 
scribed by hyperbolic systems with constraints is [58j . In the latter reference a 
slightly different notion of hyperbolicity is used, but the given examples still fit 
into our framework provide they are first reduced to symmetric hyperbolic form. 



4. Chracteristic Geometry, Causality, Domain of Dependence 

The term Lorentzian geometry refers to the study of structures induced on space- 
time manifolds by the presence of a Lorentzian metric. One example of this kind 
of structure are the cones of null vectors. In particular, it is these cones that 
determine causal relationships between points in Lorentzian spacetimes. Below 
we will be similarly interested in the geometry of characteristic covectors of a 
hyperbolic PDE system, which also form cones and also induce a causal order 
on the points of a manifold. We refer to the investigation of the geometry of 
cones of characteristic covectors and related structures as characteristic geome- 
try. Even more generally, we are interested in conal geometry, which is concerned 
with the differential topology of conal manifolds |78ll87j . Each point of a conal 
manifold is smoothly assigned an open cone of tangent or cotangent vectors. 
The study of conal manifolds, referred to below as cone bundles, though still 
an immature field, has the potential to capitalize on and then subsume much 
of the earlier work on causal order on spacetime manifolds pnH^[5nil48j . The 
generalization from Lorentzian cones to more general ones, for the purposes of 
describing causality in quantum field theory has been considered before 
193] . but not in a concrete way. 

Fix a vector bundle F M and a first order, quasilinear, symmetric hy- 
perbolic PDE system £ C J^F on it with adapted equation form (/, F*) and a 
symmetric adapted principal symbol /. We mostly follow 14 with the regards 
characteristic geometry terminology. 
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4-.1. Geometry of cone bundles. 

Definition 13. For u e F, denote by Ff C (T*M)? the set of u-spacelike, 
future oriented covectors. Similarly, denote by C (T*M)? the set of u- 
characteristic covectors. Let a covector be called cocausal, future oriented if it 
belongs to F^ = F® , where the bar denotes closure. Finally, let a covector be 
called inner u-characteristic and future oriented if it belongs to = C,* n F*^ . 
Also, denote 

r® = y r®, r* - y F*^, C* = [j Q, and c* ^ [j CI (28) 

UeF ueF u£F u&F 

It is easy to show that F® C (T*M)° is an open subset and hence a submanifold, 
as well as that dF* = dF® = C* . Moreover, each F® is a convex cone (invariant 
under multiplication by positive scalars). Note that, since the interior of an open 
convex cone is diffeomorphic to a ball, F® — ?> M is actually a smooth bundle, 
in the sense of SnOt. On the other hand, F® has much more structure than a 
generic smooth bundle. To take this structure into account we introduce another 
category. 

Definition 14. Let the category of cone bundles £58nc)[ be the subcategory of 
*8nc)[ described as follows. An object C — > M o/£58nc)[, termed a cone bundle, is 
a smooth bundle for which there exists a vector bundle E ^ M and an inclusion 
bundle morphism l: C G E, such that each fiber Cx, x € M, is an open convex 
cone in the corresponding fiber E^ (where we have implicitly identified C with 
its image i(C) G E). Given two cone bundles C ^ M and C M' , with 
corresponding enveloping vector bundles E — > M and E' — >■ M' , for every £58nO[ 
morphism % : C — )■ C" there exists a vector bundle morphism ijj: E ^ E' such 
that X = "^Ic; namely the following diagram commutes 

C — "^E 

X V (29) 

C ^4 E' . 

Thus, F® -> F is obviously a cone bundle enveloped by {F*M)^ -> F. We 
refer to F® — > _F as the cone bundle of future oriented, spacelike covectors. 
Unfortunately, the closure F* = F® is in general not expected to be a manifold, 
since its boundary, consisting of the inner characteristic covectors C*, is fiberwise 
a piecewise algebraic variety and hence can have corners and other singularities. 
However, even though C*,F* G {T*M)° are not objects of £*8nJ)[ (by lack of 
convexity or by presence of a boundary) we refer to them as cone bundles anyway, 
namely the cone bundle of future oriented, inner characteristic covectors and the 
cone bundle of future oriented, cocausal covectors, respectively. 

Next we turn to causal and related vectors. These are most easily described 
using the following geometric notion of duality that is often used in convex 
geometry 1100!. 

Definition 15. Given a finite dimensional vector space V and a convex cone 
C G V , denote its closure by C and its open interior by C. We define the 
convex dual C* G V* as the set 



C* = {u(^V* \ u-v>Q for alive C}. 



(30) 
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We define the strict convex dual C® C V* as the set 

C® ^{ueV* \u-v>0 for alive C\{0}}. (31) 

To attribute strict may be dropped from the description of C® when it is clear 
from context. It is easy to check the following 

Proposition 2. Consider a convex cone C . 

(i) The convex dual C* is always closed and also convex. In addition, C** — C. 
The strict convex dual C® is always open and convex. 

(ii) C* \ {0} is non-empty iff C is contained in a closed half space. C® is non- 
empty iff C contains no affine line (it is salientj. 

(Hi) If C is open and salient, then C®® — C. 

(iv) The inclusion of cones Ci C C2 implies the reverse inclusion of their duals, 
CI D C2* andC® D C® . 

(v) The convex dual of the intersection of closures of cones Ci and C2 is the 
convex union ( convex hull of the union) of their duals, [Ci H C2)* = CI + C2, 
where the right hand side is written as a Minkowski sum, which for cones 
coincides with the convex hull of the union. The converse identity holds as 
well, (Ci + C2)* = C^ Co - Similarly, if Ci and C2 are open and salient, 
then (Ci n C2)® = C® + C® and (Ci + C2)® = C® D C® . 

For convenience, we extend this duality to cone bundles. That is, if C — ?> A/ is a 
cone sub-bundle of a vector bundle E — M, then the convex dual cone bundle 
C* — >■ M is the cone sub-bundle of E* ^ M such that each fiber C* is the 
convex dual of the corresponding fiber Cx, for x € M. Similarly, we can define 
the strict convex dual cone bundle C® — > M . The operations of intersection (n) 
and convex union are extended to cone bundles in the same way. Clearly 
both C M and C® M are cone bundles, that is objects of £*8ni)[. On the 
other hand, C* — ?> M may not be, though again we shall refer to it as a cone 
bundle anyway. 

Definition 16. The cone bundle of future directed, timelike vectors F C [TM]'^ 
is defined to be strict convex dual of the cone bundle of future directed, spacelike 
covectors, F = (F®)® . Let F C {TM)° be_called the cone bundle of future 
directed, causal vectors. Incidentally, F = {F®)* . Let C = OF = dF C {TM)° 
be called the cone bundle of future directed, outer ray vectors The past 
directed versions are obtained by reflecting the cones through the origin. 

While the definition of F® is primitive, it is also straightforward to show that 
in fact it is the strict convex dual of F, that is, F® = (F)®. 

These definitions extend straightforwardly to hyperbolic systems with con- 
straints. Suppose the adapted equation form (/ h,F* (Bm E*) defines the 
corresponding compound system, with (f,F*) the hyperbolic and {h,E*) the 
consistency subsystems. The compound system is quasilinear, but the consis- 
tency subsystem is linear in the constraints (sections of i? — > M). Therefore the 
cones of inner characteristic and spacelike covectors, as well as outer ray and 
timelike vectors depend only on the field bundle F and not on the constraint 
bundle E. Therefore, for a hyperbolic system on — )■ ill with hyperbolically 
integrable constraints, we define the cone bundles C* , F® and F over F as those 
associated to the corresponding compound system. 
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Note that these cones arc defined intrinsically (independently of any back- 
ground metric) by the geometry of a first order, quasilinear, symmetric hy- 
perbolic PDE system, or equivalently algebraic and geometric properties of its 
principal symbol. Of course, the dependence on a background metric (if one 
is present) may appear implicitly through the principal symbol, as it does for 
linear and semi-linear wave equations. This remark and the significance of the 
geometry of F and F® in the domain of dependence theorem to be described 
later is sufficient to motivate the fact that such cones are the central objects 
of the causal structure of classical field theory. The following definition leads to 
the study of conal manifolds and conal geometry alluded to in the introduction 
to this section. Many of the results in conal geometry do not depend on the 
provenance of these cones from the characteristic or spacelike covectors of a hy- 
perbolic PDE system. However, the later results concerning causal relations in 
PDE and classical field theories do use that connection. 

Definition 17. Given a manifold M, a chronal cone bundle on M is a cone 
bundle C ^ M (in the sense o/£*Bn5[j enveloped by the tangent bundle TM ^■ 
M such that each fiber Cx, x & M, is a proper cone (non-empty, open, convex, 
salient). The elements of C are future directed, timelike vectors. The strict 
convex dual C® — > M, enveloped by the cotangent bundle T*M — > M , is the 
corresponding spacelike cone bundle. The elements of C® are future oriented, 
spacelike covectors. The corresponding cone bundles C M and C® — )• M are 
referred to, respectively, as the causal and cocausal cone bundles. 

Note that the open convex dual of a proper cone is again a proper cone. Prototyp- 
ical examples of chronal and spacelike cone bundles are the puUbacks of the cone 
bundles of timelike vectors and spacelike covectors along a section (p: M ^ F, 
C = (f)*F and C® = (j)*F®. 

Definition 18. Let the category of chronal cone bundles C()r*BnZ)[ be the sub- 
category o/€*8nO[e where each object C ^ M must be enveloped by the tangent 
bundle TM — >■ M and the morphisms must be restrictions of the natural maps 

between tangent bundles induced by base space maps. Similarly, let the category 
of spacelike cone bundles ©p*Bnc)[ be the subcategory o/£*BnJ)[e where each ob- 
ject C® — )• M must be enveloped by the cotangent bundle T*M M and the 
morphisms must be restrictions of the natural maps between cotangent bundles 
induced by base space maps. 

Note that the morphisms are restricted to be open embeddings between base 
spaces to make sure that their natural extensions to the tangent and cotangent 
bundles are defined unproblematically. 

At this point, one may recall some standard notions of Lorentzian geometry, 
as long as they are defined only in terms of timelike cones, and apply them to 
the geometry cone bundles. Many of the standard theorems translate as well, 
some directly and others with some extra effort. We restrict ourselves to those 
that are relevant to the issues at hand. 

Fix a chronal cone bundle C — ^ M on a spacetime manifold M. Note that 
any chronal cone bundle is time oriented, since by definition the fibers consist of 
single cones rather than dotible cones like in the Lorentzian case. By assumptions 
the cones are directed into the future. 
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Definition 19. A smooth curve ^ in M is called future directed, timelike if 
the tangent to 7 is everywhere contained in C. The chronal precedence relation 
7+ C M X M (also ) is defined as 

= {{x,y) £ M X M \ 37, future directed, timelike curve from x to y}. (32) 

When {x, y) G /+ , we say that x chronologically precedes y and also write x <^y. 

If we replace C by C in the above definitions, we obtain causal curves and the 
causal precedence relation J+ C M x M, denoted x < y. The inverse relations 
are written I~ and J~ . 

It is immediate that /+ is an open, transitive relation {x <^y and y <^ z implies 
x <^ z). Using this relation, we can define the usual causal hierarchy. 

Definition 20. (i) If I'^ is irreflexive (x <^ x or, equivalently, no closed timelike 

curves exist), then C is chronological. 
(a) Given an open N C M, 7^|^ and I^r\{NxN) are both relations on NxN. We 
say that N is chronologically compatible if both of these relations coincide. 
More conventionally, this means that any two points of N that can be joined 
by a timelike curve in M can also be joined by a timelike curve in N. 
(Hi) An open N C M is called chronologically convex if any timelike curve that 
joins any two points of N must also lie in N, which is a stronger condition 
than chronological compatibility, 
(iv) The chronal cone bundle C is said to be strongly chronological if it is chrono- 
logical and for every x G M and every neighborhood N C M of x, there exists 
a smaller open neighborhood L C N that is chronologically convex. 

(v) The chronal cone bundle C is said to be stably chronological if there exists 
another chronal cone bundle C that is itself chronological and an open neigh- 
borhood of the closure of C , that is, C\{Q} C C . For spacelike cone bundles, 
stable chronology is equivalent to the reverse inclusion C® \ {0} C C® . The 
full subcategories l£[)r*Bn^ l^^. and Qp^ndlg^ o/£()r*Bnt)l and 6pQ5nt)( contain 
only stably chronological cone bundles as objects. 

Each of these definitions has an obvious analog when the adjectives chrono- 
logical or timelike are replaced by causal. 

Note that stably chronological is equivalent to stably causal, so these terms will 
be used interchangeably. Moreover, the chronological chronal cone bundle C 
such that contains a stably chronological chronal cone bundle C can itself be 
chosen to be stably chronological. Next we turn from curves to surfaces. 

Definition 21. Each of the following concepts may be prefaced with C- or C®- 
to be more specific. 

(i) An oriented codim-1 surface S G M is called future oriented, spacelike if its 

oriented conormals are everywhere contained in C® . 
(ii) A codim-1 surface S C M is called achronal if it {S x S) Ci I~^ = 0, that is, 

no two points of S are connected by a timelike curve. Similarly, S is acausal 

when {S X S)r\J+ = 0. 
(Hi) A codim-1 surface S C M is called Cauchy if it is acausal and every inexten- 

sible causal curve intersects S exactly once. 
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(iv) A chronal cone bundle C ^ M is called globally hyperbolic if there exists a 
Cauchy surface S Q M . A spacelike cone bundle C® M is called globally 
hyperbolic if C ^ M is. The full subcategories £t)c*Bn£)l^ and ©pSnOljif of 
£f)r!BnO[ and ©p*8nO[ contain only globally hyperbolic cone bundles as objects. 

Next we define some commonly used domains. Fix S* C Af to be a C-acausal 
codim-1 submanifold, such that either S is closed or S* C Af is a submanifold 
with boundary. 

Definition 22. (i) The future/past domain of influence I^{N) of a subset N C 
M is the set of points y S M such that there exists x ^ N with either x y 
(+) ory-^x (-). Let I{N) = I+{N) U /"(iV). 
(ii) The domain of dependence D(S) is the largest open subset of M for which S 
is a Cauchy surface. More commonly, D{S) is the set of points y G M such 
that every inextensible timelike curve through y intersects S. Let D^{S) = 
D{S)r\L^{S). 

(Hi) An open subset L <Z M is lens-shaped with respect to S if it can be smoothly 
factored as L ("1? 1) ^ with t: L ^ 1) denoting the projection onto 
the first factor (the time function j, such that the level set t = is S and all 
other level sets are spacelike as well as share the same boundary as S in M 
(which may be empty). 

The literature in relativity and Lorentzian geometry mostly makes use of 
the notion of global hyperbolicity as given above, but specialized to Lorentzian 
cone bundles. On the other hand, the literature on symmetric (and regular) 
hyperbolic PDE systems mostly makes use of lens-shaped domains. It is a non- 
trivial fact that these two notions coincide. The argument is essentially that the 
time function of a lens-shaped domain foliates it with Cauchy surfaces. Con- 
versely, a globally hyperbolic cone bundle admits a time function and a smooth 
factorization that turns it into a lens-shaped domain (glossing over some details 
related to spatial compactness). The original argument establishing the converse 
link in Lorentzian geometry is due to Geroch [IS]. However, his argument only 
established the existence of a continuous time function. The details necessary 
to establish the smooth version of the result are due to more recent work of 
Bernal and Sanchez [I8lll7j . The very recent result by Fathi and Siconolfi [43], 
using completely different methods, established the existence of a smooth time 
function and factorization for a class of cone bundles sufficient for the current 
discussion. 

Proposition 3. An open subset D C. M is C -lens-shaped with respect to S C D 
iff it is globally C'\d -hyperbolic, with S being C\D-Cauchy. 

The last point we address in this section is the kind of maps we allow between 
chronal cone bundles. Let £()rSnc)I(M) and (3p*Bn5[(M) denote the respective 
subcategories of £f)r!Bnc)[ and (3p*Bni)[ where the objects are restricted to have 
the same base manifold M and all morphisms must restrict to the identity on 
the base M . Essentially, the only allowed morphisms between the cone bundles 
in these categories are inclusions. These inclusions, and hence the morphisms, of 
the corresponding categories, implement a partial order on chronal and spacelike 
cone bundles, ordering the bundles by speed. Consider chronal cone bundles 
Ci — M and C2 — Af, as well as the corresponding spacelike cone bundles 
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Cf -s- M and Cf M. The morphism 
Ci >C2 




slow ■ 



->■ fast 



should read as 



id 



(33) 



Fas< and slow are used in the conventional sense: C2-timelike curves can propa- 
gate signals faster than Ci-timelike ones. This ordering is reversed for spacelike 
cone bundles. The morphism 




should read as 



fast 



M- 



slow 



id 



(34) 



. 



Definition 23. Given a chronal cone bundle C — > M , a section (f>: M ^ F is 
called C-slow (also slower than C) if (j>*r C C, that is, signals propagated by 
C -timelike curves can travel at least as fast as signals propagated by (j)-timelike 
curves. Also, (f) is said to be strictly slower than C if (ff F C C. Given a convex 
dual spacelike cone bundle C® — > M, the section 4> is (strictly) C®-slow (also 
(strictly) slower than C®) if it is (strictly) C-slow. 



F, 



In terms of spacelike cone bundles, C®-slowness is equivalent to C® C cj: 

while strict C® -slowness is equivalent to C® C (f)*F®. 

Consider two spacclike cone bundles C® C C® and two sections (f>i: M 
i = 1,2, such that (pi is C®-slow. Then clearly (f>2 is also C®-slow, while (pi need 
not be C®-slow. This observation shows that slow sections can be pulled back 
to slow sections along morphisms in the 6pQ5n9[(M) category. It is elaborated 
on in the next section. 

Let us expand our attention to the larger categories £()r*Bnt)( and ©p*Bnt)(, 
where the base manifold is no longer fixed, however only open embeddings are 
allowed as morphisms between base spaces (which naturally extend to morphisms 
of the total spaces). Consider a pair of objects and a morphism between them 
in the £()rQ5n()[ category. 

Tv 
Ci 



(35) 



Ml 



-4 Mo 



It is interesting to note that the same base space morphism x '■ M\ M2 extends 
to a morphism between the convex dual cone bundles in the Sp*BnZ)l category, 



(36) 



Ml 



->M 



2 , 
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only if the two cone bundles agree exactly, that is x*^i = ^"2, or equivalently 
Cf = X*Cf ■ In this case, it is clear that the relations Iq^ and on. the 

source manifold Mi are the same. However, it may not be true that the chrono- 
logical precedence relation /^^| ^ = xi^Ci) necessarily the same as the 

restriction I^^ n {x{Mi) x x{Mi)) on the target manifold Af2- The two are only 
equal when xi^i) is a C2-chronologically compatible subset of Af2. 

Definition 24. A £t)c*Bnc)l morphism Tx- Ci — > C2 in is called crhonologically 
compatible if x{Mi) is a C2- chronologically compatible subset of M2. Similarly, 
the morphism is called chronologically convex if x(Mi) is a C'2- chronologically 
convex subset of M2 ■ 

A ©pSnOl morphism T*x'- Cf C® is called chronologically compatible 
or convex if the corresponding £f)rf8nc)[ morphism Tx- Ci — ?> C2 exists and is 
respectively chronologically compatible or convex. 

When Ci — ^ M is a globally hyperbolic chronal cone bundle, there is no distinc- 
tion between and chronally convex C[)r*Bn5l morphisms. Thus, in this context, 
these terms will sometimes be used interchangeably. 

Lemma 2. Consider chronal cone bundles Ci — >■ Mi, i — 1,2, and an open 
embedding x- Mi — >■ AI2 such that the induced morphism Tx- Ci — >■ C2 in 
£[)r*BnI)[ is chronally compatible. Then Tx is also chronally convex. 

Proof. Consider any two points p,q (z Mi that are Ci-chronally related. Then, 
by chronal compatibility, x{p)tX{q) ^ M2 are also C2-chronally related. By a 
simple application of Zorn's lemma, among all timelike curves passing through 
xip) and xil) there must exist some maximal (inextensible) curves 72. For each 
such curve, the preimage 71 = x~^{'l2 H x(Mi)) is also inextensible. Assume 
for the moment that chronal convexity fails. Then we can chose 72 such that it 
has at least one point not contained in x{M), which implies that 71 consists of 
more than one connected component, each of which is inextensible. On the other 
hand, global hyperbolicity implies that each component of 71 intersects a Ci- 
Cauchy surface Ei C M and so the same can be said about 72 and 2J2 = xi^i)- 
In other words, 72 intersects IJ2 more than once, which is impossible, since 
from chronal compatibility M2 must be achronal. Therefore, Tx must also be 
chronally convex. □ 

These conditions on cone bundle morphisms will become important in the 
formulation of the generalized classical or quantum LCFT functors. In particu- 
lar, it would be extremely convenient for the construction of an LCFT, namely 
verifying the Isotony property (cf. Sect. 17. ip . if it were possible, in a spacetime 
with a globally hyperbolic chronal cone bundle, to be able to construct a Cauchy 
surface by extending an existing acausal surface, which may not be intersected 
by all inextensible causal curves. However, this question does not seem to have 
received much attention, even in the more restricted setting of Lorentzian ge- 
ometry. Therefore, we simply formulate this property as a plausible conjecture, 
whose validity would have to be explored in future work. 

Conjecture 1. Consider globally hyperbolic spacelike cone bundles C® — M and 
C® — )■ M', together with an embedding M — )■ M' that induces a chronologi- 
cally compatible morphism T*x'- C® C®. If T C M is a C®-Cauchy surface. 
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then for any compact subset K Q S, there exists a C"®-Cauchy surface S' C M' 
such that S' agrees with x(^) on x(?7), where J7 is a neighborhood of K in S. 

We say that the C"®-Cauchy surface S' extends K . If there exists a single C'®- 
Cauchy surface S' that extends every compact C we say that S' extends 
S. It is possible that there exists no S' that extends a given E. Think of a 
lens-shaped domain in M' whose spacelike boundary is not smooth. 



4-2. Slow sections. As we have noticed above, a distinctive feature of quasilinear 
PDE systems, compared to linear ones, is the field dependence of their causal 
structure, as embodied in the geometry of the cone bundles of timclikc vectors 
and spacelike covectors F and F® , whose base space is the field bundle F instead 
of just the manifold M. This means that any two sections (pjtp- M — >■ F of the 
field bundle (solution or not) define by pullback the cone bundles (j)*F M and 
tp*F — )■ M. A priori, these cone bundles need not be related to each other in 
any way. This complicates the analysis of the causal structure in contexts where 
multiple solutions or sections have to be considered simultaneously. In particular, 
such contexts arise when considering nontrivial subsets of the space of solution 
sections in classical field theory, probability distributions on the space of solution 
sections in classical statistical mechanics, and deformation quantization of the 
space of solution sections in quantum field theory. 

It is fortunate that multiple cone bundles can be easily combined via intersec- 
tion or convex union. Thus, if we have two sections 0, ijj : M — > F, we can form a 
chronal cone bundle from their convex union. C,/,.,^^ = (d)* F + -ip* F) . Then clearly 
both (j) and tp are C^^^-slow, or equivalently C<|i,^ is faster than both (j) and tjj- 
Similarly, in terms of spacelike cone bundles, both (p and tp are C|^^-slow, where 
C®^ = <p*r® n ip*r®. The only caveat is that the combined cones may fail to 
be proper, C^.^ could fail to be salient or C® ^ could be empty. In this case, we 
call the sections (p and ip chronologically incomparable. We can generalize this 
construction to any collection of sections {(pi}iQi. Define C/ = J^iei'^i^ 
Cf = (n,e/ (i)tF®)°. If cither of these cones fails to be proper, we again call the 
collection of sections {(j>i}i^i chronologically incomparable. 

This correspondence between sections and collectively faster cones can be 
easily inverted. 

Definition 25. Given a chronal cone bundle C — > M or the corresponding space- 
like cone bundle C® — >■ M , we denote the set of all sections strictly slower than 
C ^ M by T{F,C) C T{F), or also T{F,C®). We can similarly restrict the 
space of solutions, S{F,C) = S{F,C®) = S{F)r\V{F,C). It follows immediately 
from the definition that given two chronal cone bundles Ci C C2 (equivalently 
C® C C® or C\ slower than C2 ) we have the inclusions 

T{F,Ci) <ZT{F,C2) and S{F,Ci) C S{F,C2). (37) 

Incidentally, we denote by Tsc{F) C T{F) the space of stably causal sections 
and by Th{F) C TgciF) the space of all globally hyperbolic sections. The same 
notation also extends to solution spaces: Ssc{F) = <S(F) fl FgcC-F), andSuiF) = 
S{F)nTH{F). 
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li C ^ M is stably causal, then all sections in T{F, C) are a fortiori also stably 
causal. It then follows immediately from basic properties of stable causality that 

TsciF)^ U nF,C), (38) 

stably causal C 

where the union is over all stably causal cone bundles C M. It would be 
desirable to make a similar statement about TniF). However, that would entail 
showing that global hyperbolicity is stable in a certain technical sense. This is 
a known result for Lorentzian cone bundles [ L6. , but remains to be proven in 
the general case. However, since TniF) C TsdF) (and both are topologized 
by their inclusion in r(F)), the sets Th{F,C) = TniF) D r(F,C), where C is 
stably causal, do furnish a cover of TniF). We show below that each r{F,C) 
is open in T{F) in the Whitney fine topology. Thus, in this topology, T{F, C) 
and Th{F, C), with C stably causal, furnish open covers of Tsc{F) and Th{F), 
respectively. An observation that is important for later developments in Sect. [71 
While we could continue to work with an open cover by the sets Th {F, C) , it is 
simpler and likely not to make a large difference to instead assume the following 

Conjecture 2 (Stability of global hyperbolicity). The space of globally hyperbolic 
sections Th{F) is open in the space of all sections T{F), provided with the 
Whitney fine C*^ topology. Moreover, the following identity holds 

Th{F)= y r(F,C). (39) 

globally hyperbolic C 

This conjecture also implies the stability of global hyperbolicity in a common 
Frechet topology on T{F), which coincides with the Whitney fine C°° topology, 
as well, as it is finer (has more open sets) than the Whitney fine C° topology. 

Recall that Whitney fine C" topology (also known as the Whitney strong C° 
topology or the wholly open topology) on the space of sections of the bundle 
F — )■ M is defined as follows 62,75 . Consider sets of the form U{(t),V), with 
(/>: M — > Fa section and V C F an. open neighborhood of the image (j){M), 
consisting of all sections ip: M ^ F such that il}{M) C V. These sets form a 
sub-base for the open sets of the Whitney fine topology on T{F). This topology 
should be contrasted with the compact open topology (also known as the 
Whitney weak C° topology)., where a sub-base for the open sets consists of sets of 
the form U{4>, K, V), with 4>: M ^ F a section, K C M compact and V C F\k 
an open neighborhood of (j){K), consisting of all sections ip: M ^ F such that 
i'iK) C V. These two topologies coincide only if M is compact. 

The C°° versions of these topologies are easily defined in terms of jet prolon- 
gations. As before, consider the sets Uk{4>,V), but now with V C j'^F an open 
neighborhood of j''(j){M), and Ukifj), K,V), but now with V C J^F\k an open 
neighborhood of j''(j){K). They form sub-bases for the open sets of the Whitney 
fine C*^ and compact open topologies, respectively. The index k may also 
take on the value oo. With these topologies, of which the compact open one is 
somewhat more commonly used, T{F) is a Frechet space. 

Theorem 1. Given a chronal cone bundle C — M , the subset of smooth sections 
strictly slower than C, T(F, C) C T{F), is open in the Whitney fine topology. 
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Proof. Pick a section (j> e r{F,C). We will construct a Whitney fine C° neigh- 
borhood of it that is contained entirely in r{F,C). Every element ip of such 
a neighborhood would satisfy C® C ip*r® (recall that convex duality reverses 
inclusion) . 

The fibers of the cone bundle C® are cones of compact section. We can 
collect these sections into a continuous fibration K — > M with compact fibers as 
follows. Pick an affine sub-bundle r* C T*M whose fibers are of fiber codim-1 
and intersect the fibers of C® on compact sets with non-empty interior. Such an 
affine sub-bundle is easily constructed by modeling it on the vector sub-bundle of 
T*M that is annihilated by a section of C — )■ M, which always exists, since each 
fiber is contractible. Let K be the union of the boundaries of these fiberwise cross 
sections, K = 9(r* n C®). Essentially, K M is a. topological sphere bundle 
due to convexity of the fibers of C®. From the strict C-slowness hypotheses, we 
can conclude that 4>*r® is an open neighborhood of K. 

Next we pick special open neighborhoods for each point {x,u,p) G FxmT*M 
with u = (j){x) and {x,p) € K. Let f: M ^ S'^F* ®m TM be the adapted 
principal symbol of the symmetric hyperbolic PDE system under consideration. 
Consider the following composition of continuous bundle maps 

FxmT*M {S'^F* ®mTM)xmT*M S^D* , 

{y,v,q) ^ {y,f(y,v),q) ^ (y, f{y,v) ■ q) ■ 

We note the following facts: (i) the subset of S^F* consisting of orientation 
positive definite bilinear forms in the fibers is open, (ii) the above composition 
of continuous maps sends {x, u,p) in this open set, (iii) F Xm {(t>*r®) is an open 
neighborhood oi{x,u,p). These observation allow us to conclude that there exists 
a neighborhood U of (x, u,p) that maps via the composition (1401) to orientation 
positive definite bilinear forms in S'^F*, such that tj 'Z F Xm {4'*r®). Choosing 
a local bundle trivialization l contained within the base projection of U shows 
the existence of an open neighborhood of (x, u,p) of the form i{U xVx W) C U ^ 
where U is open in Af, V is adapted to the fibers of F and W is adapted to the 
fibers of T*M. This means that for any (y, v, q) G l{U xVx W) we have that q 
is {y, i')-spacelike and future oriented. 

Next, for each x G M, we let u — (j){x) and we build special neighborhoods of 
(x, u) € F and (x, u) Xm € F XmT*M . For each (x, u,p) € (x, u) Xm K^, let 
''x,p{Ux,p X yx,p X Wx,p) be the special neighborhood described in the previous 
paragraph. By compactness of the fiber K^, we can choose finitely many pj such 
that the corresponding neighborhoods cover (x, u) Xm K^- Denote by tt: F Xm 
T*M — >■ F the canonical projection. Then, let 

{Ux,p, X 14,p^. X W^^p^), (41) 

j 

W = 7r-\V) n |Jt,,p^.((7,,p^. X V,,p^ X W,,p^). (42) 
j 

That is, V is an open neighborhood of (x, u) G F and W is an open neighborhood 
of (x, u) X M Kx Q F Xm T*M. Moreover, by construction, for each (x, v) € V 
and each (x,p) G K^, the covector p is (x, w)-spacelike, since (x, v,p) G W. 
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Next, by continuity of the section : Af — > F and the bundle K — !• M, for each 
X G M, we can find an open neighborhood Ux of x, triviahzations bx of F and l'^ of 
T*M on Ux (with Lx = ttol'^), an open neighborhood Lx{Ux x Vx) C of (x, u) G 
F, and an open neighborhood i'x{Ux x 14 x 14^) C of (x, u) x m ii'a; C T*M, 
such that (j){Ux) C ia;(C/a; X 14) and K\if^ C (.^(C/a; x 14 x M^)- By construction, 
for every {y, v, q) £ l'x{Ux x 14 x Wx), and a fortiori for every (y, g) e -ft^lc/^, we 
have that q is (jj, u)-spacehke. 

Finally, let 

^f=\J ix{Ux X 14), (43) 

which is an open neighborhood of (t>{M) in F such that, for any other section 
tp: AI F with ■ip{M) C A/", each {x,p) G K is V'-spacelike. By properties of 
convex cones, this implies that C® C ip*r®. The set U{(j),Af ) of sections whose 
images are contained in Af is then a Whitney fine C° neighborhood of cj) that is 
contained in T{F, C). □ 

This theorem will be used in the following section when discussing the geometry 
and topology of the phase space of classical field theory. 

An analogous result for Lorentzian cone bundles was first given in |82) . The 
above proof is significantly more complicated because the fibers of both C® and 
(j)*r® cone bundles need not have elliptical cross sections, like they do in the 
Lorentzian case. 

4-3. PDF theory. As we shall see in Sect [5l PDE theory can be seen as a tool 
for a non-perturbative construction of classical theory. The following theorems 
are the main workhorses in this construction: local existence, possible global 
existence, as well as uniqueness, domain of dependence and finite propagation 
speed. The details of these results can be found in the standard literature on 
hyperbolic PDEs [6911791167] . See also more in depth discussion and bibliography 
in [M]. 

Theorem 2 (Local Existence). Consider a spacelike initial data set {S,ip). 
There exists an open neighborhood U of S in M , S d U ^ M , and a solution 
section (j): U F\jj that agrees with {S,ip), that is (l)\s — if on U . 

Recall that, for a linear PDE system, cone bundle F® does not depend on 
the dynamical fields. Therefore, it can be seen as a cone bundle over the space- 
time manifold M, F® — ?> M. For linear equations, the existence result can be 
strengthened to a global one [SSJ Ch.XIII], [75], [IHl Ch.IV], [Ml Ch.7], [7], 
p8, . 

Theorem 3 (Global Existence). Consider a spacelike initial data set {S,Lp) 
and suppose that the PDE system is linear. If the spacelike cone bundle F® is 
globally hyperbolic and the surface S is F® -Cauchy, then, under generic condi- 
tions, there exists a global solution section (p: M F that agrees with {S,ip), 
that is(f>\s — ^p. 

Theorem 4 (Uniqueness). Suppose that 4>: M ^ F is a solution section, that 
the spacelike cone bundle C® — (jf F® is globally hyperbolic and that S G M is 
a C® -Cauchy surface. Then (j) is the unique solution section on M that agrees 
with the initial data {S,(f>\s). 
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There is a more local, though equivalent, version of the uniqueness theorem, that 
usually goes under a different name. 

Corollary 1 (Domain of Dependence). Let {S,(p) be an initial data set and 
(p: U ^ F\u be solution section defined on an open neighborhood U Z) S such 
that (j)\s = ip. If V Q U is a lens-shaped domain with respect to S and the 
spacelike cone bundle (f)* F® , then (p is the unique solution section on V that 
agrees with {S,ip). 

Proof. By Prop. |31 the interior of V is globally hyperbolic and S is Cauchy with 
respect to (l)*r®. The result then follows directly from Thm. H) □ 

The same result can also be interpreted as showing finite speed of propagation 
of disturbances. 

Corollary 2 (Finite Propagation Speed). Suppose (f>: M ^ F is a solution 
section such that the spacelike cone bundle C® — (j)*F® is globally hyperbolic 
and the surface S C M is C® -Cauchy. If (j)' is another solution section whose 
restriction to S differs from (j) only on an open submanifold S' C S, then 4> and 
(j)' agree in the interior of the complement of the domain of influence Ic{S'). 

In other words, if a disturbance is confined to 5', it does not propagate faster 
than allowed by the chronal cones C = (j>*F oi the undisturbed solution. 

Proof. It is a consequence of the definitions that the interior of the complement 
M \ I{S') is the domain of dependence D{S"), where S" = {S \ S')° and the 
interior is taken with the respect to the submanifold topology on S. Since D{S") 
is globally hyperbolic, the result follows directly from Thm. □ 

What the preceding theorems allow us to do is parametrize globally hyperbolic 
solution sections in terms of initial data sets. The parametrization is unfortu- 
nately not bijective. Any solution section (p on M corresponds to many different 
initial data, one for each 0-Cauchy surface S, namely (SyCpls). On the other 
hand, not all initial data sets correspond to global solutions. There may exist 
data sets {S, (p) for which the corresponding solution, which is guaranteed to ex- 
ist in a neighborhood of S develops singularities within (p and so does not extend 
to a global solution. It is also possible that there exists a solution cp that agrees 
with the given initial data, but is not globally hyperbolic on M. Let us refer 
as globally hyperbolic initial data sets to those that come from the restriction 
of a globally hyperbolic solution section to a 0-Cauchy surface S C M. Two 
globally hyperbolic initial data sets that come from the same solution section 
are said to be equivalent, which defines an equivalence relation. Thus, by the 
above uniqueness theorem, globally hyperbolic solution sections are in bijective 
correspondence with equivalence classes of globally hyperbolic initial data. 

The last paragraph refers to solutions of hyperbolic PDE systems. What 
about those with hyperbolically integrable constraints? Consider the adapted 
equation form (/ © c, F* (Bm E), with (/, F*) and (c, E) respectively the hyper- 
bolic and constraints subsystems, and the corresponding hyperbolic compound 
system (/ h, F* (Bm E*) with identity h o c = q o f . Certainly, if a section 
(p: M ^ F satisfies both f[ip] =0 and c[(p] = 0, so does its restriction (or the 
restriction of j^(p for sufficiently high k) to any <^-Cauchy surface. On the other 
hand, consider initial data {S, ip) such that any extension of to a solution p 
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of = on an open neighborhood of S satisfies the constraints c[0]|5 = 0. Is 
it necessarily true that an extension of to a solution of /[</>] =0 on all of 
M also satisfies c[(l>] — on all of M? That is indeed the case if the section </> is 
globally hyperbolic on M with respect to the compound system, which includes 
the consistency subsystem {h, E*). To see this, notice that (j) © c[(/)] gives a solu- 
tion of the compound system on M, since the consistency identity and q[Q\ —Q 
yield 

h[cm = q[fm = <7[0] = 0. (44) 

On the other hand, we know that the section (/)©c[0] restricts to the initial data 
set (S'j ® 0) for the compound system. But the consistency subsystem is linear 
and zero is always a solution, implying that (/)®0 is also a solution with the same 
initial data. Finally, if the section (j) is in fact globally hyperbolic with respect 
to the compound system, the uniqueness Thm. 2] shows that the two solutions 
are identical, that is that c[0] = on M . So, for symmetric hyperbolic systems 
with symmetric hyperbolically integrable constraints, globally hyperbolic solu- 
tion sections 4>: M ^ F are in bijective correspondence with equivalence classes 
of globally hyperbolic initial data sets (S', ip) satisfying the constraints. Note that 
the notion of global hyperbolicity is in both cases with respect to the hyperbolic 
structure of the compound system (/ h, F* ©m E*). 



4.4- Linear inhomogeneous problems. In the preceding section we have discussed 
the Cauchy problem for quasilinear systems. If the system under consideration is 
linear, we can also discuss the linear algebra of the corresponding inhomogeneous 
problem 

m=a*, (45) 

where d* is a compactly supported dual density, that is, a section of the bundle 
F* M. 

The field independent spacelike cone bundle F® defined by the principal 
symbol /^^ endows M with causal structure. Let us assume that M is globally 
hyperbolic with respect to F®. It is convenient to introduce spaces of sections 
with restricted supports, in particular in ways related to the causal structure. 

Definition 26. Consider a vector bundle V M . We define the following sub- 
spaces of the space of sections T(V): 



To{V) = {(j)e r{V) I suppcj) IS compact}, (46) 

T+{V) = {0 e r{V) I suppcj) is retarded}, (47) 

r_(y) = r{V) \ supp cj) is advanced}, (48) 

^sc{y) = r(y) I snppcj) is spacelike compact}, (49) 



where retarded support, advanced support, or spacelike compact support means, 
respectively, that supp C I'^{K), supp0 C I^{K), or supp0 C I{K) for some 
compact K <Z M . 

The corresponding subspaces of the solution space S{F) are defined as 



So,±,sc{F) ^ S{F) nro,±,sciF). 



(50) 
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4-.4-1- Duhamel's principle. It is a well known fact of classical PDE theory that 
for some equations (like the heat and wave equations) the solutions of Cauchy 
and inhomogeneous problems are equivalent. That is, a solution of one problem 
can be obtained from the other. This is usually known as Duhamel's Princi- 
ple [311132]. For completeness, we present it here in a form appropriate for our 
geometric formulation of symmetric hyperbolic systems. In addition, once we 
have the retarded and advanced Green functions for the inhomogeneous prob- 
lem, we can construct an exact sequence that conveniently parametrizes the so- 
lution space of the linear PDE. This parametrization will be useful in Sects. [573751 
and 15.3.41 for the construction and analysis of the symplectic and Poisson struc- 
tures of classical field theory. 

Consider a Cauchy surface S C M (where the Cauchy property holds with 
respect to F® — M). Let M — M be a time function, such that t\s = and 
dt is future oriented. The ability to solve the Cauchy problem on S means that 
we have access to 

Definition 27. The Cauchy Green function G^' : To{F\s) ^sc{F) is a linear 
map uniquely defined by the requirement 

cj) = GE[ip] =^ /[0]-O and = (51) 

More explicitly, its action is given by the contraction of a bitensor distribution 
on. M X S with the initial data on S: 

Gs[vnx)= J^dSGsUx,s)^''{s), (52) 

where we have used local coordinates (a;^ u") on F and (s^ , v'^) on F\s, and ds 
is the coordinate volume (n — l)-form on S. 

Consider a compactly supported dual density a* G To{F*). The ability to 
solve the inhomogeneous problem f[(j)] = d* means that we have access to 

Definition 28. The retarded/advanced Green function G± : ro(-F'*) — >■ Tsc{F) 
is a linear map defined uniquely by the requirement 

4>± — G±[a*] fVf>±] = |5* o.'i^d supp0± C /^(suppci*), (53) 

where + denotes the retarded and — the advanced the boundary condition. 

More explicitly, its action is given by the contraction of a bitensor distribution 
on M X M with the source term, 

G±[a*r{x)= I Gt{x,y)al{y)dy, (54) 

J M 

where we have used the local coordinates (a;*, u°) and {y^ ,v^) on the two copies 
of F, and aj(y) dy are the components of d* in local coordinates, with dy the 
coordinate volume form on M. We also use the notation G{x\y) = G{x, y) for the 
Cauchy or retarded/ advanced Green functions, to highlight the different roles of 
the arguments. 

Consider any section of the field bundle 0: M — >■ and define 6± = 0{zLt), 
with 0{t) the Heaviside step function, which are the characteristic functions of 
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— I^{S), the future and past of U. The main identity that ahows us to 
relate the Cauchy and retarded/advanced Green functions is 

f[e±^] = {f-de±)4, + e±m. (55) 

A calculation in local coordinates (x^m") reveals d^^ = ±5{t)dt and 

[(/ . de±)q^]a - ±s{t) dufZ,^'' = ±s{t) p^,4>\ (56) 

where for convenience, we have defined /^^ by the equation dt A — dtif^^. 
We are now ready to prove what is known in classical PDE theory as 

Lemma 3 (Converse Duhamel's Principle). In local coordinates (x*,m°) on 
{t,s^,v'') on the two copies of F, with the second restricting to {s^,v^) on F\x:, 
the Cauchy Green function can be expressed in terms of the retarded/advanced 
Green functions as follows 

± 



where the pullhack is along the inclusion l: S (Z M, corresponding to t — 0. 

Proof. Consider ip e ro(i^|2;). We will construct a solution </> agreeing with 
this initial data, (f>\s — ^p, using retarded/advanced Green functions. If such a 
solution (j) existed, then from the identity (1551) it would satisfy 

f[9±cl,]a = ±S{t) dt A = ±S{t) dt A /*,^^ (58) 

where the S{t) prefactor allows us to consider only the values oi = ■ 
Conversely, knowing only the initial data on Z", we can define the following 
solutions to distributional inhomogeneous problems: 

/[0±] =5(t)diAf -0 =^ 0± -G±[(5(t)dtAf -^l. (59) 

Clearly, the section = 0+ — </>_ is then a solution of /[0] =0 that satisfies the 
identities = ±0± and hence, by Eq. ([SS]) . agrees with initial data = p^- 
Global hyperbolicity of M implies that </> is unique. Hence we can define the 
Cauchy Green function by G-s\^ = 0- More explicitly, in local coordinates 
(a;% w") and (t, s\ v^) on the two copies of F, we find 



r (x) = V ± / Gl^{x\t, s)6{t) dt A fl,{t, s)<p\t, s) (60) 





^±G55(x|0,s)/*,(0,s) 






. ± 





wherefrom the desired result follows immediately. □ 
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To obtain the retarded/ advanced Green functions from the Cauchy ones, we need 
only set up an initial value problem whose solution determines the retarded and 
advanced responses to the point source: /[V'i] — idFS{x, y). In local coordinates 
(a;%w°) and on the two copies of F, (idi?)^ = 6^ (Kronecker delta) and 

S{x, y) — Yl- 5{x^ — y^) Ax (Dirac delta densitized on the first argument). Without 
loss of generality, we introduce local coordinates {t,r^,v^) on F such that the 
point y lies on the surface l: S (Z M , y = (0, r) with S a Cauchy surface defined 
by f = 0. We define also = ('^*/a6)"^j (n - l)-multivector field on S, 
which means that 

i'^*fL)-t = fL-{^*t)=Sa: (62) 

where t* denotes the puUback of (n — l)-forms from M to S, l^, denotes the 
pushforward of (n— l)-multivectors from E to its image in M, and the • denotes 
the contraction of the (n — l)-form indices with the corresponding (n — 1)- 
multivector indices. We then have the direct form of 

Lemma 4 (Duhamel's Principle). In local coordinates (x%u°) = {t,s\u°') 
and (y^jV^) ~ (i,r^, v'') on the two copies of F , restricting to {s^jU"") and (r^, f'') 
on the corresponding copies of F\x;, the retarded/advanced Green function can 
be expressed in terms of the Cauchy Green function as follows: 

Gf{x,y) = ±e±ix)Gs:ix,r)df ■ ffir). (63) 

Proof. The result can be verified by direct calculation. First, note that when x — 
(0, s) the Cauchy Green function satisfies Gi;^(0, s\r) df = Yij ^i^'' ~ '""') 
Let 

r\x,r)^Gs:{x,r)df-ff{r). (64) 
It is clear that (/> is a solution of /[</)] =0 with initial data 

{cl,\sr\s, r) - n - • Lfi^)- (65) 

j 

If (j)f{x,r) = ±e±{x)(f>''''{x,r), then Eq. ^ implies 

f[cb±ra{x, r) - 5{t) dt A /Ib(x)/^(t, s\r) (66) 

= 6{t) dt A /l,(x) H S{s^ - r^) df • f;;{r) (67) 
j 

= Sit) n Sis^ - r^)dt A d~s {i*fl,){r) ■ r;{r) (68) 
j 

= 5l5{x\Q,r). (69) 

Since y = (0,r), we find the desired identity from Gj''(a;, y) = G'!^{x\Q,r) — 
(p'^ix, r). For other values of t we can find G'^{x\t, r) in a similar way. □ 
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4--4-^- Causal Green function (without constraints). Now that we are sure to 
have access to the retarded/advanced green functions G± for the hnear, sym- 
metric hyperbolic system f[(j>] — 0, we can define the so-called causal Green 
function 

G = G+-G_. (70) 

This new Green function helps to conveniently parametrize the space of solutions 
Ssc (F) — ker / C Tsc (F) by featuring in the following 

Proposition 4. The sequence 

> ro(^^) — U ro(F*) TsciF) — U TsciF*) > 0, (71) 

is exact (in the sense of linear algebra). 

The proof given in [7, Thm.3.4.7] and [107[ Lem.3.2.1] (which excludes the final 
surjection) for wave equations directly carries through to the symmetric hyper- 
bolic case. The final surjection is covered by the proof in [S] Ch.3,Cor.5]. 

We can interpret the above proposition in the following way. Since Ssc{F) = 
imG, we can express any solution to the homogeneous problem as (f) = G[a*], 
where a G ro(F*) is some smooth dual density of compact support. Also, since 
Tsc{F) ^ ini/, for any dual density a* with spatially compact support, there 
exists a solution (j) with spatially compact support of the inhomogeneous problem 
m = a*. 

Definition 29. Consider one Cauchy surface S <Z M and two more Cauchy 
surfaces C M to the past and future of S, where C and let 

— Let {x+tX-} be a partition of unity adapted to the open cover 

{S^ ,3"} of M , that is, x+ + X- — ^ ^"f^d supp x± C . We call {x+ , X- } ^ 
partition of unity adapted to the Cauchy surface S . 

Lemma 5. The exact sequence of Prop, ^splits at 

To{F*)=To{F)®Ssc{F) and T sc{F) Ssc{F) ®T sc{F*). (72) 

Given a partition of unity {x+tX-} adapted to a Cauchy surface S, there exist 
(noncanonical) splitting maps 

f^: imG -> ro(i?*), fM = ±/x*M = ±/[x±0], (73) 

G^:Tsc{F*)^Tsc[F), G^[a*] = G+[x+a*] + G^[x^a*]. (74) 

Proof. We demonstrate that the maps defined in Def . [251 are precisely the split- 
ting maps needed for this Lemma. Note that the splitting maps are not canonical, 
as they depend on the choice of a Cauchy surface and a partition of unity adapted 
to it. 

When (j) G Ssc{F) = imC, the identity [5S] shows that 

4^M=/[x±0]-(/-dx±)'/' (75) 

does in fact have compact support, as supp cj) is spacelike compact while supp dx± 
<Z r\ is timelike compact. Also, since d(x+ + X-) = 0, we have + 
f^Wi = 0; which means that the map = ±./^ is well defined. On the other 
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hand, we have G±[/[x±'/']] = X±0 from the uniqueness of solutions to the in- 
homogeneous problem with retarded/advanced support. The definition of the 
causal Green function then immediately implies that G o f-^ = ±id on SsciF)- 
Also, a direct calculation shows that / o = id on TsciF*)' 

foG^[a*]=x+a*+X-a*=a*. (76) 

This concludes the proof. □ 

We conclude this section by noting a simple but important fact. 

Lemma 6 (Covariance). Consider two manifolds M and M' with globally hy- 
perbolic, linear, symmetric hyperbolic PDE systems on them, with respective 
equation forms {f,F*) and {f',F'*). Suppose that the open embedding x- M ^ 
M' is such that {f,F*) — {x* f tX* F'*) o-nd, moreover, the induced morphism of 
chronal cone bundles is chronally compatible. Then the Causal Creen function 
of the two systems agree: G = (x x x)*G'. 

Proof. Since G = G+ — G_ , the result follows trivially if we know that G± = (x 
x)*G'j_. The latter identity follows from the uniqueness of the retarded/advanced 
Green functions on M and the fact that the pullbacks (x x x)*G'^ do in fact sat- 
isfy the desired Green function identity with appropriate boundary conditions, 
which is implied by chronal compatibility. □ 

4.4^.3. Causal Creen function (with constraints). The presence of constraints 
complicates the parametrization of solution spaces to the homogeneous and in- 
homogeneous problems. Recall that a symmetric hyperbolic system with hyper- 
bolically integrable constraints consists of the hyperbolic subsystem (/, F*), a 
constraints subsystem {c,E), and a consistency subsystem (h,E*) satisfying the 
identity h o c = q o h for some differential operator q. Since in this section we 
are concerned with linear systems, we take all differential operators to be linear. 
Moreover, the causal structure is presumed to be determined by the symbol of 
the linear symmetric hyperbolic compound system {f (B h, F* Q) E*). 

We will not discuss the most general kind of constraints and restrict our 
attention only to parametrizable ones. By the term parametrizable, we mean 
that there exist an additional vector bundle E' M and additional differential 
operators h' , c' and q' , which fit into the following commutative diagram 

r{E') ) T{F) — ^ r(£;) 

(77) 

r(£;'*) — ^ T{F*) — ^ T{E*) 

such that [h' ,E'*) is symmetric hyperbolic, and that the horizontal complexes 
of differential operators are formally exact (meaning that co c' = q o q' — Q 
and that the corresponding principal symbols form an exact sequence of vec- 
tor bundle maps), that is, they form a elliptic complex [66l §XIX.4]. Since 
both {h,E*) and {h',E'*) are symmetric hyperbolic, we can define their re- 
tarded/advanced Green functions, H± and H'j_, as well as their causal Green 



h' 
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functions, H = H+ — H_ and H' = — H'_. All these operators then fit into 
the following commutative diagram: 

> To{E') To{E'*) TsciE') TsciE'*) > 

/ ' / ' 

c q c q 

> ro(F) — U ro(F*) TsciE) TsciE*) > (78) 

c q c q 

> ToiE) ToiE*) TsciE) TsciE*) > 

We call the constraints subsystem globally parametrizahle if the elliptic complexes 
constituting the columns of the above diagram are all exact (their cohomologies 
are trivial). 

Lemma 7. The retarded/ advanced inhomogeneous problem 

m^p*, (79) 

40] = 7, (80) 
with /3* e ToiE*) andj£ T±iE), is solvable for (j) £ T±iE) iff h[-f] = q[i3*]. 
Proof. Let (p — G±[/3*]. We obviously have /[</)] ~ (3* . It remains to check 

= c[G+ m = H+ [q[p*]] ^ H+ [h[^]] =7. (81) 
This concludes the proof. □ 

4-4-4- Green functions and adjoints. We conclude this section by remarking the 
identities 

(G±)*-G;, (82) 

where on the left hand side iG±)* denotes the adjoint of the retarded/advanced 
Green function G± of the equation /[(/)] = 0, and on the right hand side G^ 
denotes the advanced/retarded Green function of the adjoint equation /* [(p] = 0. 
Note that taking the adjoint flips the support between retarded and advanced. 

Definition 30. Given two differential operators f,f*'. TiE) — > TiE*) are said 
to be mutually adjoint if there exists a bilinear differential operator G: TiE) x 
TiE) 17"^i(M) such that 

/[(/.]• - 0- =dG(0,V). (83) 

for any sections (j),ip: M — > E. The in — l)-form valued bilinear differential 
operator G(0, ■(/') is called a Green form associated to f and f* f89l ^IV.5], 
§ V.1.3]. Note that Eq. (j83p defines G up to the addition of an exact form W\5^. 
G ^ G + dH. Denote by [G] the uniquely defined equivalence class modulo exact 
local bilinear forms diJ(— , — ). 
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Notice that the principal symbols of mutually adjoint, first order, linear, sym- 
metric hyperbolic^ differential operators / and /* are negatives of each other 

f = (84) 

Also, the Green form of / and /*, given local coordinates on F, has the 

following representative 

G(0, V) = • (tr /) • ^|J = (tr Dab^^', (85) 

where (tr f)ab are (n— l)-forms obtained by contracting the single contravariant 
index of the symbol / with one of its covariant n-form indices. When pulled 
back to a codim-1 surface l: S C M, the Green form forms a density that can 
be integrated over E. Let E be defined as the zero set t = of a smooth function 
t. Then a straight forward coordinate calculation shows that 

i*G(0, V) - i*{trf)abr^' = i*fLr^', (86) 

where /^j, was defined in Eq. (j56p . In particular, when S is spacelike and future 
oriented, this shows that for a symmetric hyperbolic differential operator / the 
pulled back Green form t*G(— , — ) defines an orientation positive definite, sym- 
metric bilinear form on the fibers of the restricted field bundle F\s —i' S. This 
fact will be used in Sect. 15.2.51 

If we keep the orientation on M fixed, from the relation between their princi- 
pal symbols, the symmetric hyperbolic equation forms (/, F*) and (/*, F*) define 
the same spacelike covectors, except for the future/past orientation, which gets 
fiipped. This means that the causal relations defined by /* are simply the reverse 
of those defined by /. In this section, when using ± indices to denote retarded 
or advanced support, we always refer to the notions of past and future defined 

by /. 

Recall that we may introduce a natural pairing between elements cf) G r(i^) 
and a G T{F*) given by 




The pairing is only partially defined, that is, for simplicity, only for those pairs 
of sections for which the integrand 4> ■ a* has compact support. Its properties, 
including non-degeneracy, will discussed in more detail in Sects. I5.2.4I|5IT^ It 
is easy to show that the adjoint /* coincides with the adjoint of / with respect 
to this natural pairing: (/[</i],V') = {(j), f*[ip]) ■ This natural pairing allows us to 
define adjoints for integral operators like Green functions, namely (G± [a*] , (3*) = 
(a*, (G± )*[/?*]). 

It is straight forward to check that the identities foG± — G± o / = id hold on 
r±(F) as well as characterize the retarded/advanced Green functions, and also 



^ It is at this point that it be comes convenient to have chosen the adapted equation form of 
symmetric hyperboUc systems to be densitized. With the standard, non-densitized definition, 
the adjoint operator /* is in any case densitized. So it is only in symmetric hyperbolic form 
after contraction with a nowhere vanishing degree n multivector field. 
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that the natural pairing (— , — ) is non-degenerate on the spaces T±{F) x rzp(F*). 
It is now easy to verify the adjoint identities (|82t since 



JM JM 



G±o/[0±].5;, = (<^±,a;> - / o(G±)*[<5;], (89) 

M JM 

for any (j)± e T±{F) and G ^±{F)- The causal Green functions then satisfy 
(G)* = — G*, where G* is the causal Green function for /*. 



5. Construction of the Classical Field Theory 

Classical mechanical systems, and field theories in particular, have three stan- 
dard levels of description: spacetime, phase space, and observables. At the space- 
time level, the mechanical system is specified by the underlying spacetime as a 
manifold M, by the dynamical degrees of freedom as a field bundle F ^ M , and 
by the dynamics as a PDE system 8 C J''{F, M). At the phase space level, the 
set of all possible solutions forms a (possibly infinite dimensional) symplectic 
manifold, referred to as the phase space. At the level of observables, the me- 
chanical system is associated with the Poisson algebra of smooth functions on 
the phase space, known as the algebra of observables. For locally covariant field 
theories, all of these descriptions should be associated functorially to the given 
spacetime manifold M . A discussion of these functorial aspects is delayed until 
Sect. [SI For now, we discuss individual PDE systems. 

So far, we have only discussed field theories at the spacetime level, as PDE 
systems. To move on to the phase space level, we must consider the space of so- 
lutions of the PDE system. To become a viable phase space, it must be equipped 
with symplectic structure, preferably in a spacetime-local way (which is precised 
later on). It is well known that this is possible when the PDE system has an 
equation form of Euler-Lagrange equations of a local variational principle (or 
local action principle). The corresponding spacetime-local symplectic structure 
is given by a horizontally conserved symplectic current density on the jet space 
J°°F. If this symplectic current density is provided along with the PDE system, 
its space of solutions can be directly treated as the phase space of a classical field 
theory, without the need to introduce a local action. What is less well known [7T] , 
and perhaps a little surprising, is that a local action can be recovered from a 
symplectic current density. In other words, the two ways of specifying a classical 
field theory are essentially equivalent. 

The construction of the phase space of the classical field theory is broken up 
below into three sections. Sect. 15.11 starts with a local Lagrangian and extracts 
from it a local presymplectic form, which is used the construct the symplectic 
and Poisson tensors on the space of solutions in Sect. 15.31 Because the space of 
a PDE system is in general infinite dimensional, we first establish some formal 
properties of its tangent and cotangent spaces necessary for the discussion of 
these phase space structures in Sect. 15. 2[ which constitutes the bulk of this 
Section. Making the discussion the relevant infinite dimensional geometry non- 
formal would require substantially more functional analytical detail, which would 
detract from the geometric focus of this paper. More details in this direction can 
be found in [HllTlinillll] . 
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Before proceeding, it is worth remarking that the construction of the symplec- 
tic and Poisson structures on the classical phase space more commonly carried 
out in a Hamiltonian framework. However, such an approach usually requires a 
non-canonical 3+1 decomposition of an underlying 4 dimensional spacetime (and 
similarly in higher dimensions) , which destroys manifest 4-dimensional spacetime 
covariance. On the other hand, we hold spacetime covariance as an important 
guiding principle underlying the construction of locally covariant field theories, 
cf. Sects, m and [51 Fortunately, it has been known for a long time, that the 
symplectic and Poisson structures on the phase space can be built directly from 
the Lagrangian without giving up spacetime covariance |80ll38l[UD] . In fact, it is 
known from general principles that the usual Hamiltonian formalism is subsumed 
as a special case of this Lagrangian formalism pi?[|61[fiB] . 

5.1. Variational systems. Consider a field vector bundle F ^ M over an n- 
dimensional manifold M . A local action functional of order fc on F — >■ M is a 
function S[4)\ of sections (/): M ^ F, 



where £, the Lagrangian density, is a section of the bundle {A^MY J^F 
densities, which could depend on jet coordinates of order up to fc. The Lagrangian 
density is called local because, given a section (f) and local coordinates 
on J^F, the puUback at x e ill can be written as 



which depends only on x and on the derivatives of at a: up to order k. For the 
most part, the integral over M can be considered formal, since all the necessary 
properties will be derived from C. On the other hand, the finiteness of S[4)\ or 
related quantities may be important while discussing boundary conditions in 
spacetimes with non-compact spatial extent. However, we will no discuss these 
issues below. 

Recall that Sect, [^introduces the variational bicomplex Vl^''"{F) of vertically 
and horizontally graded differential forms on J°°F. Below, we use the notation 
introduced in that section. A Lagrangian density is then an element £ S n,"'^{F) 
that can be projected to J'^F. Incidentally the usual variational derivative of 
variational calculus can be put into direct correspondence with the vertical dif- 
ferential dv on this complex, which is how the name variational bicomplex was 
estabhshed [Il[2]. 

Let (x'^jUj) be a set of adapted coordinates on the oo-jet bundle J°°F, where 
all the following calculations can be lifted. Any result that depends only on jets 
of finite order can then be projected on to the appropriate finite dimensional jet 
bundle. Using the integration by parts identity (I300p if necessary, we can always 
write the first vertical variation of the Lagrangian density as 




(90) 



{f<j>r£ix)^£ix\djrix)), 



(91) 



dv/: = ELa A dvu'^ - dh9. 



(92) 



All terms proportional to dyW/, |/| > 0, have been absorbed into dh^. In the 
course of the performing the integrations by parts, ELq can acquire dependence 
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on jets up to order 2fc, and on jets up to order 2k — 1. Note that ELa = 
are the Euler- Lagrange equations associated with the action functional S'[(/)] or 
the Lagrangian density £. We can identify the form EL^ A dyu" with a bundle 
morphism EL: ,P'^F F* . Therefore, (EL,F*) is an equation form of a PDE 
system f el C J^'^F on F of order 2k. A PDE system with an equation form given 
by Euler-Lagrange equations of a Lagrangian density is said to be variational. 
Also, the form 6 is an element of fl'^~^'^{F), projectable to J'^^~^F. It is referred 
to as the presymplectic potential current density. Applying the vertical exterior 
differential to 9 we obtain the presymplectic current density (or the presymplectic 
current density defined by C if the extra precision is necessary). 

Lu = dv6i, (93) 

with w S r2"'~'^'^(F). This terminology implies that ui can be integrated over a 
codim-1 spacetime surface to construct a presymplectic form. Sect. [^75751 which 
is then necessarily local. This method of construction a symplectic form on the 
phase space of classical field theory is sometimes referred to as the covariant 
phase space method ^8(Jll35[ H] . 

The following lemma is an easy consequence of the definition of lu. 

Lemma 8. The form lu G il"^^''^{F) defined in Eq. (|93p is both horizontally and 
vertically closed when pulled back to Loo '■ £^ ^ J^^F : 

dhCw = 0, (94) 
dvC^ = 0. (95) 

Proof. The horizontal and vertical differentials on ^el are defined by pullback 
along Loo, that is, dht|^ — t^dh and dyt^^ = t^dv Since lu = dyO is already 
vertically closed on J^^~^F, it is a fortiori vertically closed on £. The rest is a 
consequence of the nilpotence and anti-commutativity of dh and A^: 

= dj£ = dvELa A dvu" - dvdh6', (96) 
AuLu = dhdv0 = -dvELa A dvu", (97) 
dhi;^w = lI^Aulu = -t^dvELa A d^u" = 0, (98) 

since EL^ and dvEL^ generate the differential ideal in Vl*{J°"F) annihilated by 
the pullback 6^. □ 

In fact, we will promote the name presymplectic current density to any form 
satisfying these properties. 

Definition 31. Given a PDE system l: £ C J^F we call a form lu a presymplec- 
tic current density compatible with £ if lu e [F) and it is both horizontally 
and vertically closed on solutions: 

dhC^ = 0, (99) 
dvC^ = 0- (100) 

The particular form lu defined by Eq. (|93p will be referred to as the presymplectic 
current density associated to or obtained from the Lagrangian density C, if there 
is any potential confusion. 
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5.2. Formal differential geometry of solution spaces. Before describing the sym- 
plectic and Poisson structures on the space of solutions, we should say something 
about the differential geometry of the manifold of solutions of a PDE system as 
well as its tangent and cotangent spaces. As usual for infinite dimensional man- 
ifolds, there are some subtleties. 

The main goal of this section is to describe the formal tangent and formal 
cotangent spaces of the manifold arbitrary field sections and the manifold of so- 
lution sections. The adjective formal, in the last sentence, alludes to the fact that 
we avoid most technical issues of infinite dimensional analysis and concentrate 
on what would be dense subspaces of the true tangent and cotangent spaces with 
a reasonable for their topologies. Results are algebraic and (finite dimensional) 
geometric identities that would form the core of an earnest functional analyti- 
cal formulation of their non-formal versions. The formal tangent and cotangent 
spaces have a natural dual pairing, which we prove to be non-degenerate, as 
a substitute for the absence of true topological duality between them. In the 
presence of constraints, the proof is carried out under some additional sufficient 
conditions. 

We start with Sect. 15.2.11 which discusses the choices of topology on the in- 
finite dimensional spaces of field configurations and solutions. Sects. 15.2.21 and 
Sect. 15.2.31 discuss sufficient conditions on the constraints and gauge transfor- 
mations needed for later results. Sects. I5.2.4[ 15.2.51 and 15.2.61 define the formal 
tangent and cotangent spaces in the progressively more complicated cases of the 
space of field configurations, the space of solutions (without constraints), and 
the space of solutions (with constraints). Finally, Sect. I5.27fl uses the preceding 
discussion to define formal local differential forms, of which the symplectic form 
will be an example. 

5.2.1. Choice of topology. What we really want to do is describe the space of 
solutions Sh{F). However, we first start with the space Th{F) or arbitrary field 
sections, because it has a simpler structure. In fact, since F — > M is a vector 
bundle, V{F) is a vector space. It can be turned into a topological vector space 
for several reasonable choices of topology. One such choice is the compact open 
topology (or more precisely the C°° compact open topology). When the base ma- 
nifold M is not compact, the Whitney topology (also known as the Whitney C°° 
topology^ or the wholly open topology) is another natural choice. Unfortunately, 
r(F) ceases to be a topological vector spac^fl, because multiplication by scalars 
fails to be continuous. These two topologies coincide iff the base manifold M 
is compact. The Whitney topology is naturally singled out by Thm. [1] which 
shows that the sets of slow sections r(F, C) are Whitney-operQ in Th{F). For 
the purposes of the discussion in Sect. |Bl it is advantageous to consider a cover 
of the space of solutions by open sets of slow sections. On the other hand, the 
differential geometry of infinite dimensional manifolds not modeled on a topolog- 
ical vector space becomes significantly more complicated. Also, it becomes more 
difficult to make contact with the approach of [171[551[^ , who use the compact 
open topology. 

^ With this choice of topology, C°°(M) is still a topological ring, but not a topological 
algebra over R. The space of sections T{F) is then a topological module over C°°{M), where 
both have the Whitney topology. 

^ They are proved to be open in the Whitney C" topology. They are also open in the Whitney 
C°° topology, since the latter is finer (it has more open sets). 
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One can think of different ways out of this impasse. One could bite the bullet 
and consider infinite dimensional manifolds modeled on topological C°°{M)- 
modules. One could also restrict the spacetime manifold M to be the interior 
of a compact manifold with boundary and restrict r(F) to the set of sections 
that extend continuously in some way to the boundary. At the moment, we 
remain agnostic about these or other possibilities, as they do not affect the 
results presented below, and simply assume that some choice has been made 
so that the sets r{F,C) are open in TniF), or could be effectively treated as 
such. For example, the choice of the functor C°°(— ) that assigns the algebra 
of smooth functions to our infinite dimensional manifolds is insensitive to the 
difference between the compact open and Whitney topologies. Therefore, we 
presume (assuming also that Conj. [5]holds) the following hypothesis 

Hypothesis 1. The space Th{F) of globally hyperbolic sections is open in T{F), 
and hence a manifold modeled on T{F), with an open cover by T(F,C), with 
globally hyperbolic C — > M, as in Eq. p9p . 

Now, the space of solutions can be topologized as a subset S{AI) C T{F). We 
are not interested in all solutions, only in the open subset of globally hyperbolic 
ones, Sh{M) = S{M) n Th{M). This will be our classical phase space. 

5.2.2. Constraints. Recall that we are considering a symmetric hyperbolic PDE 
system with constraints defined on the field bundle F ^ M , whose equation 
form is given by the hyperbolic subsystem {f,F*), the constraint subsystem 
{c,E) and the consistency subsystem {h,E*), satisfying the identity ho c = qo f 
for some differential operator q. The formal tangent vectors on the space of 
solutions will essentially consist of solutions of the linearized version of this 
PDE system. Supposing that some dynamical linearization point (p G Sh{M) 
is held fixed, the PDE subsystems linearized about (/) will consist of {f,F*), 
{c,E) and {h,E*), satisfying hoc — go/, respectively. All of /, c, /i, q are 
now linear differential operators. The essential properties of linear, symmetric 
hyperbolic systems with constraints that we will refer to are laid out in Sect. 14.41 
Dually, the formal cotangent space will be defined using the adjoint operators 
/* : r(F) T{F*), c* : r{E*) r{F*), h* : r{E) r{E*), q* : r{E) r(F), 
which satisfy c*oh* = f*oq*. Note that (/*, F*) and {h*,E*) are also symmetric 
hyperbolic. 

When dealing with constrained systems, some results covered in this section 
will require the further sufficient condition that the constraints be parametrizable 
so that we can extend both the linearized system and its adjoint to the following 
commutative diagrams: 

> ro(£;') To{E'*) Tsc{E') TsciE'*) > 

c' q' c' q' 
> To{F) — ^ ro(F*) TsciF) — ^ TsciF*) > (101) 

c q c q 

> ToiE) ToiE*) Tsc{E) Tsc{E*) > 
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(102) 



The rows form exact sequences, while the columns form elliptic complexes, as 
described in Sect. 14.4.51 Note that the adjoint diagram also describes a symmetric 
hyperbolic system with hyperbolically integrable constraints, except that the 
role of the constraint subsystem is now played by iq'*,E') and the consistency 
subsystem is {h'* , £"*), which satisfies the consistency identity h'* oq'* — c'* o f*. 

Since these systems are linear, their causal structures are field independent. 
Since we will be discussing both the linearized system and its adjoint system, we 
will refer to the causal structure defined by the extended compound systems (/i'© 
f®h,E'*(BMF*(BME*) and ih* ® j* ®h'\E* ®m F* ®m E'*). In fact, it is easy 
to check that their causal structures coincide and in turn coincide with that of 
the dynamical linearization point (f) £ SniM) and hence are globally hyperbolic. 
Recall that adjoint symmetric hyperbolic systems have opposite notions of future 
and past. We shall always take future and past to be defined by h'(Bf(Bh, rather 
than its adjoint. 

5.2.3. Gauge transformations. Many important classical field theories exhibit 
gauge invariance, like Maxwell theory, Yang-Mills theory, and GR. A gauge 
transformation is a family of maps g^: T^F) — >■ r(_F), parametrized by sec- 
tions s e r(P) of the gauge parameter bundle P — > M, that take solutions to 
solutions, while not changing modifying a field section outside the support of (5, 
gs[<p]{x) = (j){x) a X ^ suppe, which may be compact. If we linearize about some 
pair of background section 5 — ?► S+e, we obtain a linearized gauge transformation 
ge[(f\ gE[(f\ It is another requirement on gauge transformations that the 

generator of linearized gauge transformations g: r(P) — )> TiF) is a differential 
operator, which may depend on the background sections S and (f). 

Equivalence classes of sections under gauge transformations are considered 
physically equivalent. Therefore, physical observables will consist only of those 
functions on phase space that are gauge invariant (constant on orbits of gauge 
transformations) . Equivalently, observables are annihilated by the action of lin- 
earized gauge transformations. Another way to look at it, is to consider ob- 
servables as functions on the space of gauge orbits, denoted r^f(F) for the 
space of field configurations and iSi/(F) for the space of solutions. We call the 
space of globally hyperbolic solution sections modulo gauge transformations, 
SniE) — SniF)/^, the physical phase space. 

Often it is convenient to impose subsidiary conditions on field sections, called 
gauge fixing, that restrict the choice of representatives of gauge equivalence 
classes. The gauge fixing is called full if they only allow a unique representative 
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from each equivalence class, and otherwise called partial. The gauge transforma- 
tions that are compatible with a partial gauge fixing are called residual. 

Unfortunately, PDE systems with gauge invariance cannot have a well-posed 
initial value problem, and hence cannot be hyperbolized. However, the addition 
of subsidiary conditions on field sections can make the new PDE system equiva- 
lent to a hyperbolic one, usually with constraints. In practice, many hyperbolic 
systems with constraints arise after adding such gauge fixing conditions to a 
non-hyperbolic system with gauge invariance. However, there may remain non- 
trivial residual gauge freedom. For later convenience, as we did with constraints, 
we restrict our attention to what we call recognizable gauge transformations. 
That is, given linearized gauge transformations of the form g[e] and a partially 
gauge fixed symmetric hyperbolic (/, i^*), we can fit them into the following 
commutative diagram, whose columns form elliptic complexes: 



> ro(F) ro(p*) Tsc{P) rsciP*) > o 



> ro(^^) — ^-^ ro(F*) TsciF) — U r5c(^*) > o (103) 



>ro(P')^^ro(P" 



K' 



k' 



-40 



Their adjoints fit into the adjoint diagram whose columns are also elliptic com- 
plexes: 



< TsciP*) rsciP) ro(P*) ro(F) < o 



0^ 



TsciF*) TsciF) To{F*) ToiF) i (104) 



/* 



< TsciP'*) TsciP') ToiP'*) ^ ro(P') i 



The systems (fc, P*) and (A;'*, P'*) are required to be symmetric hyperbolic and 
P' — > M is called the gauge invariant field bundle, while g' is called the operator 
of gauge invariant field combinations. 

The above commutative diagrams are formally similar to those of symmetric 
hyperbolic systems with parametrizable constrains, as described in Sect. 14.4.31 
By analogy, we say that the gauge transformations are globally recognizable if all 
the vertical elliptic complexes in the two diagrams above are exact. 



5.2.4- Formal T and T* for configurations. Here we consider a section e 
TniF) and examine the formal tangent and cotangent spaces at 0, T^F = 
T^TniF) and r*F = TIF^ (P). 
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Definition 32. We define the formal full tangent space at as the set of space- 
like compact sections We define and the formal full cotangent space at (/) as the 
set 

T^r = rsc{F) and T^F = ro(F*). (105) 
The natural pairing (—,—): T^T x T^T — > R is 

{ip,a*) = [ V ■ a*. (106) 

Lemma 9. The natural pairing between T^F and T^T is non-degenerate. 

This is essentially the fundamental lemma of the calculus of variations and the 
proof is standard [33^ §IV.3.1]. 

Since the physical phase space will be identified with the space of gauge 
orbits Sh{F) in the solution space Sh{F), given a solution section </> S Sh{F), 
the formal tangent space T^S — T^Sh^F) at the corresponding equivalence 
class [<j)] £ Sh{F) in the space of gauge orbits consists of equivalence classes of 
linearized solutions up to linearized gauge transformations. Dually, the formal 
cotangent space T^S = T^Sh{F) will consist of dual densities annihilated by 
the adjoint of infinitesimal gauge transformation generator. 

Since gauge transformations act on field configurations and not just solutions, 
it makes sense to consider all field configurations related by gauge transforma- 
tions as physically equivalent. Thus, we also introduce Th{F) as the space of 
gauge orbits of field configurations as well as its formal tangent and cotangent 
spaces. The natural pairing between them is shown to be non-degenerate under 
the condition of global recognizability, that is, the vertical elliptic complexes in 
diagrams (|103p and (jl04p are exact. We deal with field configurations first and 
delay the discussion of solutions to the next section. 

The exactness of the composition g'og = ensures that we can recognize pure 
gauge field configurations, which are of the form ip — g[e] for some ^-spacelike 
compact section e : M ^ P, precisely as those (/)-spacelike compact field sections 
ip: M ^ F that give vanishing gauge invariant field combinations g'lip] = 0. 
On the other hand, the exactness of the dual composition g* o g'* = ensures 
that we can parametrize gauge invariant, compactly supported dual densities 
a* : M ^ F*, those satisfying g*[a] = 0, precisely as the image of the differential 
operator g'* acting on compactly supported sections of P'* — > M. 

Since gauge transformations act on arbitrary sections, not just solutions, we 
introduce the formal tangent and cotangent spaces for arbitrary and solution 
sections at the same time. As before, fix a section cf) in FniF) or Sh{F), as 
needed. 

Definition 33. The formal gauge invariant full tangent space at (p is the set of 

gauge equivalence classes of (p-spacelike compact sections, 

T^f = T^fH{F) = m I ^ e Tsc{F)}, (107) 
[V-] - ^ + .g[e], with e e VgciP). (108) 

The formal gauge invariant full cotangent space at (j) is the set of compactly 
supported gauge invariant dual densities, 

T^T = T;fH{F) - {a* e ro(F*) | g*[a] = 0}. (109) 
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The natural pairing (— , — ) : T^F x T^T — > K is 

m,a*) = f ^-a*. (110) 

Lemma 10. // the gauge transformations are globally recognizable, the natural 
pairing between the gauge invariant spaces T^T and T^T is non-degenerate. 

Proof. Non-degeneracy in the second argument follows once again from the fun- 
damental lemma of the calculus of variations: {[ip],a) = {ip,a) — for all 
i/j € T^r implies that a = 0. 

Non-degeneracy in the first argument requires an appeal to the global recog- 
nizability of the gauge transformations. Suppose that {[tpj^a) = for ah a S T^f . 
We need to show that this implies tp = g[e] is pure gauge, for some (/)-spacelike 
compactly supported £ : M ^ P. From the elliptic complex property of the first 
column of (|104p . we are free to use gauge invariant dual densities of the form 
a = g'*[e'*] for arbitrary e'* G ro(F'*) and find 

m,a*) = {^,a*) = = {9'm,n- (m) 

Since e'* could be arbitrary, the vanishing of {g'[ip],i'*) implies that — 0. 
On the other hand, from the exactness of the third column (I103p . we can then 
conclude that ip — g[e] for some e G ^sc{P)- This shows that ip is pure gauge 
and completes the proof. □ 

5.2.5. Formal T and T* for solutions (without constraints). The formal tangent 
space T^S will consist of linearized solutions, that is solutions of the linearized 
constrained hyperbolic system f[ip] — and c[ip] — 0. The formal cotangent 
space will naturally consist of equivalence classes of dual densities up to the 
images of the adjoints of / and c. After giving the precise definitions below, we 
prove that that the natural pairing between these formal tangent and cotangent 
spaces is non-degenerate, provided the constraints c[V'] = are trivial. The 
discussion of the case with non-trivial constraints is deferred to the next section. 

Definition 34. We define the formal solutions tangent space at as the set of 

spacelike compact linearized solution sections, 

T^S - T^Sh{F) = {V^ e Tsc{F) I m = m = 0}. (112) 

We define the formal solutions cotangent space at as the set of equivalence 
classes of compactly supported dual densities, 

t;s = t;sh{f) = {[a*] \ a* e ro(F*)}, (113) 

[a*]^a* +f*[i]+c*[e*], with ^eToiF), e* eVoiE*). (114) 
The natural pairing (—,—): T^S x TJS — > M is 



(V,[<5*]) = / V-a*. 
Jai 



(115) 
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Before proving the non-degeneracy of the above pairing, we introduce an 
auxihary pairing and prove its non-degeneracy first. Recah that the Green form 
(Def. [50)1 associated to / and /* is orientation positive definite when pulled back 
to a future oriented, spacelike surface, or at least the representative in Eq. (|86l) 
is. 

Definition 35. Let l: E C M be a future oriented, (j)-Cauchy surface and G 
the Green form associated to f. The Green pairing on S is a bilinear form 
(-, -)g,s ■■ ro{F\sf ^ M gwen by 

{^u^2)g,s^ [ i*G((^i,^2)- / >^*fL^>2- (116) 

J E Je 

Lemma 11. The Green pairing {—,—)g,e depends only on the equivalence class 
[G] and is positive definite (hence non-degenerate) . 

Proof. Any two representatives Gi and G2 of [G] will differ by an exact term 
dH, with H{—,—) a bilinear bidifferential operator. Therefore, the integrands 
i*Gi{(pi, ip2) will differ by the exact term di* H{ipi, (^2), with necessarily compact 
support. Therefore, since E has no boundary, we can use any representative of 
[G] to evaluate the pairing. 

As already mentioned above, the representative given by Eq. ((86|) is orienta- 
tion positive, L*G{ip, <p) > at s G i7, whenever ip{s) ^ 0. It then an elementary 
conclusion that 

{^,^)g.e^ f i*G((p,^)>0 (117) 
Je 

whenever ^ (p £ Tq{F\s)- Therefore, the Green paring on E is positive definite 
and non-degenerate, which concludes the proof. □ 

Lemma 12. // the constraints c[0] = are trivial, then the natural pairing be- 
tween T^S and T^S is non-degenerate. 

Proof. Non-degeneracy in the first argument follows again from the fundamental 
lemma of the calculus of variations: {if), [a*]) — for all a* € ro(-F*) implies 
that V = 0. 

Non-degeneracy in the second argument is more tricky. Suppose we have 
(i/), [a*]) = for all (^-spatially compact linearized solutions g T^S. From this 
we need to deduce that [a*] = [0], which means a* — f*[^] for some compactly 
supported ^ G Tq{F). Let K — suppa* and t+ : E^ C M a future oriented, 
0-Cauchy surface in the future of K and not intersecting it. Let 0- be the 
characteristic function of I^{E^), the past of E^ . If E^ is defined as a regular 
zero set t = of smooth function t on M , then d6'_ ~ —S{t) dt. 

Also, let ^+ = G;^[c5;*] be the retarded solution of f *[£,+] = a*, so then 
supp^+ C I^(K). Note that the intersection I~{E~^) nI'^{K) is compact and 
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contains K. In particular, we have O^a* = (5*, from which follows 

{^,[a*]} = {yj,0-a*)^ [ ^^-V-TK+l (118) 

JM 

?-(/M-e+-v^-rK+]) = - / e-dG(V',?+) (119) 

M JM 

d9^AG{tP,£,+ )^ - [ (5(i)di AG(V,C+) (120) 

M JM 

i;G(^,e+) = -(V^U+,e+b+)G,i;+- (121) 

Since ip is allowed to be any 0-spatially compactly supported linearized solu- 
tion, its restriction to the Cauchy surface S'^ could be any element of 
Tq{F\jj+). The non-degeneracy of the Green pairing, Lem.[Tl] then implies that 
the restriction is identically zero. 

On the other hand, by global hyperbolicity, we can foliate the future of S+ 
with other 0-Cauchy surfaces and apply the same argument to each of them. 
Then ^_|_ must vanish identically in the future of so its support is contained 
in the intersection D I'^{K), which we have already noted is compact. 

Finally, we have the desired conclusion that a* = f*[£,] with ^ = e ro{F). 
□ 

In the presence of gauge symmetries, the formal tangent space consists of 
equivalence classes of linearized solutions up to gauge transformations. On the 
other hand, the formal cotangent space is restricted to equivalence represented 
by gauge invariant dual densities. After giving the precise definitions below, 
we prove that the natural pairing between these formal tangent and cotangent 
spaces is non-degenerate, provided the constraints c[(j)] = are trivial and the 
gauge transformation are globally recognizable. The discussion of the case with 
non-trivial constraints is deferred to the next section. 

Definition 36. We define the formal gauge invariant solutions tangent space 
at (j) as the set of gauge equivalence classes of (p-spacelike compact linearized 
solution sections, 

T^S = T^Sh{F) = I ^ e Tsc{F),m = O,c[0] = 0}, 

[V'] ~ V' + .9[£], withe (^Tsc{P)- (122) 

The formal gauge invariant solutions cotangent space at <f> is the set of equiva- 
lence classes of compactly supported gauge invariant dual densities, 

t;s = t;sh{f) = {[«*] | &* e ro(F*),r[51 = o}, 

[<5*] - <3* + +c*[e*], (123) 

with ^ e Tq{f), e* e ro(l;*) (124) 

andg*[f*[S]+c*[e*]]=Q. (125) 
The natural pairing (—,—): T^S x T^S — > M is 

([^],[a*]) = / i^-a*. (126) 

JM 
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Lemma 13. If the constraints c[(j}\ = are trivial and the gauge transformations 
are globally recognizable (cf. diagrams (jlOSp and (I104p ). then the natural pairing 
between T^S and T^S is non-degenerate. 

Proof. Unfortunately, we cannot directly rely on the fundamental lemma of 
the calculus of variations to prove non-degeneracy in the second argument of 
([V'], [<5*]) = (V'l'i*), since we are only allowed to use equivalence classes rep- 
resented by gauge invariant dual densities, g*[a*] = 0. Fortunately, any dual 
density of the form a* = g'*[j3'*], with unrestricted choice of /3'* e Tq{P'*), 
is an allowed representative. Then the fundamental lemma of the calculus of 
variations and 

(M, [a*]) = [/3'1> = {g'm,n = 0, (127) 

for arbitrary /3'*, imply that g'[ip\ = 0. By global recognizability of gauge trans- 
formations, or more precisely the exactness of the third column of diagram (jl03|) . 
we have tp = g[e\ for some e S Tsc[P)- That is, ip must be pure gauge, [?/;] = [0], 
which proves non-degeneracy of the natural pairing in the second argument. 

To prove non-degeneracy in the first argument, we proceed as in the proof of 
Lem.[T21 Namely, assume that we have a gauge invariant dual density g* [a*] — 0, 
such that {[iP], [a*]) for all [iP] G T^S, and let ^+ = G;[a*], so that /*[^+] = 5*. 
Also, let K = supp/3'*, C M a future oriented, (/)-Cauchy surface that 
does not intersect K itself but does intersect I^(K), and 9^ the characteristic 
function of Then 

m, [a*]) = (V^,0-rK+]) - ~{Me+,^+\s+)g.e+, (128) 

which implies that supp^+ C and hence is compact. In other words, 

with € ^o{F), we have [a*] = [f*[£,]] = [0] and the that natural pairing 

is non-degenerate in the first argument as well. □ 

5.2.6. Formal T and T* for solutions (with constraints) . When the constraints 
are non-trivial, the proof that the natural pairing between formal tangent and 
cotangent spaces is non-degenerate complicates a bit. In fact, it requires the 
sufficient condition that the constraints be globally parametrizable. We first 
handle the case without gauge invariance. 

Lemma 14. // the constraints c[</)] = are globally parametrizable (cf. dia- 
grams (jlOip and (jl02p ), then the natural pairing between T^S and T^S is non- 
degenerate. 

Proof. Non-degeneracy in the first argument follows again from the fundamental 
lemma of the calculus of variations: (V', [<5*]) — {ip,a*) = for all a E ro(F*) 
implies that = 0. 

Non-degeneracy in the second argument means that {ijj, [a*]) = for all 
tp e T^S implies [a*] = [0], which means that g*[a*] = and a* = f*[^] -I- 
c*[e*] for some ^ € ro(^) and e* G ro(-E*). Using the same approach as in 
the unconstrained case, define ^+ = G!j_[Q!*], so that = a*. Also, let 

K — suppa and S'^ C M a future oriented, (/)-Cauchy surface that does not 
intersect K itself but does intersect I'^(K). Unfortunately, we will not be able 
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to show that vanishes outside a compact set, but instead we wih show that 
vanishes outside a compact set. 
To do that, we first need the identity 

(c'[V'']U+,e+U+)G,i:+ = (V''U+,rK+]b+>HM:+, (129) 

where H' is the Green form associated to the symmetric hyperbohc operator h' 
and ip' £ TsciE') satisfies — 0, hence also f oc'[^p] — 0. Let C" denote the 

Green forms associated to the operator c'. This identity then follows from 

dG(c'[^'],e+)=-c[V/]-r[e+] (130) 

= -^'.c'*[rK+]]-dc'(^',rK+]) (131) 

= ■h'*[q'*[^+]]-dC'{,p',a*) (132) 
= dff'(7A',g'*[e+])-dC'(^',fi*) (133) 

and the fact that supp C'{ip', a*) ^ supp a* does not intersect E'^. Since we are 
certainly allowed to use ip — c' [ip'] , using the same argument as in Lem. 1121 we 
have 

{c'mda*]) = -{c'm\s+,i+\s^)G,E+ ^ ~{^'\s+,q'*[M\s+)H',s+, (134) 

which forces (?'*[<^+] = on In fact, we can also conclude that the same 
equality holds on a neighborhood 5"+ of that contains /+(Z'+), on which we 
have f *[£,+] = as well. 

Now, let 7] denote the unique solution of /*[r?] = on M such that f] ~ £,+ 
on 5*+. Since, by the global parametrizability assumption (cf. diagram (|102p ') 
the constraint q'*[r]] = is hyperbolically integrable, it is satisfied everywhere in 
M as well. Therefore, from the results of Sect. 14.4.31 it follows that there exists 
e* £ ToiE*) such that 77 = q*[li*[e*]]. 

Let 77+ = g*[H^[e*]] and ^ = 1^+ — Then, by construction, ^ has compact 
support and 

5* = r[e+] = rK] + (135) 

= /* K] + /* [q* [h; [£*]]] = /* K] + c* [h* [h; r ]]] (i36) 

^rm+c*[e*], (137) 

as desired, implying [a*] = [0], which completes the argument that the natural 
pairing is non-degenerate in the second argument. □ 

The final non-degeneracy result is for the case where both non-trivial con- 
straints and gauge transformations are present. 

Lemma 15. // the constraints are globally parametrizahle (cf. diagrams (|10ip 
and (jl02p ) and the gauge transformations are globally recognizable ( cf. dia- 
grams p03p and (|104p ), then the natural pairing between T^S and T^S is non- 
degenerate. 

The proof follows straight forwardly by combining the methods of the preceding 
non-degeneracy lemmas, namely Lem. [T3] and Lem. [HI 
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5.2.7. Formal local differential forms. In this section, we use differential forms 
on J°°F to construct formal differential forms on r^f(F) or Sh{F). 

Consider any form a E rt'^''"{F) and suppose that (f>: M — > F is any section. 
We already know that the pullback along the jet prolongation j°°4' intertwines 
the horizontal differential dh with the de Rham differential on M, d{j°°(p)*a = 
(j°°</')*dha. However, if a has non-zero vertical degree, v ^ 0, then automatically 
{j°"(p)*a = from the defining properties of vertical forms. Though, if a is first 
contracted with v evolutionary vector fields tpj before effecting the pullback, the 
result 

(138) 

being the pullback of vertical degree form, is not-necessarily zero. Recall that, 
a section ip: M F defines a corresponding evolutionary vector field i/j that 
has the action ■ijj{x'^) = 0, ijj{uj) — diip'^ in local coordinates {x'\u'^). 

In fact, if we were concerned with calculus on Th(F) or Sh(F), the infinite 
dimensional manifolds of sections of F — >■ M and of solution sections of a PDE 
system, we could use vertical differential forms on J°°F to define local differential 
forms on the manifold of sections. A point local form A depends only on a specific 
point X G M. That is, evaluated at G Th{F), A depends only on (}){x) and 
finitely many of its derivatives at the same point. A local form is an integral over 
an M-family of point local forms, with respect to some smooth density on M. A 
degree v form on Th{F) or Sh{F) would be a functional of v tangent vectors ipj, 
which would correspond to sections of F — > M or linearized solution sections. 
Since we are not delving into the details of how these infinite dimensional objects 
arc defined, we are content to formally associate a local differential form ^jv[0] 
of degree v on Th{F) to each differential form a of degree {h,v) on J°°F as 
follows: 

An{^i, . . .,i^vM = I {r^rV^^ ■ ■ ■ i^a\, (139) 

J N 

where A'^ C M is a submanifold of dimension h. 

Recall that, given local coordinates [x^ ,u°-) on F, an evaluation functional is 
a function of Th{F) of the form <P°;.[d] = c6°(.t). A vector field 'F on Vh(,F) with 
{C^'F'^)[<j)\ depending only on 4){x) and finitely many derivatives at the same 
point is called a local vector field. We can formally associate a local vector ^ 
field to any evolutionary vector field '0 on J°°F, by defining its action on formal 
local 0-forms on r//(F) by 

{UAmM = I ir^'TCra. (140) 
Jn 

This is enough to define a formal exterior derivative 5 on Th{F), whose action 
on formal local forms can be formally given by the standard expressions with 
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contractions and Lie derivatives. It boils down to 



{SAn){>I'o 



(141) 




■■■L.a) (142) 



(143) 



where " means that the marked symbol is omitted. The last formula gives the 
real significance of the vertical differential. To derive the first equality, one must 
make use of the formal identification (|139l) and the identity (|299p . The second 
equality results from simplifications in the usual differential calculus on J°° F. 



5.3. Symplectic and Poisson structure. In this section, we endow the space of 
solutions Sh{F) of a variational PDE system with the structures of both a 
symplectic and a Poisson manifold (or rather formal versions of these structures), 
really turning it into the phase space of classical field theory. First, a variational 
system has to be put into symmetric hyperbolic forms, possibly with constraints. 
If the constraints are globally parametrizable, the gauge transformations globally 
recognizable and the hyperbolization of the variational Euler-Lagrange system 
satisfies an extra sufficient condition (which identifies part of the constraints 
with gauge fixing conditions), then we can use a generalized Peierls formula 
to construct the (formal) Poisson bivector and hence show that the covariant 
presymplectic structure is invcrtible and hence actually symplectic. 

5.3.1. Constraints and gauge fixing. Consider a Lagrangian density £ € f2"'°(i^) 
of order k and the corresponding Euler-Lagrange PDE system Sel C J^^F, 
with equation form (EL,F*). More often then not, the EL equation form is 
not directly an adapted equation form for a symmetric hyperbolic system. That 
could be true for several reasons. It is possible that ^el is in fact symmetric 
hyperbolic and that the EL equations need only be hyperbolized. However, if 
the EL equations are not of first order, then would need to be first reduced to 
first order form (this may be done already at the level of the Lagrangian, by 
introducing auxiliary fields). On the other hand, the EL equations may only 
admit a hyperbolization after several prolongations that identify all the relevant 
integrability conditions. Finally, if the Lagrangian is singular (it has non-trivial 
gauge invariance), the EL system cannot be hyperbolized at all, since systems 
with gauge invariance do not have a well-posed initial value problem. In that 
case, one would first have to introduce gauge fixing conditions, whose conjunction 
with the EL system is hyperbolizable (Sect. 13. Sp . 

Below, we postulate a condition on the hyperbolization of a gauge fixed vari- 
ational EL system that in the next section will be shown to be sufficient to 
construct certain Green functions needed for the Peierls formula for the Poisson 
bracket. By a hyperbolization we mean an equivalence between the gauge fixed 
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EL equations for (EL (B Cg, F* (B Eg) and a constrained symmetric hyperbolic 
system (/ (B c,F* Q) E). Here, Cg : T{F) — T{Eg) is a differential operator con- 
stituting the gauge fixing condition. The equivalence is presumed to have the 
following form 



EL = i? o (/ e c) 



Cg ^ RgOC 




(144) 



where the R, R, Rg and Rg are (possibly non-linear) differential operators. 

Since we will be concerned with symplectic and Poisson tensors, we can work 
with individual (formal) tangent and cotangent spaces at a fixed dynamical 
linearization point (p £ Sh{F), which is a globally hyperbolic solution of the 
PDE system defined by the equation form (/ (B c, F* (B E). In other words, we 
need to work with the linearized versions of each of the hyperbolic system, the 
constraints, the Euler-Lagrange system, the gauge fixing conditions, the gauge 
transformations, as well as the hyperbolization. As before, we denote the adapted 
equation form of the linearized hyperbolic system (/,_?*). The linearized con- 
straints are presumed to be globally parametrizable and fit into the commutative 
diagrams (jlOip and (|102p . The linearized gauge transformations are presumed to 
be globally recognizable and fit into the commutative diagrams (|103|) and (|104[) . 
The linearized EL equations are denoted (J,F*) and are also called the Jacobi 
system^ with J the Jacobi operator [38], while the linearized gauge fixing condi- 
tions are denoted by the equation form (cg. Eg). In local coordinates (a;*, u") on 
F, the components of the Jacobi operator satisfy the identity 

3i, A dvuj = dvELa. (145) 

The equivalence of the linearized systems takes the following form: 



r o f + Tc o c J f = ro3 + rcOCg 

TgO C [ C = f. I O J + O Cg 



(146) 



If the operator fj is non- vanishing, it means that part of the constraints consist 
of integrability conditions of the Jacobi system. 

Note that, strictly speaking, the r- and f- differential operators effecting the 
equivalence are not inverses of each other. Their compositions may differ from the 
identity by some differential operator that factors through a differential identity, 
that is, q o / — /i o c = or 5* o J = 0. In other words, we must have 

r o f + Tc o r J = id + pj o (7*, f o r = id + o (147) 

r o fc + fg — f o Tc + Tc ° fg = —pf o h, (148) 

rgorj^pgog*, fjor=PcOq, (149) 

'^g°^g ^ id, Tg o rg + fj o rc ^ id ~ pc o h, (150) 

for some differential operators pj, pf, Pg andpc- Also, the identity qo J —hoc — 0, 
when expressed in terms of the J and Cg operators, is identically satisfied when 

qof-hofj^qjog*, (151) 
q o fc — h o fg — 0. (152) 
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It is worth noting that the above relations involving the r- and r-operators follow 
from the equivalence (jl46|) only when J and Cg satisfy no additional differential 
identities. However, we will simply presume that they hold as needed sufficient 
conditions for the derivation of the Peierls formulas later in Sect. 15.3.41 

Finally, to make sure that the condition Cg ['0] = in fact constitutes a gauge 
fixing condition, we require the following compatibility between the gauge trans- 
formation operator and the constraints that we shall refer to as the gauge fixing 
compatibility condition: 

s' o fcO Cg = (153) 

This condition connects constraints (represented by Cg) and gauge transforma- 
tions (represented by s'). Essentially, this condition says that the part of the 
constraints €["0] = that comes from Cg[^] = is sufficient, when adjoined to 
J[-!/'] =0 to make the gauge fixed system hyperbolizable. 

To summarize, we require that the PDE system of interest satisfies the fol- 
lowing conditions: (a) variationality, (b) hyperbolizability, (c) global parametriz- 
ability of constraints, (d) global recognizability of gauge transformations, and (e) 
gauge fixing compatibility. The consequences of imposing these conditions are 
explored in the following section. We stress that these conditions are sufficient 
for our purposes and can in fact satisfied by most fundamental physical field 
theories, but some of the same results could also hold under weaker conditions. 

5.3.2. Causal Green functions. The goal of this section is to use the gauge fixing 
compatibility condition (jl53l) to construct a causal Green function for the gauge 
fixed Jacobi system. 

First, recall that the Jacobi system, due to its variational character is easily 
shown to be self-adjoint: 

J* = J. (154) 
Also, gauge invariance and Nother's second theorem imply the identities 

J o p = and o J = 0. (155) 

The equivalence of the gauge fixed Jacobi system and a constrained hyperbolic 
system postulated in (|146|) then gives 

Sok — fog^fcOCgOg. (156) 

Suppose that ip e Tsc{F) such that €(,[0] — and — g[e'] for some e' G 
^sc{P)- Then s o k[e'] = f ^ o Cg o g[e'] = 0, or k[e'] = (3*, where /3* e kers C 
Tsc{P*)- Let {x±} be a partition of unity adapted to a Cauchy surface (Def. 
and recall the associated splitting maps (Lem. [S|) , which lead to the identity 

± 

where s^[f3*] = ±s[x±f3*] G To{F*). Note that s'[s^[f3*]] = 0, which by global 
recognizability implies that there exists a 7* S ro{P*) such that s['y*] = s^[[3*]. 
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Note that 7* cannot be just x+P* or — because their supports need not 
be compact. Now, let e = e' — K^[/3*] + K[7*] and notice that k[e] = as well as 



= V. 



(158) 
(159) 
(160) 



We have just proven 



Lemma 16. // there exists tp € ^sc{F) such that Cg[ip] = and € im there 
exists £ e 7sc{P) such that k[e\ —Q and tf^ = g[e\. 

Equivalence with a constrained hyperbolic system now allows us to solve the 
inhomogencous problem 



(161) 



where the source must necessarily satisfy the gauge invariance condition g* [a*] — 
0. The equivalent inhomogeneous problem in symmetric hyperbolic form is 



m 



(162) 



Recall from Lem. [7]that this system is solvable iff the sources satisfy the consis- 
tency identity: 

q[f[a*]] - = pj o g*[a*] = o, (i63) 



which is obviously satisfied, after using identity (jl5ip . for any gauge invariant 
source. The retarded and advanced solutions to this inhomogeneous problem are 
then il}± — G±[r[a*]]. This means that = c[G±[f[Q;*]] and, in particu- 

lar, Cg[il}\ = 0. Motivated by this formula, we introduce the following retarded, 
advanced and causal Green functions for the gauge fixed Jacobi system: 



E± = G± o r and E = E+ - E_ = G o f. 



(164) 



One can immediately check that -0 — E[q!*] satisfies both /[■)/)] = and c[V'] ~ 0, 
whenever a* is a gauge invariant dual density. By the equivalence (I146p . the 
same solution also satisfies J['0] =0 and Cg[ip] = 0. 

Theorem 5. Provided the gauge fixed Jacobi system J[ip] — 0, Cglip] = is 
equivalent to a constrained hyperbolic system obeying the gauge fixing compati- 
bility condition (|153p . as described in Sect. \5.!^7J[ we can define the Jacobi causal 
Green function E, Eq. (|164p . such that the following diagram 



ro(P) 







TsciP) 







ro(^^) ro(F*) Tsc{F) Tsc{F*) > (165) 



ro(F*) 



rsciP*) 
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becomes an exact sequence after taking the vertical cohomologies. Moreover, we 
have the following splittings at ro(-F*) and Vsc{F)' 

Yerg* =mi3®Ssc{F) and coker 5 = 55c(i^) ® ker 5*, (166) 

where a (j)-Cauchy surface S d M and a partition of unity {x±} adapted to it 
define the splitting maps ( cf. Def. dP] and Lem. 0) 

3^:Ssc{F)^To{F*), J^[V.] = ±J[x±^] = r o /^[^] + o c^M, (167) 

E^:TsciP'*)^'^sc{F), E^[a*] = Y,G±[fW*[x±a*]], (168) 

± 

where Ssc{F) denotes the spacelike compactly supported solutions o//[0] = 0, 
dotc[(j>] — 0, and c^lip] — c[±x±V']- 

Proof. The proof consists of checking the desired conclusions at each object of 
horizontal sequence in the above diagram. 

1. Ifijje To{F) and J[^p] = 0, then ip = g[e] for some e G To{P). 
Let a* = /[V^]. From the equivalence (|146l) . we have 

a* = f o J[-0] + fc o Cg[il>] = fco Cg[ip]. (169) 

Then the gauge fixing compatibility condition (jl53p implies s'[a*] = 0. Now, 
note the identity tp = G+[q;*], which holds from the uniqueness of solutions 
with retarded support and from the fact that supp ip is automatically retarded 
by virtue of being compact. Then 

5'[^]=5'[G+[a*]]-KV[5'[<i*]]-0, (170) 

which, together with global recognizability of the gauge transformations al- 
lows us to conclude that -0 = g[£], for some e e ro(P). 

2. Every a* G ro(F*), such that E[a*] = g[e], with e e ro(P), and g*[a] = 0, 
can he written as a* — ^[^] for some tp G ro(-F). We also have imEoJ C img. 
The condition 

= G[f[&*]] = g[e] (171) 

implies that sofc[e] = fog[s] = 0, which in turn implies that g[e] — g[K[j*]] = 
G[s[7*]], for some 7* e ro(P*). We then have 

G[r[a*] - sir]] - 0, (172) 

which implies that r[a*] — s[7*] = /[V"], where ip e To{F) is unique and given 

by 

^ = G+[f[a*] - s[r]] - G_[r[5*] - sir]], (173) 
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since compact support is both retarded and advanced. Then the following 
calculation gives the desired conclusion: 

= r o f[ip] + Tc o c[iIj] (174) 
= r o f[a*] + Vc o CO G+ o f[a*] (175) 

~ r o s[7*] — 7'c o c o G-i- o s[7*] 
= (id + pj o g* — fj)[cy*] + r^o H+ o qo r[Q*] (176) 

- r o s[7*] - Tc o (fj o J + 7"g o Cg) o go K+[7*] 

= a* +rcO (-fj +H+ o /lo fj +H+ o qj o g*)[di*] (177) 

— r o s[7*] + r o (fc o Cg o g) o K+[7*] 

= a*, (178) 

where we have used the definition (|164l) and identities (|146D , (I147D , p5ip and 
the commutativity of diagrams (jlOip and (I103|) . Now, consider the following 
composition of operators applied to any ^ e ro(_F): 

5'[Eo J[?/;]] =.g'oGofo J[7/>] (179) 
= g' o (G o /)[V7] - K' o {s' o fe o (180) 
= 0, (181) 

where the last two terms vanish by the exactness of the causal Green func- 
tion sequence ([7T|) and the gauge fixing compatibility condition (I153P . Global 
recognizability then dictates that there exists a section e G ro(P) such that 
E o J[?/)] — g[e], which is the desired result. 
3. Every solution ^ e ^sc{F) of Jlip] = is of the form ijj = E[q;*] + g[e] for 
some e € Tsc{P) and a* G ro(F*), with g*[a*] = 0. 

Let S C M be a (/)-Cauchy surface and {x±} a partition of unity adapted to 
it. Let a* = ixi'^] ^'^'^ ^' = where now J['0'] = and €[-0'] = 0. Also, 

denote Cg;^[V'] = Cg[±x±'0]. It remains to show that '0 — -0' is pure gauge: 



g'[il;'] = g' o E[a*] = .g' o G o f[&*] (182) 

= .g' o G o f[r o /^[0] + o c^M] (183) 

= .9'°G[/xM+^corgOCx[V;]] (184) 

= (?'[^]+K'o(s'ofeOCgx)[0] (185) 

= 5'[0'], (186) 



where we have appealed to identities (|147p . (jl48p and the gauge fixing com- 
patibility condition (|153p . Now, we have (7'[-0— -0'] — 0, which implies, together 
with global recognizability, that V' — "0' = for some e € r5c'(-P)- 
4. Every a* € Tsc{F*), such that g*[a*] — 0, can be written as a* — J[0] for 
some ip € Tsc (-?")• 

Global recognizability implies that there exists a e'* e rs(7(P'*) such that 
a* = g'*[e'*]. Let = Ej(,[£'*]. Then the following calculation gives the desired 
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conclusion: 

m = 5](r o / + o c) o G± o f[g'*[x±e'*]] (187) 
± 

= rof[a*] +rcOf,-i[a*] = (id + p,j o g*)[5;*] (188) 
= a*. (189) 

5. We have E o J^^ = id (mod img) on SsciF)- 
Let ip G Ssc{F)- Direct calculation shows 

(?'[EoJ^[^]] =5'oGofoJ^[V'] (190) 
= 9' [Go f^] + H o (s' o o CgM (191) 
= g'[^] (192) 

where we have used a splitting identity from Lem. [S] and the gauge fixing 
compatibility condition p53p . It follows that — E o g im g is pure gauge. 

6. We have J o E^ = g'* on TsciP'*)- 

This was already checked in item H] above. □ 

It is easy to see from its variational nature that the Jacobi operator is self- 
adjoint (J)* — J. If it were directly invertible, the Green functions E± would 
satisfy the same relation with their adjoints as shown in Sect. (|4.4.4p . making 
the causal Green function anti-self-adjoint, E* = — E. However, due to gauge 
invariance the relation of the gauge fixed Green functions to their adjoints is 
more complicated. 

Lemma 17. When restricted to act on gauge invariant dual densities, the causal 
Green function of the gauge fixed Jacobi system is anti-self-adjoint up to gauge: 

(E)* = E (mod img). (193) 

Proof. First, note that from identities (jl46D and (I147p we have 

JoE±[<5*] =^(ro/ + rcOc)oG±of[5*] (194) 
± 

= r o f[a*] + rc o f][a*] (195) 
= (id+pjog*)[5*]. (196) 

It the follows that 

(E^)* o J o E± = (E^)* o (J o E±) = (E^)* + (E^)* o pj o 5*, (197) 
(Ezp)* o J o E± = ((Ezp)* o J*) o E± = E± + g o p* o E±, (198) 
and hence E± (Ezp)* + (Ezp)* o pj o 5* - g o p* o E±. (199) 

Given that E — E+ — E_, we then have 

E = -(E)*-5op*oE-(E)*op*og*, (200) 

which gives the desired conclusion. □ 



62 



Igor Khavkine 



As noted earlier in this section, there is considerable gauge freedom left that 
is compatible with the partial gauge fixing considered above. We call this kind 
of gauge fixing purely hyperbolic. Any further gauge fixing conditions arc then 
called residual. We leave the consideration of residual gauge fixing to future 
work. A principal difficulty in dealing with residual gauge fixing conditions is 
that the resulting constraints are no longer parametrizable (such as operators 
that are elliptic on a family of spatial slices). Thus, the kernel of the gauge 
fixing conditions may contain very few, if any solutions with spacelike compact 
support, which would be difficult to fit into the current formal framework for 
tangent and cotangent spaces to the space of solutions. 

5.3.3. Formal symplectic structure. In this section we define a formal presym- 
plectic structure on a slow patch S{F, C) of the solution space Sh{F); using the 
language of Sect. 15.2771 it consists of a formal local differential 2-forni fix: that is 
formally closed, 5i2s = 0. In a later section. Sect. 15.3.41 we shall see that some 
sufficient conditions ensure that is invertible and hence symplectic. Thus, 
the solution space becomes the phase space of classical field theory, as it can be 
endowed with both formal symplectic and Poisson structures. In the presence 
of gauge symmetries, this presymplectic form projects to an actual symplectic 
form on the space of gauge equivalence classes of solutions, that is, the physical 
phase space, or rather a slow patch S{F^ C) thereof. The symplectic forms will 
also be seen in a later section. Sect. 15.3.61 to agree on overlapping slow patches 
and hence constitute symplectic structure on the global phase space Sh{F). 

Fix a globally hyperbolic chronal cone bundle C — M , as well as a C-Cauchy 
surface S c M. 

Definition 37. Given a presymplectic potential current density co G 57"^^'^(i^) 
compatible with a PDF system l: £ C J^F and a closed codim-1 surface S C M , 
we can define the following local differential 2-form on the solution space S{F, C) : 



Recall that two formal tangent vectors £ T^S correspond to linearized 
solutions with spatially compact support, while V' and ^ are evolutionary vector 
fields on J°°F defined by them. 

Ideally, we would then prove that is a smooth differential 2-form on 

the infinite dimensional manifold S{F,C). However, that would necessitate a 
level of functional analytic detail going beyond the formal approach we have 
adopted. Instead, we note that the dependence of f2s[4>] on the solution section 
(j) is entirely through the ordinary differential form ui on the jet bundle J°°F 
(or even one of its finite dimensional projections J^F). Since w is smooth (by 
hypothesis), we also declare fls formally smooth. In addition to 0, the form 
i7i:[«?] depends on the surface S. For every e S{F,C), the supports of the 
linearized solutions and ^ (formal tangent vectors) will have compact intersec- 
tion with any surface that is 0-Cauchy. A fortiori, since S is 0-Cauchy for every 
e S{F, C), the integral defining flx:{ij,0[4'] converges. Moreover, because the 
presymplectic current density uj is horizontally closed on solutions, the value of 
the integral remains well defined and does not change as S is varied within the 
same homology class. 

It is now straightforward to prove the following 




(201) 
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Lemma 18. The formal local 2-form is formally closed, Sfix: = 0, and it 
depends only on the homology class of S. 

Proof. Fix a background solution cf) G Sh{F) and consider linearized solutions 
X, ^, '0 G T^S. It is straight forward to check that Qx: is vertically closed: 



{Sf2r)ix,i,m]= / (j°°0)*[i^t|t^dvH- / O-°°0)*C[^^^^^xdv^] (202) 



where we have used the vertical closure t^w, from the presymplectic hypothesis. 

Independence of the representative of the homology class of 2J follows if we 
can show that the integrand of i^s{x, £,) [0] = /i;O°°0)* ['-^'-x'^l Rham closed 
on M. This follows directly from the horizontal closure of i^w, again from the 
presymplectic hypothesis: 



d[(i°°0)*/,-i^H = (/»*dh/.*k-/.x^] = (jV)*'..-tx[dht*H = 0. (204) 



Note that, strictly speaking, on needs to apply Stokes' theorem to establish 
independence of the representative 2J, which requires the convergence of all in- 
termediate integrals. This property will actually hold in the main application of 
this result, which is Def. [551 below, so we do not discuss it in more detail here. 

In both cases, we have introduced the puUback l^^ because is a solution 
and are linearized solutions, so that the puUback (j°°0)* factors through 

Definition 38. Given a (gauge fixed, if necessary) variational PDE system on 
the field bundle F — )■ M, with Lagrangian density C, that is equivalent to a 
quasilinear, symmetric hyperbolic system on F ^ M with constraints and a 
globally hyperbolic chronal cone bundle C M , we define the formal variational 
presymplectic form on the space S{F,C) of C -slow solutions as Q — C2e, where 
we use the presymplectic current density u associated to C and S is C-Cauchy. 

Using the smooth splitting of M induced by a globally hyperbolic cone bundle, 
it is straight forward to show that all C-Cauchy surfaces belong to the same 
homologjG class. That, together with the fact that a C-Cauchy surface will 
compactly intersect the support of any 0-spacelike compact linearized solution, 
with (j) E S{F,C), shows that J7 is well defined, independently of the choice of 

A bilinear form defines a linear map from a vector space to its algebraic dual. 
A similar statement holds for a continuous bilinear form and to topological dual 
space. However, our formal cotangent space T^S is neither the algebraic nor the 
topological dual of the formal tangent space T^S. Thus we have to check this 
property for f2 by hand. This will be accomplished using an analog of a splitting 
map for the Jacobi system from Thm. [5l which is analogous Lem. [5] for purely 
hyperbolic systems. The argument in the proof was inspired by |46) . 

* The appropriate homology theory here should correspond to a variant of locally finite 
Borel-Moore homology, where one considers only chains whose intersection with every C- 
spatially compact set is compact. This variant does not appear to have gotten any attention 
in the literature and thus deserves further study. 




(203) 
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Lemma 19. The presymplectic form f2 defines the following map from the for- 
mal tangent space to the formal cotangent space: 

f2:T^S^T*S (205) 
tp^[a*], (206) 
with a* = = ±3[x±^], (207) 

where J: r(F) r{F*) the Jacobi differential operator at the dynamical lin- 
earization point (j) e Sh{F) and {x±} is a partition of unity adapted to a C- 
Cauchy surface S. 

Proof. Using the adapted partition of unity, we can write any solution of f[tp] = 
as tp = ip+ + V^_, with tp± = XiV*- If the solution also satisfies the purely 
hyperbolic gauge fixing condition s* o r*['ip] — 0, by the equivalence discussed 
earlier in this section it also satisfies the Jacobi equation = 0. Hence J[V'+ + 
ip-] = or J[^+] = — J[^/^_] = Jx[V']- Note that the support of J[tp±] is compact, 
since V'± satisfy the Jacobi equation away from the intersection S~^r\S~ flsupp tp, 
which is by hypothesis compact. 

Next, we want to find a dual density a* that satisfies f2[(j)]{^,ip) = {^,a*) 
for any ^ € T^S, which in particular satisfies J[^] = 0. Recall that an adapted 
partition of unity also depends on two additional C-Cauchy surfaces S"^ C 
/^(Z") and the supports of the partition are contained in suppx± 

mi^)[<l>]= I w(e,v)[0] = E / ^i^^^±M (208) 

Je _|_ Je 

= E/ w[<A](e,v±) + E± / MCM[<t>] (209) 

_l_ JET _|_ Js±nlT(E) 

= V±/ (dha;)(e,V±)['^] (210) 
{-d,ELaAd,u''){t^±M (211) 



IT(E) 



E^/ [{Ji,djem-{3Udjyji}a (212) 

e-J[V'+]-/ ^-Jl^-] (213) 

I~{E) Jl+(E) 



= [ e-JxM-/ ^(-J> 

Jl-(E) Jl+{E) 
J A 



(214) 



e-JxM (215) 



X 

M 



After the integration by parts, the boundary integrals over were dropped 
since they did not intersect the support of their integrands. Since the supp'^± C 
S^, the integration over 5=^ n P{S) could be extended to all of I'^iS). The 
term ip± ■ J[^] was dropped since ^ is a linearized solution. 
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Recall that we are not interested in the dual density a* = JxtV"] specifi- 
cally, which explicitly depends on the adapted partition of unity, but rather the 
equivalence class [a*] € T^S, defined modulo imJ. Consider another adapted 
partition of unity {x±}- Because they both provide splitting maps, we have 
^j; = E o 3 ^[tp] = Eo J^/[-0]. Then E[J^'[-0] — J^[V']] = 0. By exactness, and 
3y [ip] must differ by an element of im J; in other words, they represent the same 
equivalence class in T^S. 

To complete the proof, we use the non-degeneracy of the natural pairing 
between T^S and T^S to define the operator f? by the formula 

(e, 12^) = m (e, ^) = (c, «*) = (c, [5*]) , (216) 

so that 12-0 = [a*] G T^S, with a* = J^lip]. □ 

Notice that in the presence of gauge symmetries (the gauge fixing is only 
partial), the form i? is degenerate, since every pure gauge solution lies in its 
kernel: 

n{g[e],^M = (5[£],±J[x±V]) = (e,±r° JfXiV-]) =0 (217) 
for any ip, since Nother's second theorem implies that g* o J = [50] , 

Corollary 3. The 2- form Q on T^S projects to a 2- form D on T^S and hence 
defines a map 

f2:T^S^T;S (218) 

m ^ (219) 

with a* = J^lip]. (220) 

In other words, formally, the quotient projection to the physical phase space 
effects a presymplectic reduction {Sh{F), fl) — ?> {Sh{F), We shaU see later 
on that Q is non-degenerate and hence symplectic. 

5.3.4. Formal Poisson hivector, Peierls formula. In this section we show that 
the formal symplectic form Q defined above is invertible and that its inverse, 
the formal Poisson hivector 77, is given by the Peierls formula 

77 = E, (221) 

where E is the causal Green function of the Jacobi operator J as defined in 
Sect. 15.3.11 To show that 77 is indeed a Poisson bivector, it suffices to show that 
(a) it is an antisymmetric bilinear form on the formal cotangent space, (b) it 
defines a map from the formal cotangent space to the formal tangent space and 
(c) it is a two-sided inverse of Q defined in Cor. |31 The fact that 77 defines 
a Poisson bracket, with its Leibniz and Jacobi identities, then formally follows 
from standard arguments. 

Lemma 20. The Peierls formula specifies a map from the formal cotangent 



space to the formal tangent space: 

n-.T^S^T^S (222) 

[a*]^[iP\, with g*[a*]=Q (223) 

and V = E[a*]. (224) 
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Proof. The challenge is to show that 77 maps equivalence classes to equivalence 
classes (Def. [36]). That is, that any representative a* +c*[7*] of an equiv- 

alence class [a*] e T^S, with = 0, gets mapped to the same equivalence 

class in T^iS. By linearity, it suffices to check that [0] £ T^S is mapped to 
[0] e T^S. Recall that any solution representing [0] G T^S is pure gauge g[e]. 
Note that the equivalence of the {f(Bc,F*(BE) and {J®Cg,F*®Eg) equa- 
tion forms, together with the self-adjointness of the Jacobi operator (J)* = J, 
allows us to rewrite any representative of [0] G T^S as J[^] -|- c*[7*], for some 

C € ro(^) and 7* € To (7?*). This representative wih also satisfy the identity 

r°c;[7l=r[JK]+c; [711= 0. (225) 
Direct calculation then shows that 

mm + c;[r]] = e o m + e o cur] (226) 

= g[e]- (cg o E)*[7*] - .g o o E o d^r] (227) 

= g[e - q*j o (H)* o r;[r] ~P*joEo c;[r]] (228) 

is pure gauge. We have used the identity that E o J[^] = g[e] for some e G ro(P) 
(Thm. [S]), the anti-self-adjointness identity (^00]) . that (Eqs. (fUSl) and (fT5T|) ) 

Cg o E = Tg o (c o G) o f = o (H o q o r) (229) 
= {rgoRoqj)og* (230) 

and the identity (P^ . 

Therefore, we can conclude that if [a*] — [0], then [E[q;*]] = [0]. □ 

Lemma 21. The Peierls formula defines an antisymmetric bilinear form on the 



formal cotangent space: 

n{[a*], ri) = (7T[ci1, [/31) = -77([r], [«*]), (231) 
for any [a*],[^*] e T;S. 

Proof. Recall that the representatives always satisfy (7*[a*] = = 0. Ap- 

pealing directly to the anti-self-adjointness identity (j200p we have 

7T([5*],r]) = (77[5*],r]) = (E[a*],/3*) - ((E)*r],5*) (232) 

= -((E[r] + g o p* o Er] + (E)* o p, o r[r]), «*) (233) 

= -([Er]],[a*]) = -(77[r],[a*]) (234) 

= -7T(r],[a*]). □ (235) 



Theorem 6. The Peierls formula gives a two-sided inverse to the formal sym- 
plectic form, fin = id on T^S and II f2 = id onT^S . 

The argument in the proof below was inspired by |46| and |107l Lem.3.2.1]. 
However, the argument has been generalized to handle hyperbolic systems with 
constraints and gauge invariance (presuming purely hyperbolic gauge fixing). 
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Proof. The proof uses in an essential way the splitting identities of Thm. [SJ 
Consider any [?/;] £ T^pS and [a*] £ T^S. To use these splitting identities, we 
introduce a (/)-Cauchy surface S C M and a partition of unity {x±} adapted to 
it. Then 

{nn[^],[&*]) = {Eo3^[ij],a*} (236) 
^{^ + 9[e],a*) (for some e G TsciP)) (237) 
= (238) 

Therefore, from the non-degeneracy of the natural pairing between T^S and 
T^S, we concluded that USl — id. Similarly, 

m,nn[a*]) = {^,3^oE[a*]). (239) 

But then 

E[Jx o E[<5*]] = (E o J^) o E[a*] = E[&*] + g[e], (240) 

for some e G ro(P). But, by the cohomological exactness of Thm. [SJ this means 
that J;^ o E[q;*] — a* = J[i^] for some ^ G ^o{F). In other words, 

{[^],nn[a*]) = {^j,a* + m) = (MJ^*]). (24i) 

Therefore, from the non-degeneracy of the natural pairing between T^iS and 
T^S, we concluded that f2n = id. □ 

5.3.5. Algebra of observables, locality. This section and the next complete the 
construction of a classical field theory by constructing, at least formally, the 
Poisson algebra of observables associated with it. Along the way, several impor- 
tant points are discussed, (i) Since the focus of this work is more geometric than 
analytical, some parts of the construction are kept formal, (ii) In the physics lit- 
erature, the given construction is often called on-shell. We discuss its relation to 
and trade offs with respect to the off-shell formalism, which has received consid- 
erable attention in the recent literature on perturbative algebraic quantum field 
theory (pAQFT). (iii) We define the notion of spacetime support, corresponding 
to the spacetime localization of observables, and use it to prove a classical ver- 
sion of microcausality that is generalized to field theories with dynamical causal 
structure. 

We will use the Peierls formula for the Poisson bivector U to construct a 
formal Poisson bracket on the formal algebra of observables on the phase space 
{Sh{F), ST) of a classical field theory. This section will actually concentrate on 
a single slow patch S{F,C) C Sh{F) of full phase space. Globalization to the 
full phase space Sh{F) by covering it with slow patches will be considered in 
the following section. Sect. 15.3.61 

In the last paragraph, the word formal was used in two different ways. The 
algebra that we will consider is the polynomial algebra generated by local func- 
tionals with compact spacetime support. This algebra is much smaller than a 
reasonable space of smooth functions on the infinite dimensional manifold of 
solutions SniF), and hence is only a formal substitute. However, this polyno- 
mial algebra is much less complicated from the analytical point of view and 
is expected to be dense in the larger, ultimately desired algebra C°^{SH{Fj), 
however it is rigorously defined. Also, the Poisson bracket of two elements of 
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this polynomial algebra, if defined using the bivector 77, is not an element of the 
same algebra, though still a well defined function on Sh{F). However, just as the 
polynomial algebra is expected to be dense in C°°(5_f/(F)), the Poisson bracket 
is expected to be continuous, so nested applications of the Poisson bracket can be 
considered by replacing each evaluation thereof by an approximating sequence 
within the polynomial algebra. In this sense, the definition of the Poisson bracket 
is given only formally. The hypotheses needed to justify its non-formal use are 
summarized in Hyp. [51 

Before moving on to the technical part of this section, briefly discuss the 
relation between our approach to the algebra of observables and the one taken 
in the recent literature on pAQFT ,64..23 „47.,95..24j . That work is in the so-called 
ojf-shell formalism, while we work in the on-shell one. The space of all globally 
hyperbolic solution sections SniF) C TniF) is a submanifolcfl of the space of all 
globally hyperbolic sections and is sometimes referred to as the shell. Thus, on- 
shell field configurations correspond to solutions, while the off-shell ones are not 
so restricted. Dually, on the algebraic side, the algebra of smooth functions on 
Sh{F) is a quotient C°°(cSjy(f )) ^ C°° {T h{F)) / h by the ideal h generated by 
the equations of motion (an equation form of the PDE system J^F D £" — M). 
Correspondingly, C°°{Sh{F)) is called the on-shell algebra while C°°(r/f(i^)) 
is called the off-shell algebra. The on-shell formalism attempts to construct the 
on-shell algebra directly as the algebra of functions on the space of solutions (or 
gauge equivalence classes thereof), while the ofF-shcll formalism rather constructs 
it as the above quotient of the off-shell algebra (with a further quotient to factor 
out gauge). 

The off-shell formalism is advantageous in perturbative algebraic quantum 
field theory (pAQFT) as it dramatically simplifies the perturbative renormal- 
ization of interacting theories (64.,23.,47.„95j . On the other hand, the on-shell 
formalism connects more directly with the mathematical PDE literature. It is 
the theorems on the well-posedness of hyperbolic PDE systems, as quoted in 
Sec. 14.31 that will allow us to establish the validity of the (generalized) Causality 
and Time Slice axioms in Sec. [T] I am not aware of any methods that can be 
used to establish these results directly in the off-shell formalism (particularly for 
quasilinear systems) without first establishing its equivalence with the on-shell 
construction. Since understanding causality in field theories with quasilinear 
equations of motion is the main motivation for this work, we feel our choice to 
work in the on-shell formalism is justified. Of course, it is fully expected that 
future work will realize an exact duality between the two approaches. 

First suppose that that the gauge transformations are trivial. Consider a 
globally hyperbolic chronal cone bundle C — ?> M and a local functional A[(j)] 
(Sect. I5.2.7P on the slow patch T{F, C) of the space of field configurations. Fix 
some (f) S T{F, C) and consider the formal exterior differential and denote 

its value in T^T by a*. Then 5A[(/)](?/;) = (-0, a*) for any tangent vector € T^T. 

Definition 39. We refer to the support of the dual density a* as the local 
spacetime support of A[(j)\ at (j), supp^^ ^ A = suppa*. The global spacetime 



^ This discussion should be considered formal and analogous to the situation where all 
manifolds are finite dimensional. The justification of its conclusions is a matter of ongoing 
work in this field. 
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support of A[(f)] is the closure of the union of all local spacetime supports, 
suppm a = U,^er(F,c) supPm,,^ A. 

Wc may sometimes speak simply of spacetime support of A when the precise 
version of the notion is clear from context. 

More than functionals on field configurations, we are interested in functionals 
on the slow patch S{F,C) of the space of solutions Sh{F). So a functional 
on <S(F, C) is called local if it is the pullback of a local functional on field 
configurations along the inclusion S{F, C) C T{F, C). The definition of spacetime 
support for a functional on solutions is complicated by the fact that many local 
functionals on field configurations may pullback to the same functional A[(j)\ 
on solutions. This difficulty is resolved by deciding that A[(j)\ admits multiple 
spacetime supports. 

Definition 40. Consider a local functional A[(j)] on S{F,C) and any local func- 
tional A'[(f)] on r{F,C), which pulls back to A[(f>] on solutions. Then we say that 

suppj(^ A' is a local spacetime support of A at (j) and that supp^ A' is a global 
spacetime support of A. The collections of all local and global spacetime supports 
of A are denoted supp^ ^ A = [suppj^^ ^ A'] and suppj^^^ = [supp^i^ A'] , where 
the square brackets enclose a representative. 

Essentially, this definition states that suppj,^ ^ is the collection of supports of the 
representatives a* of the value of the formal exterior differential 6A[(p] = [a*] 
in TJS. When referring to the properties of the spacetime support of A[(/>], we 
are free to refer to any representative. For instance we say that a local func- 
tional A[(f)] on S{F, C) has a compact local spacetime support at (j) if supp^ ^ A 
contains a compact representative. The same goes for global spacetime support. 
A subtly, but importantly different notion is of globally compact local spacetime 
support, which means that the local spacetime support suppjy^ ^A has a com- 
pact representative for every (p G S{F,C). Also two such local functionals A[<j)\ 
and B[(j)\ have spacelike separated local spacetime supports at (j) if suppjy^ ,^^ 
and suppf^^B have representatives that are <;^spacelike separated. If A[(l)] and 
B[(j)] have spacelike separated local spacetime supports at every (j) € S{F,C), 
we say that A[(f)] and B[(j)] have globally spacelike separated local spacetime sup- 
ports. The same terminology also works for functionals on field configurations. 
On the other hand, if we say that A[<l)] and B[<j)] have C -spacelike separated 
global spacetime supports if supp^j^ A and suppj,^ B have representatives that 
arc C-spacclike separatcid. Of course, C-spacelike separation of global spacetime 
supports implies global spacelike separation of local spacetime supports, but 
the converse is does not hold. In fact, one can imagine local functionals A[<l)] 
and B[(j)] with globally spacelike separated local spacetime supports but with 
suppjv^ A = suppjv^ B = [M], because of the way suppj^^ ^ A and supp^j ^ B vary 
as a function of (j). 

In the presence of non-trivial gauge transformations, the gauge invariant 
cotangent spaces, T^S, are distinct from the ordinary ones, T^S. 

Definition 41. A local functional A'[(p] on T[F,C) is called gauge invariant if 
its formal exterior derivative SA'[(I)] takes values in the gauge invariant cotan- 
gent space T^T C T^F. Similarly, a local functional A[(i)\ on S{F,C) is gauge 
invariant if its formal exterior derivative dA[(p] takes values in T1S C TfS. 
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In particular, a local functional A[(j)\ (whether on solutions or field configura- 
tions) satisfies (5v4[(/)](g[e]) for any linearized gauge transformation at (j) with 
gauge parameter section e G ro(P). The notions of compact spacetime sup- 
port and C-spacelike separated supports specialize straightforwardly to gauge 
invariant local functionals. 

We are now ready to define the Poisson bracket of two gauge invariant local 
functionals A[(j)\ and B[(j)\, which we write simply as 

{A,B} = n{SA,6B). (242) 

In more detail, let </> e S{F, C) and denote the values of the formal exterior 
differentials at </> by 5A[(i)] = [a*] and = [^*], both in r|5. We then have 

{A,Bm = ^([«*], W\) = (E[ci1,/3*), (243) 

where we have invoked Lem. [5l]and E is the causal Green function of the Jacobi 
system at the dynamical linearization point 0. Incidentally, the last formula in 
Eq. (j243l) allows the extension of the Poisson bracket off-shell to gauge invariant 
functionals on field configurations jMllS^lHOll^ , though we will not discuss this 
possibility at the moment. 

A problem with this definition, and the reason we announced at the beginning 
of this section that the Poisson bracket would be defined only formally, is that 
{A, is in general not a local functional. That is because the Green function 

E is a non-local integral operator. Consider two local coordinate charts (x^ ^u"-) 
and {y^,u^) on F so that the dual densities a* and /3* have the coordinate 
expressions a*(a:)di and /3^(?/)dy, while E has the integral kernel E'^''{x,y). 
Then 

{A,Bm= [ E''''{x,y)a:{x)f3;{y)dxAdy. (244) 

As E''''{x,y) can be seen as a vector valued distribution on M x M, the above 
integral is well defined, but cannot be re-expressed as a local functional unless 
F,'^^{x,y) is supported only on the total diagonal x — y, which it obviously is 
not. It should be also clear that the product of two local functionals A[(/)]_B[(/)] is 
also not a local functional. 

To address both of the above problems, we expand the set of functionals 
under consideration to polynomials in local functionals. 

Definition 42. Let hoc{F, C) denote the space of local Junctionals on T{F, C) 
with compact local spacetime support and PolyLoc(i^, C) the polynomial alge- 
bra generated by LiOc{F,C). Similarly, Loc(i^) and Loc{F) refer to functionals 
on TniF). The elements o/ PolyLoc(F, C) and PolyLoc(i^) are referred to as 
polylocal functionals. 

As before for local functionals, polylocal functionals may be defined on field con- 
figurations as well as solutions. Recall that we did not give a precise definition of 
infinite dimensional manifolds in our category SOtan nor of the algebra of smooth 
functions on them. However, whatever precise definitions that are ultimately 
chosen, we require the following 
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Hypothesis 2. The functor of smooth functions C°° : VJlan -> £2l[g is such 
that C°^(T(F,C)) is a topological algebra containing PolyLoc(i^, C) as a dense 
sub-algebra, which also contains the Poisson brackets {A, B} of local functionals 
and such that {— ,— } extends continuously to all of C°°{T{F,C)). Analogous 
statements hold for C°°{Th{F)) and PolyLoc(-F') and upon replacement of field 
configurations with solutions. 

Remark 3. It is straight forward to generalize the formal exterior derivative from 
local functionals to polylocal ones, as well as to more general C°° functionals. 
The same goes for the notions of local spacetime support and global spacetime 
support. We silently make use of this generalization from now on. 

So, we will content ourselves with approximating the Poisson bracket {A, i?} 
by a sequence (or net or filter) of polylocal functionals Ci\4>\, with i from some 
index set /. This concludes the formal definition of the algebra of observables 
and the Poisson bracket on it. 

Definition 43. We call the algebra of gauge invariant smooth functions on the 
slow patch S{F, C) together with the Poisson bracket (|243p the Poisson algebra 
of observables of the classical field theory and denote it 

J-(F,C)-(C-(5(^^,C)),{}), (245) 

where the subscript est means that we only take functionals with compact local 
spacetime support. (Note that we have identified the gauge invariant functions 
on S{F,C) with functions on S{F,C), the space of gauge equivalence classes of 
solutions.) Denote the underlying commutative algebra of T{F,C) by A{F,C) = 
C^^{S{F,C)), which will also be referred to as the algebra of observables. 

Of course, any results concerning J-{F, C) and A{F, C) that we can establish in 
this paper will only be at the level of PolyLoc(F, C), which will be strengthened 
byHyp.m 

Finally, we are ready to state and prove the main result of this section. In 
quantum field theory, it usually goes under the name of microcausality: observ- 
ables that depend on fields in spacelike separated regions have vanishing com- 
mutators. The classical version replaces the commutator by the Poisson bracket. 
Note, however, that the result below is an important generalization of the known 
microcausality property of field theories with semilinear equations of motion, 
i.e., those with a non-dynamical causal structure. Usually, only the global space- 
time support is considered and, even if the local spacetime support were to be 
considered, the notion of space-like separation is independent of the dynamical 
fields and so does not change when the observables are linearized about different 
solutions (for instance, in perturbation theory). For field theories with quasi- 
linear equations of motion, and hence dynamical causal structures, the notion 
of spacelike separation is intrinsically changes from solution to solution on the 
phase space. Thus, to even have a suitable generalization of the spacelike sep- 
aration hypothesis for microcausality, we are forced to consider local spacetime 
support. 

Theorem 7 (Classical microcausaUty) . Consider two observables A and B 
belonging to the Poisson algebra of observables J-{F,C). If their local spacetime 
supports are spacelike separated, then their Poisson bracket vanishes identically. 



{A,B)=Q. 



(246) 
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Proof. Consider any solution (f> e S{F, C). Let the dual densities a* , (3* e ro(_F*) 
represent the formal exterior derivatives of the observables, SA[(j>] = [a*] and 
5i3[(/>] = [/3*]. Pick a pair of local coordinates and {y^,u'') on F such that 

the dual densities have the components a*{x)dx and l3^{y)dy. By hypothesis, 
supp a* and supp are 0-spacelike separated. 

From the Peierls formula for the Poisson bracket (I243P 

{A,Bm^ f E''''{x,y)a:{x)l3;iy)dxAdy, (247) 

J MxM 

where E°''(a;, y) is the integral kernel of the causal Green function E of the Jacobi 
system linearized at (j>. However, by Eq. (|164p . E = G o f , where G is the causal 
Green function of a linear symmetric hyperbolic system with adapted equation 
form (/, F*), whose principal symbol precisely determines the causal structure of 
the solution cj) (Sect. [44]). Recall that the causal Green function is the difference 
of the retarded and advanced Green functions, G — G+ — G_ . Thus, the support 
of G, and hence E and Fi°'^{x,y), since f is a differential operator and does 
not increase supports, excludes any pair of points (x, y) that are (/)-spacelike 
separated. 

It is now an easy conclusion that, according to the formula (|247l) . the Poisson 
bracket {^4, vanishes at € S{F, C). Since the solution was arbitrary, it 

follows that the Poisson bracket vanishes identically, {A, B} = 0. □ 

A coarser version of the above result, though stated in the more familiar 
terms of global spacetime support follows immediately once the definitions are 
unwound. 

Corollary 4. Consider two observables A and B belonging to the Poisson al- 
gebra of observables T{F,C). If their global spacetime supports are C -spacelike 
separated, then their Poisson bracket vanishes identically, {A,B} = 0. 

Remark ^. The notion of spacetime support in previous work on the functional 
approach to classical and quantum field theory [23ll47ll95] corresponds to what 
we call global spacetime support. The notion of local spacetime support and of 
spacelike separation of local spacetime supports seems to have been overlooked 
until now, with the notable exception of ^illW] . The distinction between local 
and global spacetime support is particularly important for gauge theories and 
theories with dynamical causal structures, both exemplified by GR. It is well 
known that there are no gauge invariant observables in GR with compact global 
spacetime support, a significant technical complication. However, the door is 
still open to discover a large, technically convenient class of observables with 
compact local spacetime support. 

Another important but simple result is 

Lemma 22 (Covariance). Consider two globally hyperbolic chronal cone bun- 
dles C — > M' and C — > M , where M' C M and C is faster than C. If we 
restrict our Lagrangian and PDF system to M' , there are induced maps of so- 
lution spaces S{F,C) S{F',C') and algebras A{F\C") AiF,C), where 
F' = F\m' is the field bundle over M' . These maps induce morphisms of Pois- 
son manifolds {S{F,C), 11) — >■ {S{F,C'), 11') and a Poisson homomorphism 
{A{F' ,C'),{Y) = F{F\C') T{F,C), where 77' is also defined by the Peierls 
formula and in turn defines {—,—}'. 
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Proof. The agreement of the Poisson bivectors 77' and 77 follows directly from 
the Peierls formula and the covariance lemma for the causal Green function, 
Lem [SI This agreement immediately implies that the induced algebra homo- 
morphism A{F' , C) — !• ^(7^, C) is automatically a Poisson homomorphism with 
respect to the brackets {—,—}' and {—,—}. □ 

5.3.6. The global phase space. Here we take the final step in the construction 
of the phase space of the classical field theory and the Poisson algebra of ob- 
servables on it. We appeal to the results obtained for each slow patch S{F, C) 
and recall that these patches form an open cover of the total space of solutions 
Sh{F), according to Hyp. [1] Moreover, each of the open patches carries a for- 
mally smooth presymplectic form r2c[0] (which becomes a symplectic one Dc[4>\) 
when projected to S{F,C), Sects. [573.31 and [5.3.41 The symplectic structure is 
translated to a Poisson structure 77c [0] by pointwise-inversion (smoothness and 
Jacobi identity are presumed to follow from the corresponding properties of the 
symplectic form). Sect. 15. 3T^ It remains only to check that the symplectic and 
Poisson structures agree as these patches are glued together. In a later section. 
Sect. 16.31 this construction via slow patches will be interpreted in categorical 
terms as a colimit, which will play an important role in the generalized Causal- 
ity property of classical field theory with quasilinear equations of motion. 

Consider two chronally comparable globally hyperbolic chronal cone bundles 
Ci — >■ M and C2 — >■ M, as well as their intersection Ci n C2 = C3 — >■ M, also 
a globally hyperbolic cone bundle. The each we can associate the corresponding 
patch of Ci-slow solutions 5(7^, Ci) modulo gauge transformations. On each such 
patch, we have the symplectic and Poisson tensors f^Ci [<P] and Ud [<j>] ■ The slow 
patches intersect as S{F,Ci) n 5(F, C2) = S{F,C3). Therefore, we must check 
whether >!?Ci[0] slW agree when restricted to iS(F, C3). Recall that each of the 
symplectic forms was defined with the help of a Ci-Cauchy surface Si C M. 
Each of the three symplectic forms is defined by exactly the same formula. 
Cor. [21 with the exception that the integration surface S is replaced by the 
corresponding Si. By construction. Si and S2 are both Ca-Cauchy and so lie 
in the same Ca-spacelike homology class. Since Dc^ [4>\ depends only on that 
homology class, we can directly establish that 

f2c^[^] = f2cM ^ ^cM (248) 

on iS(7^, C3). The Poisson tensors Ud [(j)] necessarily satisfy the same property. 

From the above discussion, we can conclude that the global space of glob- 
ally hyperbolic solutions SniF) is endowed with both symplectic and Poisson 
structures, and 77[(/)], defined by their agreement with the corresponding 
tensors on the slow patches S{F,C) C Sh{F). We call the symplectic ma- 
nifold Vh[F) = {Sh{F), fi) the global phase space of the classical field the- 
ory under consideration. The Poisson algebra of smooth functions on it, with 
Poisson structure defined by 77 [</)], is called the global algebra of observables 
Fh{F) = (ylff (7^), {}), where Ah{F) ^ C^^{SH{F)) is the underlying commu- 
tative algebra. 

The notions of local spacetime supports of observables and their global space- 
like separation, as in Defs. [5niand[TO translate directly from the slow patches to 
the global phase space. Therefore, as a direct consequence of Thm. [T] we have 
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Corollary 5 (Global classical microcausality). Consider two observables A 
and B belonging to the global Poisson algebra of observables Th{F). If their local 
spacetime supports are globally spacelike separated, then their Poisson bracket 
vanishes identically, 

{A,B}= 0. (249) 

6. Natural Variational PDE Systems and Functoriality 

In this section we take account of the various categories and functors introduced 
so far and introduce a few more. The main goal is to make explicit the functorial 
nature of all the steps in the construction of a classical field theory. For reference, 
some basic information about categories and functors can be found in |1091I21] . 

6.1. Natural bundles. So far, we have looked at a fixed spacetime manifold M , 
a fixed field vector bundle over F — >■ M, and a fixed PDE system £ — >■ M on 
it, without concerning ourselves how the field bundle of the PDE system should 
change when the spacetime manifold is changed. To make a connection with 
covariant field theory, which functorially assigns algebras to spacetime manifolds, 
it is helpful to look at systems of partial differential equations that can also be 
assigned to spacetime manifolds functorially. As before, we presume that all fields 
are sections of vector bundles. The discussion below can be easily adapted to 
arbitrary smooth bundles. There are two issues that need to be addressed here: 
vector bundles and differential operators on them. Given an open embedding 
of manifolds M — > M', there is a priori no natural map from vector bundles 
F M and F' M, even if the fibers of F and F' are modeled on the same 
vector space. The absence of such maps is an obstruction to relating sections of 
F —i' M to sections of F' — M. The needed bundle maps must to be specified 
along with the typical fiber and the total space topology, which leads to the 
notion of a natural bundle. It is common practice to use a connection to define 
differential operators acting on sections of a vector bundle (that practice was 
adopted in |50| . instance). But, once again, for an arbitrary vector bundle, there 
is no natural choice of connection. It then becomes convenient to work with jet 
bundles instead, as we have elected to do here. The switch between these points 
of view is straightforward [751 §17]. 

A field configuration corresponds to a section of a vector bundle F M over 
a spacetime manifold M . To implement covariance with respect to arbitrary base 
space morphisms x - M M' we need to be able to turn sections of F — )■ M into 
sections of the corresponding bundle F' — M', restricted to the image x(M). 
This requires a uniquely determined corresponding bundle morphism, F ~> F' 
that fits in the commutative diagram 

F ¥F' 

(250) 

M — ^M' . 

In essence, we cannot allow F and F' to be arbitrary vector bundles; we need 
them to be associated functorially to AI and M' respectively, which motivates 
the following definition. 
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Definition 44. A natural vector bundle F over n- dimensional manifolds is 
a functor from the category of manifolds to the category of vector bundles, 
F: 9Jtan" — > 2J*Bn£)l, that is right-inverse to the forgetful base space functor, 
U: {A^ M) ^ M, that is, H o = id. 

Examples of natural vector bundles are the tangent bundle T, the cotangent 
T* . Operations on natural vector bundles like direct sums, direct products, linear 
duals, tensor products and jet extensions also define natural bundles. These and 
other examples are discussed extensively in |73) . Given a natural vector bundle 
F: man 2J<Bn0l, its fc-jet bundle M ^ J^F{M) is also natural. 

As we have seen in Sect. 13.11 to specify a fc-th order PDE system on a bundle 
F, we need only specify a bundle map from J^F to an equation bundle E. 
If F and E are natural bundles, then a natural transformation between J^F 
and E defines precisely these kinds of morphisms, J^F{M) — >■ E{M) in a way 
covariant with diffeomorphisms of M. Recall that it is sometimes convenient to 

see M ^ M, a trivial bundle with zero dimensional fiber, as an object in !Bnt)( 

and id: M i-)- {M ^ M) as the natural identity bundle. 

Definition 45. A natural bundle map / between natural (^ndl-valued) bundles 
F and E is a natural transformation f:F^E. 

A natural section f of a natural bundle F, is a natural bundle map f : M F , 
where the domain is treated as the identity bundle. 

Two examples of natural sections, given say a natural vector bundle F , are the 
zero section and the identity section of F* ® F. An example of a natural bundle 
map, is the dual pairing (—,—): F* ® F — > R, where M is the trivial real line 
bundle. Other examples are provided by smooth real functions / : M — >■ R, which 
can be promoted to a natural bundle map by applying it fiberwise to the natural 
trivial bundle M x R. 

We now have the necessary concepts to introduce natural variational hyper- 
bolic PDE systems. 

Definition 46. Let F, B , V = F x B and E be respectively natural dynamical 
field, background field, total field and constraint vector bundles. 

A natural first order, quasilincar, symmetric hyperbolic PDE system with 
constraints f ® c is a natural bundle map f (B c: J^F x J°°B F* (B E, such 
that the bundle maps fM, for any M in 9Jlan", satisfy Def. \l(A for some natural, 
non-empty subset of each fiber of J°° B{M). 

A natural Lagrangian density is a natural bundle map C: J^F x J°°B — > 

yl"M. 

A natural symmetric hyperbolic variational PDE system is a pair (/ ® c, £) 
where C is a natural Lagrangian density and f is a natural first order, quasilin- 
car, symmetric hyperbolic PDE system that is naturally equivalent (Def. to 
gauge fixed Euler- Lagrange equations of C (Sect. [5.3.2]) . 

The introduction of background fields is sometimes necessary to make the 
Lagrangian natural. One may also need extra background fields to put the corre- 
sponding Euler-Lagrange equations into symmetric hyperbolic form in a natural 
way, as discussed for example in l&O] . The example of a scalar wave equation on 
a curved background is discussed in some detail at the end of the next section. 

A field section ij; on M will from now on always be a section of the total 
field bundle ip — {ip, /3) : M^V = VxB. As a. consequence of the dependence 
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of the principal symbol / of / on the value of background fields, implies that 
the naturally defined chronal and spacelike cone bundles r{M) and r®{M) are 
cone bundles over F x J°°B{M). 

An unfortunate complication in the presence of background fields is that not 
all background field configurations are compatible either with the intrinsic inte- 
grability conditions of the PDE system (like the restriction to Einstein vacuum 
backgrounds in the Rarita-Schwinger system [5 7) ) or with symmetric hyper- 
bolicity (like the restriction to Lorentzian signature of the background metric 
for wave equations). With this in mind, we say that a background field sec- 
tion f3: M B{M) is admissible when (a) there exists a total field section 
i} — {(j), /?) : A/ — > that is globally hyperbolic and solves the PDE system 
f — 0, (b) for any ^/i-Cauchy surface S C M, there exists an open neighborhood 
of D 17 such that the initial data ^lu of the space of globally hyperbolic solu- 
tions of the form (^, /3) on N forms an open neighborhood of (pls in the space of 
initial data allowed by the constraints c = (without taking further integrability 
conditions into account). 



6.2. Functoriality. In this section we (a) systematically summarize some of the 
concepts and notations used throughout this paper, (b) systematically summa- 
rized various constructions that appear as steps in the construction of a classical 
field theory, and (c) remark on the categorical and functorial properties of these 
objects and constructions. This section will serve as a reference for the later 
discussion of how LCFT axioms, including causality, translate to field theories 
with dynamical causal structure. 

1. DJlan, SOTan", D}lan°°: Category of smooth manifolds, subcategories of n-di- 
mensional and infinite dimensional manifolds. Subscript e stands for restric- 
tion of morphisms to open embeddings. 

2. Q3nc)[: Category of smooth bundles, a subcategory of DJlan. Fibered over DJlan 
with the base space functoiF"! 11 : (i7 — > M) ^ M. 

3. 2J*BnOl: Category of vector bundles, a subcategory of SnOl. 

4. £*BnD[, £f)t*Bn£)l, 6pQ3nO[: Subcategories of *8nD[e of cone bundles, chronal 
cone bundles, spacelike cone bundles. Subscripts sc and h indicate stable 
chronality and global hyperbolicity and superscript " denotes the restriction 
to n-dimensional base manifolds. Morphisms in £f)r*8nc)[ and 6pQ3ni)[ must 
be chronally convex. As functors from 931an, the images of M i— )■ £()r*Bnc)[(M) 
and M ©p*8 nO[(M) indicate the respective subcategories of cone bundles 
over M. 

5. J'^ : 2J*8ni)[ — ^ 2J*8n5[: Functor of jet prolongation of a vector bundle. 

6. r, Fo : 2J*Bn0l — > S[Han°°: Contravariant functors of sections, sections with 
compact support. 

7. rsc,+ - ■ 23*8n5[ Xn Gp^ndl^^ 9Jlan°°: Contravariant functors of sections 
with spacelike compact, retarded (-I-), or advanced (— ) supports. The objects 
of the domain category can be identified with pairs {F, C® ) , with F M 
a vector bundle and C® — )■ M is a spacelike cone bundle over the same M. 
Throughout the paper we have also used the notation Tsc+.-iFiC), where 
C — ?► Af is a chronal cone bundle. These notations can be used interchangeably 



We will re-usc the symbol EI to generically denote a forgetful functor. 
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because of the duality between chronal and spacelike cone bundles. However, 
these sections functors are covariant only if applied to spacelike cone bundles. 

8. B,F,V,E: 9Jlon" 2J*Bnc)l: Natural bundle functors of background, dynam- 
ical, and total {V = B ® F) field bundles as well as the equation bundle. 

9. C: J°°F X J°°B — >• A"M: Natural bundle section, defining a covariant La- 
grangian density. It can also be interpreted as an differential form £ € 

® B) in the variational bicomplex. 

10. /®c: J^Fx J°°B — >• Natural bundle section defining a covariant first 
order, quasilinear, symmetric hyperbolic PDE system, which is equivalent to 
gauge fixed Euler-Lagrange equations of C. 

11. *8tgr: Category of manifolds augmented by collections of admissible back- 
ground field configurations. Objects are pairs (M, B) with M an n-manifold 
and B C T{B{M)). Each element (3 G B must satisfy all the necessary integra- 
bility conditions and be compatible with the existence of globally hyperbolic 
dynamical field solutions. A morphism x- {M,B) — ?> {M',B') determines an 
open embedding x - M ^ M' compatible with background fields, x*{^') ^ ^ 
being a closed embedding (to play nice with C°° , as discussed in item [T^ 
below). *8tgr is fibered over 97lon" with respect to the forgetful functor to 
H: (M,e) ^ M. 

12. 6p*B{gr — (6p*8n5[^ Xn *8lgr)c: Category of admissible background field 
configurations, equipped with a globally hyperbolic spacelike cone bundle. As 
functor from ^Bfigr, the image of M i— 6p*Blgr(Af, B) indicates the subcate- 
gory of different ways of equipping (Af, B) with a spacelike cone bundle. The 
objects of 6p*B{gr will be denoted a.s M = (M, C®,B). The forgetful functor 
is again denoted H: 6p<8«gr ^ <B«gr, with U{M) = U{M,C®,B) = {M,B). 
The subscript c indicates that morphisms are restricted to those that are 
chronally compatible fPef. 

13. ©r^mp: The generalized category of symplectic manifolds, a subcategory of 
DJlan. It includes symplectic manifolds as well as manifolds foliated by sym- 
plectic leaves (which is actually the same as the category of regular Poisson 
manifolds). The morphisms are leaf preserving symplectomorphisms. The for- 
getful functor to manifolds is 11: ©tjmp — DJlan. 

14. £2tlg, *Poiss: Categories of commutative algebras and Poisson algebras, with 
the forgetful functor 11: Cpoiss £2t[g. Though, keeping away from functional 
analytic details, we implicitly treat these categories as though the objects 
are equipped with topology and the homomorphisms are continuous. Such 
continuity requirements also identify a tensor product (identified with the 
categorical coproduct), which is compatible with the categorical product in 
Dylan, as specified in Hyp. [3] below. The extension of the tensor product to 
*Poiss is the independent subsystems tensor product. 

15. C°°, C^(, POISS: Functors of smooth functions of ordinary manifolds, that 
restricted to compact spacetime support, and that equipped with Poisson 
bracket on symplectically foliated manifolds, which fit in the following com- 
mutative diagram: 



6t}mp 



POISS 



n 



n 



(251) 



c: 



est 



man 
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The definition of C^^{M) is unambiguous wlien M is a finite dimensional 
manifold. On the other hand, its definition requires some functional analytical 
detail not tackled here. However, we have already postulated that, whatever 
ultimate definition is adopted for it, it will have to satisfy Hyp. [21 In addi- 
tion to that, let us make explicit additional hypotheses, of a more category 
theoretic nature, that we presume it would satisfy. 

Hypothesis 3. The functor gives a sheaf of commutative algebras Jlllf 
on a manifold M when applied to the open subsets of M . Moreover, maps 
surjections of manifolds to injective homomorphisms ("left exactness") and 
maps closed embeddings of manifolds to surjective homomorphisms ("right 
exactness" or smooth Urysohn lemma). Products of manifolds are taken to 
tensor products of algebras, C^t{M x N) ^ C^t{M) (S) C^^tW ■ "Transverse" 
pullbacks of manifolds are taken to pushouts of algebras. 

Remark 5. The statement about transverse pullbacks is standard in differ- 
ential geometry of finite dimensional manifolds [Sni Thm.2.8]. However, the 
notion of transversality is much more subtle for infinite dimensional mani- 
folds [501I771I55] . In fact, the notions of transversality discussed in these ref- 
erences may not be general enough to encompass the situation where it is 
needed in the formulation of the generalized Time Slice property in Sect. 17.21 
At this point, we have no choice but to leave the existence of the appropriate 
notion of transversality as a conjecture. 

16. S : 6p*8tgr OTan, V : &p^iQV^ ©tjmp: The contravariant functor assign- 
ing S{A4), the set of all C*-slow solutions of the gauge fixed Euler-Lagrange 
equations (equivalently of the corresponding constrained symmetric hyper- 
bolic system) compatible with the background field configurations in B, to 
the object M = {M,C*,B) of 6p®l0c. The contravariant functor V aug- 
ments the images of S with symplectic structure defined by the Lagrangian 
C. A fixed element /3 G B singles out a symplectic leaf in S{Ai). It is the 
dependence of solutions on background fields that forces us to consider gen- 
eralized symplectic manifolds. Note, however, that a symplectically foliated 
manifold must have at least dimension two. On the other hand, many of the 
S{Ai) spaces will be empty sets, hence cannot be the underlying manifolds 
of images of V. For that reason, the domain category of V is restricted to 
©pCBSgr^, where the subscript c indicates the largest subcategory such that 
no object in the image 5(6p!B{grc) is empty. 

17. y{M,B): &pmgv{M,B) Man, 3^{M,B): &\)^igx^{M,B) ©omp: The 
contravariant functor whose image in 9Han, resp. (5r)mp, consists of the di- 
agram S{M) —T' S{M'), resp. V{M) -> V{M'), of ah slow patches on the 
total solution space Sh{M, B), resp. phase space Vh{M, B), with the natural 
inclusions between them. 

18. Sh ■ *B{gr — > 9Jtan, Vh ■ ^SJgr^ — > ©tjmp: The contravariant functor of all 
globally hyperbolic solutions that are compatible with specified background 
field configurations. The contravariant functor Vh augments the images of 
Sh with symplectic structure, as above, where we have defined *B£gr^ = 
H(©p*B?grj,). Each S{A4) is an open slow patch of the total solution space 
Sh{M,B). The natural isomorphisms Sh — limoj5^ and Vh — fimo^ are 
verified in Sect. 16.31 These functors fit into the following commutative dia- 
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gram: 



> Man . (252) 



19. T = POISSo-p, A = C^t"^ = no J": Covariant functors of the Poisson alge- 
bra of observables on S{M) and its underlying commutative algebra, specified 
by a slow patch M = (M, C®,B). 

20. ^(M,i3): 6p<B«0r(M,e) ^ CSllg, ^(M,6): 6p<8«flr^(M, 6) ^ <Poiss: The 
covariant functor whose image in CSllg, resp. ^Poiss, consists of the diagram 
A{M) A{M'), resp. J^iM) F{M'), of all algebras of observables asso- 
ciated to slow patches in the total algebra Ah{M,B), resp. Poisson algebra 
Fh{M,B), with the natural projections between them. 

21. Fh ^ POISS oPh, Ah = C^t o Sh = Ho Th- The covariant functors of 
the Poisson algebra of observables and its underlying commutative algebra 
of the total phase space Sh{M,B). The natural isomorphisms — ^m oJF 
and Ah = lim o.s/ are verified in Sect. [^751 These functors fit in the following 
commutative diagram: 



> eraifl . (253) 



22. Stgr^, *B?0r*, *84gr': Subcategories of *8tgr. As defined in a preceding item, 
Q3tgCc — II{&p^tgtJ. We define S^gr* as the subcategory where the collec- 
tions of admissible background configurations are singletons, B = {/?}, and 
QStgr' is the image of the category QStgr* , where the auxiliary components of 
the background fields have been thrown away. 

With all the categorical notions summarized above, it should now be clear 
how to recover the standard formulation of the LCFT axioms given in Sect. [21 
Consider a semilinear wave equation on a curved background with Lagrangian 
density 

C=-^{g''''d^^d,cj, + V{cf>))V^dx. (254) 

The trivial dynamical field bundle F: A/ i-> R x M is clearly natural. However, 
the Lagrangian density C only becomes natural if we include the metric tensor 
among the background fields, that is we set B: M ^ S'^T*M, which is also 




«B«0r, > Q5«0r 



epStgr^ -^-^ 6p«8«0r 




n *poiss 




Ah 



^tgv, — > <Bt0c 
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clearly natural. The Euler-Lagrange equations automatically give a natural hy- 
perbolic PDE system. (This system is normally hyperbolic [71124) . For simplicity 
we do not put this system into symmetric hyperbolic form, which would simply 
require extending the dynamical and background field bundles [50].) 

It is simple to check that restricting to the subcategory Stgt* of Stflc^, 
where objects {M, B) are equipped with a single background field configuration 
B = {g}, forces the allowed morphisms to be only causal isometrics (cf. Def. [2l 
Dcf.[Ml Lem.[5]and items [TTl [T^ above) . This means that we have an equivalence 
of categories ©tgr* = &lobS)t)p^. By construction then the covariant functor 

Th ■■ *8«0r* ^ ©[obijtjp^ «Poiss (255) 

has the potential to be a classical LCFT. It remains only to check the axioms of 
Def. [SJ These checks are carried out at the end of Sect. 17.41 



6.3. Limits, colimits and the global phase space. We have seen earlier, more 
specifically in Sect. 15. that the slow patches S{A4) constitute an open cover 
of the global solution space Sh{M,B). How can this relation be restated for 
the algebras Ah{M,B) and A{MtB), or the Poisson algebras J^h{M,B) and 
J-{A4)7 The answer is not immediately obvious. However, the situation becomes 
more clear when one realizes that the desired constructions can be formulated 
in terms of the categorical notions of limit and colimit. Categorical limits (also 
projective or inverse limits) and colimits (also inductive or direct limits) are 
briefly introduced in Sect. [BJ 

A great advantage of considering the Poisson algebras J-iM) is the simpli- 
fied microcausality property. Cor. 21 while only the more refined microcausality 
property of Cor. IS] survives for Fh{M,B). On the other hand, the J-{M.) for 
different cone bundles C ^ M are included as subalgebras of Fh{M, B), though 
not as independent ones. The situation is strongly parallel to the fact that the 
slow patches S{M) are open subsets of Sh{M, B), though with non-trivial over- 
laps. The parallel is made precise by recognizing that the globally hyperbolic 
slow patches S{A4), as well as their symplectic analogs 'P{Ai), and the inclu- 
sions between them constitute the diagram S^{M, B) in CDTan, and respectively 
,^(M, B). Similarly, the algebras A{M) and y-{M) and the projections between 
them constitute diagrams (M, B) and ^(Af, B) in €S\Xq and *Poiss, respec- 
tively. The open cover of Sh {M, B) by slow patches and the agreement of sym- 
plectic structure on their overlaps. Sect. I5.3.6[ ensures the following limit and 
colimit identities hold 

SHiM,B) '^\\^.y{M,B) and VHiM,B) ^lw^^{M,B). (256) 

Similarly, in the algebraic categories, the sheaf property of ensures that the 
following colimit identities hold 

AHiM,B)=]^£/{M,B) and THiM,B) = ]^^{M,B). (257) 

The same category ©p*84gv(A'f, B) contravariantly indexes the geometric colimit 
and covariantly indexes the algebraic limit in the left column, while the same 
category &p^tQV^{M, B) contravariantly indexes the geometric colimit and co- 
variantly indexes the algebraic limit in the right column. Recall that objects 
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M — {M,C®,B) of these indexing categories specify a spacelike cone bundle 
C® — i' M. If a chronal cone bundle were to be specified instead, the above lim- 
its and colimits would have been swapped, since 6pQ3 nO[(M) is the opposite 
category of £f)r*BnO((M). However, it is more convenient to use spacelike cone 
bundles, in light of the considerations below. 

Remark 6. Colimits appear in the framework of LCFT already at the level of 
spacetimes. Namely, consider a globally hyperbolic Lorentzian manifold (M, {g}) 
in the category ©lobf)r)pg. Let ^ be the diagram whose image in 0[obior)pg con- 
sists of all globally hyperbolic Lorentzian submanifolds (Mi,{g}i) of (Af, {(/}), 
with inclusions as morphisms. The colimit lim^ is isomorphic to {M, {g}) it- 
self. In fact one recovers (M, {g}) as the colimit even if one restricts the objects 
of ^ to be simple in particular ways, like being topologically trivial, having 
bounded radius, being the interior of a causal diamond, etc. Each kind of simplic- 
ity translates to a correspondingly simple property of the algebra of observables 
T{Mi, {g}i) associated by a LCFT. Thus, just as the non-trivial topological or 
geometric structure of (M, {g}) is encoded in the diagram ^ in ©[obi^tjp^, the 
non-trivial algebraic structure of T{M, {g}) = lii^ o ^ is encoded in the di- 
agram o ^ in *PoisSj . (We use the colimit in both categories because is a 
covariant functor.) This approach underlies the investigation of superselection 
sectors of LCQFT ^56j Ch.IV], j26j and the quantization of vector field theories 
on manifolds of non-trivial topology [BH Apx.A], [36] . 

It is sometimes convenient to simultaneously display both the limit and col- 
imit identities: 



S„{M,B) = lim(limo^ o ^), Vh{M, B) = lim(limo^ o .^,), (258) 



where we define ^ {^c) is the diagram in QStgt (Stgc^) that consists of all sub- 
objects {M',B') — {M,B). Similar ideas about constructing the global phase 
space and algebra of observables, using an open cover of the spacetime and 
a cover of the space of solutions also recently appeared in [24]. Though, in 
distinction, the covers on SniM, B) are considered there more general than ours. 
The cover elements are not required to be open (only Gs) and could be much 
more refined than our slow patches. 

7. Causality Axiom in Locally Covariant Classical Field Theory 

The axioms in the standard definition of LCFT, Def. [5j make reference to a 
fixed causal structure induced by the future timelike cones of a background 
Lorentzian metric. It is clear that they cannot be translated to direct conditions 
on the functor Fr (or Vh), which can be constructed for theories without a 
fixed background causal structure. On the other hand, these axioms can be 
straightforwardly translated to conditions on the functor T (or 7^), since objects 
in its domain category (5p*B{gr all have an externally fixed causal structure given 
by a spacelike cone bundle. Below, we prove theorems about the properties of 
the functor F, which lead naturally to generalizations of the LCFT axioms. The 
functor J-H = limo^ inherits these properties through the categorical limit 
construction. 




(259) 
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7.1. Isotony. Given a morphism x- ^ A4' in &p^tgx, is the corresponding 
morphism x* ■ S{A4) — S{M') surjective? The answer is not always positive; it 
depends on the properties of the field theory in question. In terms of solutions, 
X* is surjective precisely when every C®-slow solution on M is the pullback along 
X of a C"®-slow solution on M' (equivalently, can be extended with respect to 
X to a C"®-slow solution on M', where there is no presumption of uniqueness of 
the extension). Broadly speaking, the extension of a particular solution on M 
fails if the would be extension develops singularities in M'. 

There are different kinds of singularities, including topological, smooth, ge- 
ometric and analytical. (T) For instance, depending on the structure of the 
equations of motion, there may be a topological obstruction in extending a so- 
lution from M to M', if M is not contractible (S) Also, it is possible that a 
solution becomes unbounded or non-smooth near the boundary of x(-^) 
and hence cannot be extended continuously to M' . (A) The solution may even 
be completely regular on x(M) and its boundary, yet be forced by the dynamics 
of the equations of motion to blow up when extended to M'. The formation of 
shocks in fluid dynamics and of black hole singularities in GR are prime exam- 
ples of this phenomenon. (C) On the other hand, it is possible for a solution to 
extend from xi^) to M' , yet the fixed cone bundle C* is not fast enough to be 
compatible with the extension or simply that the image x{M) is not chronally 
convex in M' fPef. lM)) . 

It is of course desirable to formulate a set of necessary and sufficient conditions 
for extensibility to hold. Unfortunately, that is in general a very difficult problem. 
In particular, the study of blow ups of type (A) constitute a major active branch 
of modern PDE theory. Moreover, the existence singularities of type (T) and 
(A) is determined by the structure of the PDEs, which is fixed by the relevant 
physics. As such, we may just have to learn to live with them. Inextensibility 
of type (S) depends on whether we allow asymptotically irregular solutions in 
S{M). This choice is related to the issues of the choices of the structure of the 
domains M and topologies on the space of solutions, which were discussed in 
Sect. 14.2] in relation to the openness of slow patches. It is also possible that the 
failure of surjectivity of x* due to blow up of type (S) is precisely masked by 
the application of the functor that constructs the algebra of observables on 
the given slow patch of Sh (M, B) [TTII^ . As in this preceding discussion, we do 
not address it here directly. Finally, blow up of type (C) is entirely within our 
control. Since the problematic causal cone bundle C* plays only an auxiliary role 
in the ultimate construction of Sh{M', B'), it is no loss to introduce a faster cone 
bundle, say C[* that might be fast enough to be compatible with an extension 
of every solution in S{M). 

Unfortunately, as can be seen from the above discussion we cannot prove 
a general theorem about the surjectivity of the morphism x*- Therefore, we 
introduce this property as an additional hypothesis. 

Definition 47. A morphism x in ©p^Btgr, is said to he extensible if this mor- 
phism factors through 



M 




(260) 



Characteristics, Conal Geometry and Causality in Locally Covariant Field Theory 



83 



such that M' = {M',C'®,B'), M[ = {M',C[®,B'), and the morphism xl = 
S{xi) is surjective. 

Provided this hypothesis holds, there are no other obstacles in proving the fol- 
lowing generalization of the Isotony property for the functor J-. 

Theorem 8 (Generalized Isotony). Consider a morphism x in ©p^Btgr. If x 

is extensible, the morphism A{xi) is an injective homomorphism. Moreover, if 
S{M) is not empty, x ^.nd Xi olso morphisms in &p^tQV^, the morphisms 
S{x) and S{xi) lift to V{x) and V{xi) in &^mp, and A{x) and A{xi) lift to 
morphisms in *Poiss. 

Proof. Since x is extensible, the induced morphism S^M'i) — > S{Jv[) is sur- 
jective. By Hyp. [31 takes this surjection to the injective homomorphism 
A{M.) — A{M.'i). If the corresponding solution spaces are not empty, the co- 
variance lemma, Lem. [22l implies that the appropriate morphisms lift to 6t)mp 
and *Poiss. □. 



7.2. Time Slice. Given a morphism x'- ^ in ©p*B?0r, is the correspond- 
ing morphism S{x)'. S{M) — )■ S{M.') injective? If the answer were positive, it 
would mean that every C"®-slow solution on Af is uniquely determined by its 
restriction to x(M) C M' . Clearly, there is no hope that this is the case unless 
x(M) contains a C"®-Cauchy surface in M' . 

Definition 48. A morphism x in Sp*B{gr is called Cauchy surjective or a 
Cauchy surjection if for the corresponding morphism in (3p*8nc)[^, 

(J® y (ji® 

(261) 

M — ^M' , 

there exists a C® -Cauchy surface U d M whose image U' — %(-£') C M' is 
a C® -Cauchy surface. (Recall that all cone bundles in 6p*8nc)[^ are globally 
hyperbolic.) 

Unfortunately, injectivity of S{x) does not guarantee surjectivity of A(x)j 
which would imply that every smooth function on the image of S{A4') extends 
to a smooth function on S{A4). However, this property is known to fail unless 
S{x) is a closed embeddincF^. If S{x) is also surjective (x is extensible), this 
is automatic. However, if S{x) fails to be injective, then the complement of its 
image (the solutions that blow up when extended from xi^'I) to -^^') is rarely 
an open set. 

Fortunately, a slightly more complicated diagram in 6pQ3tgc does guarantee 
a closed embedding. 



^-^ A simple illustration is the inclusion (0, 1) C M and any function that is unbounded near 
an end of the open interval. 
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Definition 49. Consider morphisms xi 0,1^^ X2 in SpStgr that fit into a pull- 
back-pushout diagram in (5p*Btgr, 

M3 > M2 

X2 (262) 

That is, the images Xi{^h) form an open cover of M' and A/3 is isomorphic to 
their intersection. In terms of the spacelike cones, we have C® — Cf + C® and 
C® = cf n C® ( cf. Sect. \4-l% where for brevity we have omitted the pullbacks 
or pushforwards with respect to the appropriate morphisms, while S3 = Si U B2 
and B' ~ Bi D B2 on M3, with the same shorthand. We call the diagram (|262p 
a Cauchy pushout for M' if all morphisms are Cauchy surjections. 

Given a Cauchy pushout for M' , the maps S{M') — >■ S{Mi) are injective by 
Cauchy surjectivity. However, they are not guaranteed to be closed embeddings. 
On the other hand, consider the following diagram of spaces of field configu- 
rations in DJlan, where Vi = V{Mi) and V = V{M') are total field bundles 
naturally associated to their base manifolds. 



^{V3)< T{V2) 




(263) 



r(i/i) X r{V2) , 

which is most definitely a pullback diagram, with the canonical dotted line 
morphism a closed embedding. The set theoretic pull back property is assured 
because any pair of smooth sections on Mi and M2 that agree on M3 must 
glue togetheiQ to a smooth section on M3. Now, consider a point (V'i,V'2) in 
r(Vi) X r(V2) that is not in the image of T{V'), that is, it determines two sec- 
tions, ■01 on Ml and ip2 on M2, that do not agree on M3. It is then simple to 
construct neighborhoods (in either the compact open or Whitney fine topologies) 
of the graphs of these sections in V3 that do not intersect over at least one point 
in M3, which define neighborhoods of ipi and ip2 and hence a neighborhood of 
{tfji,'4'2) that does not intersect the image of T{V'). Therefore, the complement 
of the image of ^{V) is open. In other words, the image of T{V') in the product 
space TiVi) x r(y2) is closed. 

The above discussion will be sufficient to establish a closed embedding of 
the solution space S{M') in natural ambient space determined by the Cauchy 
pushout. From that, we can construct a surjective homomorphism onto the al- 
gebra A{M') from a natural algebra also determined by the Cauchy pushout. 
Unfortunately, that is not quite satisfactory, as it is irresistible to try to for- 
mulate the generalized Time Slice property as a stronger result, namely that 



This is sometimes known as the sheaf descent property of smooth sections. 
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the Cauchy pushout (in the spacetime category) is taken to a pullback (in the 
phase space category) and again a pushout (in the algebraic category) by the 
appropriate successive contravariant functors P and T (or S and A). However, 
to achieve that using an appeal to Hyp. |3] requires an additional transversality 
hypothesis on the interaction of the solution spaces Si^MC) and 5(7W). It is not 
clear whether this notion has already appeared in the standard PDE literature. 

Definition 50. A Cauchy pushout as in Eq. (|262l) is said to satisfy transverse 
descent if the smooth maps ai : S{Mi) x S{M2) S{Mi) S{M3), i ~ 1,2, 
are transverse to each other. 

Remark 7. Note that the images of the maps Xi '■ S{Mi) — ?> S{M3) intersect pre- 
cisely on a subspace of 5(A^3) that can be identified with S{M'). This property 
follows from the uniqueness of the Cauchy problem and the sheaf descent prop- 
erty inherited from the sheaf of smooth functions. As in the remark following 
Hyp. 13 we do not precisely specify the desired notion of transversality. However, 
we conjecture that such a notion can be formulated so that the maps ai agree 
precisely on the subset of S{Mi) x S{M2) identified with S{M') and, more- 
over, that a corresponding implicit function theorem establishes that S{A4') is 
in fact a submanifold. The reason the standard notions of infinite dimensional 
transversality [501I771I55] may not be applicable here is that (precisely due to 
the possible blow up of solutions) the images of the tangent maps Tai may not 
be closed subspaces of the tangent space TS{M3). This phenomenon appears 
already for linear PDE systems, such as, for example, the Cauchy-Riemann equa- 
tions. Though, admittedly, that is not an example of a hyperbolic PDE system. 

Lemma 23. Consider a Cauchy pushout diagram satisfying transverse descent, 
Defs. \^9\ and \5(K The following diagram contains a pullback square for S{Ai') 
and the canonical dotted line morphism is a closed embedding of manifolds: 

S{M3)i SiM2) 




(264) 



S{Mi) X S{M2) ■ 

If none of S{Mi) or S{M') are empty, then the above pullback square lifts to 
the following pullback square in &X)mp: 

V{M3)< V{M2) 

V(X2) (265) 

ViMi)^ ViM') . 

P(xi) 

Proof. The set theoretic pullback property holds again from the gluing prop- 
erty of solution sections, with which the background fields do not interfere. 
Now, recall that all section spaces are topologized as subspaces of r(— ) with 



86 



Igor Khavkine 



appropriate vector bundle argument. In particular this means that the subset 
S{Mi) X S{M2) n r{V') C T{Vi) x r(y2) is closed in S{Mi) x S{M2)- On the 
other hand, the gluing property ensures that S{M') = S{Mi) x 5(A^2) H r(F'). 
(Again, the obvious pushforwards and puUbacks have been omitted for brevity.) 
This shows that the canonical dotted line morphism is a closed embedding of 
topological spaces. This much we can establish without appeal to the transverse 
descent property, on the other hand we must appeal to it to establish that it is 
a manifold embedding, and hence a closed embedding. The covariance lemma, 
Lem. [221 guarantees that the morphisms lift to Sr)mp whenever the correspond- 
ing solution spaces are not empty. □ 

With the above discussion in mind, we can formulate the following general- 
ization of the Time Slice property for the functor 

Theorem 9 (Generalized Time Slice). Consider a Cauchy pushout diagram 
that satisfies transverse descent, Defs. \49\ and \50i Then the following is a pushout 
diagram in £2l[g.' 

AiMs) >AiM2) 

-A(X2) (266) 
AiM,)—^A{M') . 

-4(xi) 

Moreover, if none of the S{A4i) or S{A4') are empty, the above diagram also 
lifts to a pushout diagram in ^oiss: 

T{Mz) yT{M2) 

Hx2) (267) 
HMi)—^T{M') . 

Proof. According to Lem. [531 we can realize S{Ai') as a closed submanifold of 
S{Mi) X 5(7^2), with the aid of the transverse pullback diagram (|264p . On the 
other hand. Hyp. [3] immediately implies that this transverse manifold pullback 
is taken by to the algebraic pushout diagram (j266p . Finally, whenever the 
diagram (|264p in DJlan lifts to diagram (|265l) in ©tjmp, then the diagram (|266p 
in £2t[g automatically lifts to diagram (j267l) in *Poiss. □ 



7.3. Causality. Given two morphisms Xi - ■^i ~^ * = 1,2, in ©p^Btgr, do 
the corresponding morphisms surjectively factor through the product iS(A^i) x 
5(A^2)? Do the corresponding morphisms factor injectively through the tensor 
product A{A4i)^A{A42)'^ Finally, when do the corresponding morphisms factor 
through the independent subsystems tensor product J- {Mi) ® T{Ai2)- 

There are several obstacles to positive answers to the above questions. First, 
let us consider surjectivity for the spaces of solutions and, correspondingly, injec- 
tivity for the algebras. If we allow disconnected manifolds as base manifolds in 
*84gr and 6p*8tgr, while making sure that morphisms in 6p*8tgr treat discon- 
nected components as spacelike separated, we can introduce a tensor product. 
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referred to as the disjoint union: 

Mi\JM2 = (Ml U M2, Cf U Cf,Bi X B2). (268) 

Note that this is not a coproduct in the categorical sense, since dotted line 
morphisms in the diagram 

Ml >Mi M2 (269) 

M1UM2 

exists only if the images Xi(Mi) and X2{M2) are C*-spacelike separated in M. 
However, when this morphism exists, it is canonical and is denoted by Xi U X2- 
On the other hand, Tlan does have a categorical product, which means that the 
following diagram always exists: 

Ml (^^^^^^^^^ 1 y M2-^-—^ M2 (270) 
S{Mi) X S{M2) . 

Since there is no obstacle to specifying a solution independently in each con- 
nected component, the canonical dotted line morphism is in fact an isomorphism, 
S{Mi UM2) = S{Mi) X S{M2)- In the ©tjtnp category, we can also define the 
independent subsystems product given by 

{Ni,Qi) X (TVs, r?2) = (A^i X N2, Qi O2). (271) 

When Ni = S{Mi), the corresponding Poisson bivector 11 on the product is 
characterized by 11 (x, y) = ni{x, y) when x, y € Mi and 11 (x, y) = when x and 
y belong to different connected components. Note that this is not a categorical 
product in 6t)mp, for reasons very similar to why U is not a coproduct in ©p*Bfiflr. 
Ifowever, given the above information, we do have the identity V{Mi U M2) = 
V{Mi) X V{M2). 

Now, given the Sp^Btflt diagram 

Ml yMi M2 , (272) 

M1UM2 

we always obtain the diagram in ^)Ran below 

S{Mi)^^^ SiM) ^^"""^ ) SiM2) , (273) 

S{Mi) X 5(7^2) 
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where dotted lines denote the canonical morphisms. Note that the canonical 
morphism in OJlan always exists, even if xi U X2 does not exist in 6p*B6flr. 
The surjectivity of the canonical dotted line morphism in the last diagram is 
equivalent the following; given a pair of solution sections, ijji on A^i and -02 on 
M2, there exists a solution section ip on M that restricts to ipi on A^i and 
V'2 on M.2- Already in the earlier discussion of the Isotony property, we have 
noticed this surjectivity property fails when at least one of the Xi niorphisms 
fails to be extensible. However, canonical dotted line morphism may fail to be 
surjective even if both xi and X2 are extensible. Namely, there may exist a pair 
(V'i)V'2) € S{A4i) X <S(M2) such that both tpi and V2 may be extended to M. 
individually, but not jointly. 

The most common reason for that to happen is that the images xi(Mi) and 
X2(A#2) are C®-causally related in M, so that, for example, the value of tpi in 
Ml is not consistent with the causal influence of tp2 on M2. Thus, to even have a 
hope of canonical factorization, we should require the images xi U X2 does exist 
in ©pStgr, that is the images xi{Mi) and X2(Af2) are C*-spacelike separated in 
M. However, even with spacelikc separation, surjectivity can still fail. Intuitively, 
this happens when the solution data specified by Vi and V2 always produces a 
singularity while scattering in a joint extension to A4, so that the canonical 
morphism xi U X2 is itself not extensible. 

Definition 51. We say that two extensible morphisms Xi o-nd X2 in 6p5Slgt 
have regular scattering if they fit in the commutative diagram 



Ml >M< M2 , (274) 

Ml U M2 




such that the canonical dotted line morphism xi U X2 exists and is extensible. 

Once the surjectivity in StHan and (3t)mp is taken care of, the injectivity in the 
algebraic categories follow in the manner indicated in the earlier discussion of 
the Isotony property. With the above discussion in mind we can formulate the 
following generalization of the Causality property for the functor 

Theorem 10 (Generalized Causality). Consider two morphisms Xi - Mi 
M, i = 1,2, in ©pSfiflt such that xi U X2 exists with regular scattering. Then 
there exists an object M' = {M,C'®,B) in SpQSIflt such that Xi o,nd xi LI X2 
factor according to the diagram 
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the canonical dotted line morphism in the corresponding diagram in 9Jlan is 
surjective, 

S{M) 



S{Mi) i 




(276) 



S{Mi) X S{M 

and the canonical dotted line morphism in the corresponding diagram in £2l[0 is 
injective, 

A{M) 



A{Mi) 




AiM2 



(277) 



AiMi)®A{M2) 



Moreover, if none of the objects in diagram (j276p are empty, the diagram (|275l) 
lifts to ©pSJgr^, the diagram (|276p lifts to 6r)mp via V, and the diagram (|277p 
lifts to Cpoiss via T . 

Proof. Regular scattering hypothesis implies that both x^, i ~ 1,2, morphisras 
are extensible, using the same auxiliary object A^', where the images of x'i ^re 
also spacelike separated. The existence of the canonical dotted line morphism 
then follows from the definition of U. Applying the S functor to diagram (|275D 



■oo 
est 



produces diagram (|276l) . where the surjectivity of the dotted line morphism 
again, follows from the regular scattering hypothesis. Finally, applying the C, 
functor to diagram (|276|) gives diagram (|277|) . where the canonical dotted line 
morphism is injective because, by Hyp. [31 C^f. takes surjections to injective 
honiomorphisms . 

If the diagram (|275p can be lifted to ©pStgr^, then it follows from definitions 
and the covariance of the Poisson structures, Lem. [511 that the diagrams (|276p 
and (|277p lift to ©tjmp and *Poiss, respectively. □ 



7.4- Connection with standard LCFT. The above results can be applied to a 
semilinear, well-posed PDE system to construct a LCFT functor satisfying the 
standard axioms given in Def. [SI Consider a natural, variational, semilinear PDE 
system, whose background fields include a globally hyperbolic, time oriented 
Lorentzian metric g and such that its cone bundle of future directed timelike 
vectors Fm — V{M) coincides with the bundle of future directed timelike cones 
of g. Let S: ©pStgc — >• Man and V: ©pQStgc^ — > 6r)mp be its solutions and 
phase space functors. We will allow other background fields but only such that 
they influence neither the solution space nor the symplectic structure on it. The 
reason is that to achieve a natural symmetric hyperbolic form extra background 
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fields may need to be introduced, as has already been discussed and can be seen 
explicitly from the examples in 50|- However, since that is their only purpose, 
they do not influence the dynamics. We call such background fields auxiliary. 
For the purposes of comparing with a standard LCFT, we can safely ignore 
them. Recall that we introduced the category ©pStgr*, with objects that are 
endowed with background field collections consisting of a single section. Using 
this category to construct a standard LCFT is slightly problematic, because of 
the presence of auxiliary background fields. An easy solution is to introduce the 
forgetful functors ©p*84gr* — ^ Sp^BJgr' and Stgr* *B?0r' that throw away the 
auxiliary background fields, where ©p*B£gr' and *B{gr' now has only the metric 
as background field. The subscript c retains the same meaning as before. By 
stipulation, the functors S and V project to functors S : ©p*B5gr' — >■ VJlan and 
V: ©pSfgr^ — >■ ©t)mp, and similarly with Sh and Vh, up to isomorphisms in 
the target categories. Let a natural PDF system that satisfies all the conditions 
in the preceding discussion be called a semilinear geometric wave equation. 

We conclude this section by showing that the tools assembled so far in this 
paper allow us to construct a standard classical LCFT from a well-posed semi- 
linear geometric wave equation. Unfortunate, a sufficiently detailed notion of 
well-posedness that would allow us to directly construct a standard LCFT func- 
tor based on the results of the preceding section is somewhat difficult to for- 
mulate precisely. One obstacle is the fact the properties of the functor 
that constructs the algebra of observables as the algebra of smooth functions 
on an infinite dimensional phase space have so far been mostly hypothesized 
(cf. Hyp. [2] and Hyp. [3]) rather than constructively proven. Another is related to 
the boundary conditions that one is expected to impose on solutions in space- 
times with non-compact Cauchy surfaces (cf. Scct. l5.2Tl]) and their interplay with 
the standard notions of well-posedness in the PDF literature. As such, we make 
no special attempt to be precise and merely formulate the following hypothesis, 
whose application is made clear in the proof of the theorem below. 

Hypothesis 4 (Well-posedness). We consider our classical field theory de- 
fined by a natural, variational, semilinear PDE system to be globally well-posed 
in the following sense. All regular initial Cauchy data (where regularity includes 
any necessary boundary conditions) extend to the entire spacetime manifold (no 
finite time blow up). Moreover, in the presence of constraints and gauge sym- 
metry, the PDE system can be augmented with purely hyperbolic gauge fixing 
(Sects. [5.3.1\ and \5.3.2]) such that global parametrizability (Sect. V5'.2.2\) and global 
recognizability ( Sect. 1 5. 2. 3\) are satisfied. Finally, the functor constructs the 
algebras of observables in a way that is not sensitive to the behavior of solutions 
near the open ends of the spacetime manifold. 

To substantiate the above hypothesis, let us note that there exist global well- 
posedness results for what we consider fundamental bosonic fields (namely Yang- 
Mills 02 and non-linear scalar fields) other than gravity (GR definitely 
exhibits finite time blow up, as exemplified by the singularities of black hole and 
cosmological solutions). On the other hand, the algebras of observables consisting 
of microlocal functionals of compact spacetime support [Sl IMllTfll^ provide a 
candidate for the functor C^j, which seems to have the desired properties, at 
least when applied to the solution spaces of linear theories. Also, the sufficient 
conditions imposed on the constraints and gauge transformations are actually 
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very similar to those considered in previous treatments, as exemplified in (50) 
and 158]. 




Theorem 11. Consider a semilinear geometric wave equation that is well-posed 
(Hyp.^. Then we have an equivalence of categories Stgr^ = ©lob^^OPc '^'^^ 
corresponding algebra of observables functor J-r '■ ^Sfigr^ — >■ *Poiss satisfies the 
standard axioms of classical LCFT given in Def. O 

Proof. The equivalence of categories ^tgr'^ = ©[obiot^p^ is straightforward to 
establish. The objects are identical and consist of pairs (M, g) (or equivalently 
(M, {(7})), where 5 is a globally hyperbolic Lorentzian metric on the manifold 
M. There exists a morphism {M,g) — >■ {M',g') in *B£grJ, iff there is a morphism 
M ^ M' in 6pS{0r, between objects that project to n(A^) = {M,g) and 
n(A^') = {M',g') and neither solution patch S{M) or S{M') is empty. For 
semilinear equations, a solution (/) € S{A4) iff the corresponding spacelike cone 
bundle C® — )■ M is strictly faster than the bundle of future oriented Lorentzian 
cones. Now, there exists a morphism A4 — > M' only if the underlying open 
embedding M — )■ M' satisfies x*9' = 9 (is an isometry) and is chronally 
compatible with respect to C® and C®, which a fortiori implies that x is also 
causally compatible with respect to g and g'. This concludes the proof of the 
equivalence of the two categories. 

Next, we establish a crucial property that helps us prove the remaining con- 
clusions. For any object Ai with n(A^) = {M,g) and spacelike cone bundle 
C® — >■ M that is faster than the bundle of future oriented Lorentzian cones, we 
have the identity SH{M,g) = S{M). That is, all solutions are already C®-slow. 
This tells us that the image of the functor S^{M,g) in dJlan, consisting of all 
objects and morphism S{M) -> S{M') with n(A^) = U{M') = {M,g), where 
each S{A4) is either empty (if the cone bundle C® is slower than the Lorentzian 
cones) or S{A4) = Sh{M, g) (if the cone bundle C® is faster than the Lorentzian 
one). The same goes for the respective objects in other categories, P(A^), A{M), 
and J^{Ai), except that the P(A^) are never empty. In other words, the limit 
construction FniM^g) = lim^(M, 5) becomes somewhat superfluous, since the 
image of the diagram JF(M,g) already contains the object that is equal to the 
limit. 

With the above simplification, which means that we can simply replace the 
Poisson algebra FniM^g) with the algebra ^(A^) for some appropriate A^, 
the diagrams defining the generalized Isotony (Thm. |5]), Time Slice (Thni. [S]) 
and Causality (Thm. ITUl) properties collapse to the usual notions thereof as 
stated in the standard definition of LCFT, Def. [5j It remains only to show 
that the key hypotheses needed for the above theorems, namely extensibility of 
morphisms in (5p*B£gr^, transverse descent for Cauchy pushouts in ©p^BJgr^., 
and regular scattering for morphisms in Qp'^tQt^, all follow from global well- 
posedness, Hyp. H) 

Given a morphism — > A^', extensibility fPef. |47)) implies Isotony. Ex- 
tensibility almost follows from well-posedness. Namely, if there exist Cauchy 
surfaces S <Z M and Z" C M' such that S' extends S and we have initial data 
ip on S such that its pushforward x*'f smoothly extends to some initial data 
<f' on Z", then the corresponding solution extends from x(M) to A/'. However, 
even if Conj.[l] holds, the best we can expect is the existence of a Cauchy surface 
S' C M' that extends a compact subset K CI S. Also, even if we can extend all 
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of S to U', the initial data could behave near an open end of in a way that is 
not smoothly extendible to E'. So to recover Isotony at the algebraic level, we 
must appeal to the part of Hyp. 2] according to which the algebra of observables 
constructed by applying is not sensitive to the behavior of solutions in the 
neighborhood of the open ends of M and, in particular, near the open ends of 
E. In that case, using Conj.[T]and global well-posedness as above, every smooth 
solution (j) can be extended from an arbitrarily large compact subset K C S 
to all of M. The details of this argument would have to await a more precise 
formulation of HI and [1] and of the choices discussed in Sect. I5.2TT] 

Given a Cauchy surjective morphism x- -M Ai' (Def. HB]) . it could always 
be completed to a Cauchy pushout for A4' (Def. [351) ■ Given global well-posedness, 
no solution blows up in finite time. This means that all the Cauchy surjections, 
including induce isomorphisms in the diagrams (j265p and (I267p . In other words, 
transverse descent holds trivially. The generalized Time Slice property then re- 
duces to its standard version from Def. [5] 

Finally, since global well-posedness (in particular the absence of any finite 
time blow up) also implies regular scattering, Def. [5T1 the generalized Causality 
property holds, Thm. fTUl and reduces to the corresponding standard notion from 
Def. El □ 

The theorems presented in this section may be seen as the translation of the 
well-posedness properties of a natural variational PDE system into the algebraic 
setting given by the corresponding Poisson algebras of observables. One may 
even promote some of these properties to an axiom system generalizing the 
existing LCFT axioms, which are only applicable to semilinear wave systems, to 
the more general class of quasilinear systems, which includes GR. We do not do 
so immediately because several aspects of the current formalism, to be discussed 
in Sect. (HI leave room for substantial improvement, which could facilitate an 
optimal axiomatization. 

8. Causality Axiom in Locally Covariant Quantum Field Theory 

This section contains some remarks on extending the results presented so far 
to quantum field theory. The remarks essential consist of the conjecture that 
the categorical colimit construction of the algebra of observables will persist 
through perturbative deformation quantization of the classical field theory. Since 
rigorous, non-perturbative construction of quantum field theory is in general not 
yet possible, that is the best we can hope for at the moment. 

Recall that a classical mechanical system essentially consists of a real Pois- 
son algebra of observables J" = (^, {}), where A = C^((P) on a symplectic 
manifold P (the phase space), whose states are normed positive linear func- 
tionals on A, which can be expressed as probability measures on P. On the 
other hand, a quantum mechanical system essentially consists of a complex, 
associate, non-commutative *-algebra J- = (^, *,*), whose states are normed 
positive linear functionals on A. Strong physical and mathematical arguments, 
which were already eloquently expressed in the original papers |lllll2lfT5] , indi- 
cate that deformation quantization is the right way to define an ?i-parametrized 
family of quantum systems that quantize a given classical one. Briefly, a de- 
formation quantization of a classical mechanical system (A, {}) is a family of 
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quantum mechanical systems (^^, *"), where the underlying vector spaces 
are isomorphic to C (8) ^ = An, the ^-involution reduces to standard complex 
conjugation, A*'^ — ;> A, in the limit ?i — > 0, and the non-commutative product 
commutator satisfies the classical limit, ^[^, i?]ri ^ {A, B} as ?i — > 0. The non- 
commutative product *?i is referred to as the -product. The details of the kind 
of dependence on h is allowed and the sense in which the limits are taken can 
be found in the standard literature |76j . 

From the constructive point of view, deformation quantization has a num- 
ber of successes to its name. The Poisson algebra of any symplectic [37,45l|31] 
and even any Poisson manifold [74j possesses a formal deformation quantiza- 
tion. The deformations are rigid and physically inequivalent classes of deforma- 
tions are controlled by the second de Rham cohomology group of the underlying 
manifold |55|. Constructions of formal deformation quantization commute with 
reduction by quotienting out gauge symmetries, via the BV-BRST method, ob- 
structed only by well defined anomalies P^HTTll^ . Formal deformations are 
possible also for symplectic supermanifolds and infinite dimensional symplectic 
manifolds jl05) (and references therein). The formal deformations of standard 
M^" symplectic spaces can be made strict via the Wigner-Weyl-Moyal ^-product 
formula [76^. Similar results are available for some other special classes of finite 
dimensional symplectic manifolds |98j . 

As discussed in Sect. [H a classical field theory assigns classical mechanical 
systems to spacetimes in a way coherent over spacetime embeddings, tentatively 
summarized in the axioms for a locally covariant field theory (LCFT) functor, 
Def. [5l and similarly for a quantum field theory (QFT or LCQFT). By a quan- 
tization of a classical field theory, we mean an LCQFT functor that reduces 
to a classical LCFT functor in the classical limit ?i — > 0, in the sense of defor- 
mation quantization, also in a way coherent with spacetime embeddings. The 
recent literature on pAQFT deals specifically with the perturbative quantiza- 
tion of classical field theory in the above sense. From the relevant literature, it is 
clear that the formal deformation quantization point of view is compatible with 
perturbative renormalization of quantum field theories [39l[64ll47lf95] . 

Note, however, that the field theories, whose quantization has been studied in 
the above way, have only been of semilinear type or whose non-linearities have 
been treated perturbatively as well. Thus, their Green functions (of appropriately 
linearized equations) and their Poisson brackets (via the Peierls formula) have 
causal supports fixed by external background fields (essentially a Lorentzian 
metric) , as reflected in the Causality axiom of Def. [5] The quantum version of 
the Causality axiom uses ^-product commutators, which obey a similar causal 
support condition. The precise form of this condition has been known for a long 
time (in fact going back to the original axioms of Haag and Kastler) and has been 
checked to hold at each perturbative ?i-order of the deformation quantization |39l 

On the other hand, the situation for quasilinear flcld theories, which have a 
dynamical, field-dependent causal structure (GR being a prominent example), 
has been much less clear. The classical notion of causality, relying on the causal 
support of the Green functions of the linearized equations about a fixed back- 
ground solution (hence, via the Peierls formula, also of the Poisson brackets), 
is unproblematic and has been understood for a long time. However, it relies 
the notion of a fixed background solution or, in other words, of a fixed point in 
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the classical phase space. Since the standard formalism of quantum mechanics 
does not make use of a phase space, the question of translation of the notion of 
causality to the quantization of a quasilinear field theory has been generally con- 
sidered open, especially outside the perturbative context, where it is sometimes 
referred to as quantum fluctuation/smearing of light cones [721181] . 

In this paper, we have attempted to precisely formulate the analog of the stan- 
dard Causality axiom for quasilinear classical field theories. This formulation has 
been given in terms of localizing the algebra of observables to open subsets of 
the phase space (the space of all solutions modulo gauge transformations) whose 
elements (individual solutions) all define causal structures that are compatible 
with some externally specified chronal cone bundle (in a way that is standard 
for symmetric hyperbolic, quasilinear PDE systems. Sects. [3] and A cover of 
the total phase space by such slow patches, can be interpreted in categorical 
terms as the identification of the total phase space with the pushout of a dia- 
gram in the category of symplectic manifolds corresponding to the slow patches 
and the intersections between them. This categorical formulation makes it obvi- 
ous that the corresponding algebraic formulation of the analog of the Causality 
axiom is the identification of the total algebra of observables with a limit of the 
diagram in the category of Poisson algebras corresponding to the algebras of 
observables of the slow phase space patches, with the Poisson bracket respecting 
the causal structure of the algebra of observables of each individual patch. The 
details of this formulation occupy Sect. l7.3l and culminate in Thm.fTOl Similarly, 
an algebraic formulation is also given to the analog of the Time Slice axiom, 
Thm. [5] in Sect. 17.21 An important check on these generalized Causality and 
Time Slice properties is that they reduce to the corresponding standard notions 
when specialized to semilinear field theories, Thm. 1111 

Finally, the above purely algebraic formulation of a generalized Causality 
property of classical LCFT motivates the following conjecture for an LCQFT 
deformation quantization of a classical LCFT: 

Conjecture 3. Let *-2t(g be the category of non-commutative, associative, *- 
algebras, while Stgrj, and SpQStgr^ are the categories of spacetime manifolds 
endowed admissible background fields and, respectively, also with a globally 
hyperbolic spacelike cone bundles (as in Sect. 16.21) . An generalized LCQFT 
is a covariant functor J^h '■ ^tfl'^c ~^ *-2ll0, such that there exists a functor 
JF: ©pSIgr^ — > fitting into the following commutative diagram 




(278) 



and satisfying the identity Fh{M,B) — ^im^(Af, S), where the limit is indexed 
by the subcategory &p'^lQx:^{M,B) C ©pSSgr^, whose objects are of the form 
M = {M,C®,B). This LCQFT also satisfies the following Causality property. 
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The existence of the diagram 




Ml 



M1UM2 

in ©p^Btgr^ iinphes the existence of the diagram 

U 



(279) 



(280) 



P{Mi)®P{M2) 

in *-2l[g. Recall that U denotes the spacelike separated disjoint union (Sect. [721), 
while ® denotes the independent subsystems tensor product (Def. 

Similarly, the Time Slice property can be straightforwardly translated to 
LCQFT as well. Though the above formulation is rather abstract, we can see how 
it looks in rather concrete terms in the example of GR. Consider a gauge fixed 
version of GR. It may include fermionic ghost fields, if necessary, though whose 
properties we do not discuss at the moment (that is left to future work [70]). 
At the moment, we also ignore the issue of difFeomorphism invariant observables 
and simply consider components of the quantized, gauge fixed metric field gab{x) 
as observable. That is, once smeared, they constitute elements of the algebra 
J^h{M, {/?}), for a spacetime manifold M endowed with admissible background 
field /3. What can we say about the ^-product commutator [gab{x), (jcdiy)]'^ As is, 
we cannot say very much, since we do not have complete solution for the theory 
and due to dynamical causal structure, we cannot appeal to the standard Causal- 
ity axiom. However, since J"jy(M,{/3}) = lim J"(7W), with M = (M,C®, {/?}), 

there exist canonical projections ^//(M, {/3}) — >• T{A4) for each such M. De- 
note the images of the metric fields under these projections by gab{x) 1— > 3^(2;). 
Then, if we select a spacelike cone bundle C® — M such that two given points 
x,y G M are C®-spacelike separated, our conjecture about the Causality prop- 
erty in LCQFT implies that 



-0, 



(281) 



in the usual distributional sense. 

As we have seen, the Causality property is a theorem of classical field the- 
ory. On the other hand, its version in quantum field theory is only a conjecture. 
Unfortunately, we cannot do better at the moment, due to the general difficulty 
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of non-perturbative quantization of physically relevant field theories ^S] . More- 
over, even in semilinear field theories, the quantum Causality axiom has only 
been verified in free field theories, in four dimensions. On the other hand, as 
discussed above, it has been verified perturbatively in standard physics calcu- 
lations as well as in the pAQFT literature. For interacting field theories, the 
perturbative construction involves a double formal series, in h and in A (the 
collection of all coupling constants). Work is currently under way to verify it 
(as part of a more general quantization program) at a level that is perturbative 
in h but non perturbative in A, but extending Fedosov's deformation quantiza- 
tion construction to LCFT. It is likely that the above conjectured generalized 
Causality property for quasilinear field theories may be verified explicitly in the 
same setting The strategy of these approaches is, as in this paper, is to 

use PDE theory to explicitly construct the interacting (nonlinear) classical field 
theories before applying Fedosov's deformation quantization construction to the 
resulting phase spaces. 

9. Discussion 

In this work, we have in a sense taken up a thread left loose in the previous 
work of Geroch [50] on the application of the PDE theory of symmetric hyper- 
bolic systems to classical field theory. This theory provides an unambiguous, 
intrinsic notion of causality not tied a priori to a Lorentzian metric. This notion 
of causality is especially useful for theories with quasilinear (rather than linear 
or semilinear) equations of motions (of which General Relativity is a prominent 
example), whose causal structures are field dependent. This is discussed in detail 
in Sects. [3] and m See also [51 [M[[551[5T| for related earlier ideas. 

Moreover, it provides the theorems cited in Sect. 14.3] as constructive tools to 
build the solution space of the theory as an infinite dimensional manifold (the 
infinite dimensional differential geometry aspect, as discussed in Set. 15.21 was 
treated only formally here, though is tackled in earnest in the recent works [47l 
inSJUl]). Geroch's casting of most classical relativistic field theories in symmet- 
ric hyperbolic form was then supplemented by constructing (formal) symplectic 
and Poisson structures fSect. [575)) . turning the space of solutions into the genuine 
phase space of classical field theory. This was done by making use of the well 
known covariant phase space method and the Peierls formula, that was gener- 
alized to treat field theories that possess gauge invariance and constraints when 
cast into symmetric hyperbolic form. 

Together with some reasonable technical hypotheses on the structure of these 
infinite dimensional phase spaces and algebras of smooth functions on them 
(HypslU ID and 131), the generalized Isotony, Time Slice and Causality proper- 
ties were established in Sect. [71 As an important check, these notions specialize 
(Thm. [TTl) to the synonymous axioms of classical relativistic field theory (Def . O . 
In formulating the generalized properties of these axioms, which are essentially 
algebraic in nature, we identified strong parallels between them and aspects of 
the notion of well-posedness in classical PDE theory, which are more geometric 
in nature. In particular, we have highlighted the importance of the on-shell or 
phase space point of view that is dual to the off-shell point of view adopted 
in the recent works [T7ll95[[^ . While the off-shell approach has proven invalu- 
able in the modern formulation of perturbative renormalization of quantum field 
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theories [S^I^HTll^ , the on-sheh approach together with the above mentioned 
results from PDE theory is currently the only constructive method of building 
non-perturbatively interacting classical field theories. The parallels between the 
axioms of classical relativistic field theory and well-posedness also show that it 
is sometimes too much to hope for the axioms to hold in every physically rea- 
sonable example, especially in the presence of finite time blow up singularities 
(which occur in GR, for example). In particular the standard Isotony property 
fails in the presence of finite time blow (failure of global well-posedness) up 
(Sect. 17.11) . On the other hand, local well-posedness is sufficient for the Time 
Slice property to hold, unless a more technical condition of transverse descent 
is also violated (Sect. [721) ■ the other hand, the Causality property holds as 
long as a weaker version of global well-posedness holds that we called regular 
scattering (Sect. [73]) . 

The formulation of the generalized Causality property relies on the key result 
Thm. [3 which could be called a generalized Microcausality property and follows 
from the Peierls formula for the Poisson bracket (Sect. I5.3.4p and the domain of 
dependence theorem for hyperbolic PDE systems (Cor.[T]). This connection can 
be seen as one of the main results of this work when coupled with the follow- 
ing speculation. If a quantum field theory is constructed from the deformation 
quantization of a classical field theory, then the Poisson algebra formulation 
of classical causality captured by the generalized Causality property should be 
replaced by a non-commutative *-algebra version of the same (Conj.lS]). If defor- 
mation quantization is carried out using Fedosov's method, it is expected that 
the conjectured formulation of quantum causality at the very least holds for- 
mally at each order in h ^71155] . Finally, if this conjecture holds, it provides a 
concrete answer (Eq. (j281|) ) to this old question in quantum gravity: What can 
we say about the commutator of two metric field operators in quantum GR? 

One of the motivations for this work had been to translate the notion of well- 
posedness from PDE theory to the algebraic setting for classical field theory, 
so that this translation could serve as a basis of a more realistic axiomatization 
thereof (see preceding paragraph for a discussion of why the axioms in Def. [S]are 
not completely adequate). However, such a translation attempted in this paper 
has lead to a number of necessary but somewhat tangential ideas, conjectures 
and hypotheses. We believe that it is essential to address and clarify some of 
them before attempting to formulate a more realistic axiomatization. 

For example, we have chosen the notion of symmetric hyperbolicity because 
of the large existing literature on this subject, making it easy to use to leverage 
the existing results of PDE theory in the non-perturbative construction of clas- 
sical field theories (Sec. l4.3| ). However, reducing the Euler-Lagrange equations of 
a field theory to (constrained) symmetric hyperbolic form is not always an ob- 
vious task, and, even if so, can be rather laborious and require the introduction 
of auxiliary background fields to keep everything natural. On the other hand, 
the notion of regular hyperbolicity |31) seems to be much better adapted to La- 
grangian field theories, but is not yet currently general enough to handle all cases 
of interest (such as for example higher order theories) and has a substantially 
smaller literature devoted to it. However, since both notions rely (underneath 
the hood of their respective well-posedness theorems) on so-called energy meth- 
ods, it is likely that it would be possible to subsume both notions under a more 
general one that uses intrinsic information about the characteristic cohomology 
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of a PDE system [2511511^ to identify approximate conservation laws and auto- 
mate their use in an energy method that would establish local well-posedness. 
Pursuing such a comprehensive notion of hyperbolicity is a promising avenue of 
investigation. 

Another set of ideas naturally follows the notion of a conal manifold, some 
of which have already been explored in Sect. 01 Conal manifolds abstract the 
notion of causal order from Lorentzian geometry, but remain within the realm 
of differential geometry. It is interesting to see which results generalize from 
Lorentzian to conal manifolds. One important result that already exists is the 
splitting theorem for globally hyperbolic cone bundles (Prop. IS]) • This result is es- 
sentially differential topological in nature, very reminiscent of results that have 
been obtained as instances of the so-called h-principle [M1I103) . On the other 
hand, the known proofs of the classic and the (very recent) generalized [43] 
result use methods rather removed from differential topology (in the former case 
appealing to non-continuous measures and in the latter to Weak KAM theory). 
A notable variation on the classic result uses purely order-theoretic and topo- 
logical methods [85]. It would be very interesting to formulate a proof using 
the methods of differential topology. In particular, such methods could then be 
adapted to the PL (piecewise-linear) and topological manifolds. The latter are 
of interest, for example, in the study of the causal structure of rough Lorentzian 
metrics [37. Another important set of results that should be generalized con- 
cern the addition of boundaries consistent with causal structure (in particular, 
causal compactifications) [15] . There is ample evidence from conformal compact- 
ifications in relativity that understanding the structure of a causal boundary is 
an important ingredient in understanding the boundary conditions and asymp- 
totic behavior of solutions to hyperbolic equations. On a more elementary side, 
it seems fruitful to investigate a version of de Rham cohomology with space- 
like compact supports, as well as the dual homology, whose classes seem to be 
naturally represented by Cauchy surfaces. Such theories may give us a better un- 
derstanding of the topological and order-theoretic ambiguities in the covariant 
phase space method (Sect. [^331) and even give us a better understanding of the 
structure of the space of non-globally hyperbolic solutions of a hyperbolic PDE 
system. Finally, it is important to verify the validity of Conjs. [2] and \T\ given 
their importance for the topological structure of the phase space (Sect. 15.2. ip 
and for the generalized Isotony property (Sect. 17.11) . 

The refinement of the notion of spacetime support to that of local spacetime 
support (Defs. [321 and HO)) invites us to reconsider the problem of local observ- 
ables in GR. Without this refinement, it is a well known fact that there are 
no gauge (diffeomorphism) invariant observables in GR with compact spacetime 
support [47[l95j . On the other hand, the possibility of gauge invariant observables 
with globally compact local spacetime support is yet to be considered in detail. 

There is still a substantial number of not-completely-resolved questions on the 
topology and differential geometry of the infinite dimensional space of solutions 
of a hyperbolic PDE system and the algebra of smooth functions on it. Given 
the geometric, rather than the functional analytical focus of this work, we have 
avoided most of such details and instead have simply made several reasonable 
hypotheses (Hyp. [1] [2] and [3]). To complete the program of constructive classical 
field theory, it is crucial to identify the precise theorems that would replace these 
hypotheses. 
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Finally, note that our discussion of classical field theory has been restricted to 
bosonic fields. On the other hand, to take fermionic matter into account as well as 
ghost fields that feature in the BV-BRST formulation of gauge theories, we also 
need to consider fermi fields. We believe that the best way to do that is shift from 
the setting of manifolds to supermanifolds |101] , so that fermi fields are sections 
of odd vector bundles over spacetime (which is equivalent to more common 
but less precise phrasing that classical fermi fields are Grassmann valued) . This 
approach is more geometric than and complementary to the algebraic approach 
adopted in [Mlf551H7[l95) . The fermionic field theories considered so far in this 
setting have only been of semilinear type. An extension of this formulation to 
also encompass quasilinear field theories (along the lines of the older work [291) 
and their dynamical causal structure will be reported elsewhere |70| . 
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A. Jet bundles and the variational bicomplex 

In this section we briefly introduce jet bundles and fix the relevant notation. For 
simplicity, we restrict ourselves to fields taking values in vector bundles. However, 
the discussion could be straightforwardly generalized to general smooth bundles. 

We briefly introduce fc-jets, mostly to recall some basic facts and fix no- 
tation. More details, as well as a coordinate independent definition, can be 
found in the standard literature P51I73[[TU5] . Fix a vector bundle F — Af, with 
dim A/ = n, with fibers modeled on a vector space U, and consider an adapted 
coordinate patch M" x U, with coordinates (x^w"). Extend this patch to a k- 
jet patch x U x [/"'' by adding extra copies of U, with new coordinates 
(a;% u", uf, ufj,. . . , w° ...j^), which formally denote the derivatives of di^i^.-.tp'^^x) 
of a section (j) at x. To keep track of all the derivatives, we introduce multi-index 
notation. A multi-index / = iii2 ■ ■ - ik replaces the corresponding set of sym- 
metric covariant coordinate indices (the multi-index does not change when the 
defining i's are permuted). The order of this multi-index is given by |/| = k, 
with |0| — 0. To augment a multi-index by adding another index, we use the 
notation Ij — jl — ii---ikj- Thus we can write higher order derivatives as 
di^...i^(j){x) = di<j){x), the higher order jet coordinates as uf^...^^ — Uj and the 
total set of coordinates on a fc-jet patch as {x\uj), \I\ < k. In particular the 
empty multi- index 1 — corresponds to — u°'. 

Since the higher derivatives are symmetric in all indices, the number of extra 
coordinates is given by rik = X]f=i dim S''''K", with S'^ denoting the symmet- 
ric tensor product. Given two different coordinate patches on F, we define the 
transition maps between the corresponding fc-jet patches according to the usual 
calculus chain rule applied to higher order derivatives. These fc-jet patches can 
be glued together into the total space of the k-jet bundle J^F — ^ Af, which 
includes J°F ^ F. 
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Since F —i' M is a vector bundle, so is J'^F — > M. It is isomorphic to F (Bm 
{F®mS^T*M)®m- ■ ■®m{F®mS^T*M), but not naturall£3. Jet bundles come 
with natural projections J^F J^~^F, which simply discard all derivatives of 
order k. This projection gives J^F the structure of an afRne bundle over the base 
J^-'^F, with fibers modeled on the vector bundle {F (S>m S''T* Mf''^ J^-'^F. 
The bundle J^F J'^^^F is afFine because, in general, bundle morphisms of 
jkp _^ jkp jijjjuced by vector bundle automorphisms of F are not linear, but 
are affine. 

Given a vector bundle E M it can be pulled back to the fc-jet bundle along 
the projection J^F — > M . We introduce a convenient notation for this puUback. 

Definition 52. We denote by {Ef J^'F the pullback of E M to J'^F, 
which then fits into the pullback commutative square 

{Ef > E 

(282) 

Jkp ^ M ^ 

Any smooth section 0: M F automatically gives rise to its k-jet prolonga- 
tion or k -prolongation j'^cp: M — > j'^F. Namely j'^cj) is a section of the bundle 
Jkp ^ M that is defined in a local adapted coordinate patch as 

jV(a^) = = {x\dicj,'^{x)), \I\ < k. (283) 

One can think of the fc-prolongation symbol as a differential operator 

/ : r(i^) ^ r( J'^F) (284) 

of order k. In fact, any (not necessarily linear) differential operator of order k, 

f:TiF)^TiE), f:<P^f[<P], (285) 

can be written as a composition of j'^' with an order (not necessarily linear) op- 
erator /: J^F ^ E, such that /[(/i] — f{i^(t>). Note that we are slightly abusing 
notation by denoting both the differential operator and the bundle morphism by 
the same symbol /. 

Further, we can define an Z-prolongation of a differential operator / of order 

p^j-jk+ip^jip^ (286) 

which is then a differential operator of order k + I, by composing with j': 
P°°f[4'] = Prolongation is discussed briefly using coordinate-wise op- 

erations in Sect. 13.31 The fc-jet prolongation can now be thought of as a 
special case of bundle morphisms, that is, — p'^cj), where on the right hand 

Though both the fc-jet and the "direct sum" bundles can be constructed by applying a 
functor QJ33na( — > QJIBnOl to a bundle F — > M, a vector bundle automorphism x- ^ ^ ^ 
induces a vector bundle automorphism J*(x): J^'F — >■ J^'F that need not be block diagonal 
in a basis adapted to the "direct sum" bundle, while the bundle automorphism induced by 
the "direct sum" bundle is. Therefore, the isomorphism of the "direct sum" and fc-jet bundles 
cannot depend on the object F ^ M alone, and hence cannot bo chosen naturally. 
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side we interpret (f> as the base fixing bundle morphism to F ^ M from the 
trivial 0-dimensional bundle id : M — >■ M. 

M — ^ F (287) 

id 

M >M . 

Given the sequence of projections fc-jet bundles over M, 

> J'^F ^ J^F ^ J^F^ F, (288) 

it is convenient to introduce the infinite jet order (or oo-jet) bundle J°°F defined 
as the projective limit over the jet order k 

J°°F = '\^J''F. (289) 

This limit implicitly defines J°°F as an infinite dimensional smooth manifold. 
The main advantage of working with oo-jets is that any function or tensor on 
j'^F for finite k can be pulled back to J°°F. Conversely, any smooth function or 
tensor on J°°F depends only on jets up to some finite order, say k, and can be 
faithfully projected to J^F. Another major convenience of working on J°°F is 
the ability to decompose the usual de Rham differential into its horizontal and 
vertical parts 

d = dh + dv. (290) 

The defining property of dh is the following. Given a section <p: M ^ F, we 
must have the identity 

(j°°0)*dha = d(j°°0)*a, (291) 

where a is any differential form on J°°F and d is the usual de Rham differential 

on M. On the other hand, dy is characterized by the fact that its image is 
annihilated by the puUback to M along any section <p, 

{]°°^)*d,a = 0. (292) 

It can be checked that the horizontal and vertical differentials anti-commute and 
are separately nilpotent, 

d^ = = d^, dhdv + dvdh = 0. (293) 

Note that, to apply dv or dh to forms defined on a finite order jet bundle J'^F, 
the puUback and projection operations mentioned above will often be applied 
implicitly. Thus the application of say dh to a differential form on J'^F may yield 
that a differential form that projects to J^'^^F but not to J^F. In local coordi- 
nates (a:%u°) on F, and the induced coordinates (a;',Uj) on J°°F, a convenient 
basis for differential forms is 

dhx' = dx\ dvu? = dw? - d^u] = dw? - u%dx\ (294) 

We can also define two special kinds of vector fields. A vector field ^ is horizontal 
if its action in local coordinates is 



(295) 
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for some ^* = ^*(a;,w°). In particular, the vector field dj, with = Sj, is hor- 
izontal. Note that [di,dj] = 0. A vector field ■0 is evolutionary if its action in 
local coordinates is 

7A(x^) = 0, ^j{u])=dj{r), (296) 

for some i/j'^ = '4)°-{x, Uj), where = di^ {di^ (• • • di^ (/) • • • )) for multi-index 

/ = «i«2 (the order of application of these vector fields does not matter since 
they commute) . Note that the tp"" can be seen as the fiber coordinate components 
of a section of the bundle (F)°° — J°°F. These definitions can be checked to be 
coordinate independent. 

One can show that for a horizontal vector field ^ on J°°F there exists a vector 
field ^0 on M such that their actions on scalar functions are intertwined by the 
puUback along the jet prolongation of a section Af — > F, 

e(/)(j°°0) = e0(/(r0)), (297) 

for any scalar function / on J°°F. Namely, in local coordinates, — ^^9^ with 
= (t|da;*)(j°°0) — £,{x^){j°°<p) — On the other hand, evolutionary 

vector fields tp satisfy the identities 



i^(dha)+dh(i^a) = 0, (298) 

jo-Arcpya 



[j°°i^ + e^)]*a^i^d,a^C^a, (299) 

e=0 



de 



for any form a G il,*{J°°F) and section ip: M F. Actually, ip could be a 
section of (F)*"' -> J'^F, that is, it could depend on 0°(x) and its derivatives 
and not only on x € M. The only corresponding change in the above formula 
would be to replace £-0 by e{j''(t))*ip. Ostensibly, should stand for the Lie 
derivative on the infinite dimensional manifold of sections of -F — > M, where the 
section ip is identified with the vector field whose action on local coordinates is 
Ctp4>°'{x) = ip°'{x). However, since we do not delve into the differential geometry 
of infinite dimensional manifolds here, we keep the symbol C^{j°°4))* primitive 
and defined as above. 

Integrations or differentiations by parts are carried out using the following 
basic identity 

d^u^ji A dx' A a ^ d^{u%dx') A a (300) 
= (dvdhM?) A a (301) 
= -(dhdvu?) Aa (302) 
= -dvU? Adha-dh(dvU?Aa). (303) 

This split of the de Rham differential into horizontal and vertical differen- 
tials also splits the de Rham complex rt*{J°°F) of differential forms on J°°F 
into a hicomplex [HII]. Since the horizontal and vertical 1-forms generate the 
graded commutative algebra of differential forms, any form A € ^*{J°°F) can 
be uniquely written as 

A = ^Aft,„, (304) 
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where < h < n and < v are respectively the horizontal and vertical de- 
grees form degrees. We have thus turned the differential forms into a bigraded 
complex r2*(J°°F) — 0,^ „ with the dh differential increasing /i by 1 

and the dv differential increasing u by 1. This complex is called the variational 
bicomplex As with any bicomplex, we can consider its cohomology with 

respect to either or any combination of the two differentials. The horizontal 
cohomology is -ff''''"(dh) = H {il* {J°° F) , du) in degrees {h,v). The vertical co- 
homology is i?^'"(d„) = H{n*{J°°F),d^) in degrees {h,v). Both (i7'''*(dh), dv) 
and (-ff*'''(dv), dh) still form complexes, therefor we can also consider their co- 
homologies. The relative cohomologies are H^ *{dy\A\^) = H{H^'* {du),dv) and 
H*'^{du\d,)^H{H*'^{d,),du). 

B. Limits and colimits in category theory 

Some basic information on category theory can be found in |21| and more specif- 
ically about categorical limits and colimits in [2111110) . 

Consider categories 3 and £ and a functor S: 3 —>■ The image S{3) (also 
denoted by just S) forms a diagram in £ with index category 3. A cone to the 
diagram S is an object C of £, the vertex, together with a morphism ci: C ^ 
Di = S{i) for each object i of 3, such that diagram 



C 




(305) 



> 

is commutative for each morphism i — j in 3. If it exists, the limit (also inverse 
or projective limit) is an object lim ^ of £ that is the vertex of a cone of canonical 
morphisms Ui : ^m S ~> Di such that the following universal property holds: 
any cone (C, c) to S factors through it with a unique mediating morphism (also 
canonical morphism) uc, which makes the following diagram commute: 



C-uc^^^^D, . (306) 

For an illustrative example, consider the index category 3 consisting of only the 
objects and morphisms 1 — > 3 ^ 2, while ^ is a functor from 3 to the category 
of sets. The diagram the limit and the cone of canonical morphisms fit into 
the following commutative diagram 

lim^ ) D2 

(307) 
Di > . 

The universality of the limit identifies the set ^im S with the subset of the 
Cartesian product Di x D2 consisting of pairs ((^1,^2) that get mapped to the 
same element in D^, di 1— > da ^ ^2- If the object 3 and the morphisms to it 
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wore absent from the index category 3, we would simply have ^im ^ = Di x D2. 
In other categories, we may keep the same intuition, but replace the Cartesian 
product X by the categorical product. 

The notion of a colimit is dual. As before, consider a diagram : 3 — > C A co- 
cone of the diagram ^ is an object C of £, the vertex, together with a morphism 
Ci : ^{i) = Di ^ C for each object i of 3. If it exists, the colimit (also direct or 
inductive limit) is an object lir^ of £ that is the vertex of a cone of canonical 
m,orphisms Ui: Di — > lim S) such that the following imivcrsal property holds: 
any co-cone (C, c) of & factors through it with a unique m,ediating morphism 
(also canonical morphism) uc, which makes the following diagram commute: 

Di^^^P^^^C . (308) 

For an illustrative example, consider the index category J consisting of only the 
objects and morphisms 1 ^ 3 — 2, while ^ is a functor from 3 to the category 
of sets. The diagram Si, the colimit and the co-cone of canonical morphisms fit 
into the following commutative diagram 

r>3 > D2 (309) 

£>i^^lii^^ . 

The universality of the limit identifies the set lii^^ with the quotient of the 
disjoint union Di U D2 by the equivalence relation that identifies elements di ^ 
d2, respectively of Di and D2, provided they are images of the same element 
^3 e -D3, di -fH ^3 I-)- ^2- If the object 3 and the morphisms from it were absent 
from the index category 3, we would simply have l^m ^ = Di U D2. In other 
categories, we may keep the same intuition, but replace the disjoint union U by 
the categorical coproduct. 
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